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VORWORT 



Die 10. Jahrestagung der Deutschen Gesellschaft für Unternehmensforschung 
zeichnete sich in verschiedenen Hinsichten besonders aus: Sie war nicht nur 
eine Jubiläumstagung, sondern sie war gleichzeitig die letzte Jahrestagung 
der DGU. Man könnte geneigt sein, diese Tatsache zu bedauern, wenn nicht 
gleichzeitig mit der Auflösung der DGU etwas geschehen wäre, was viele 
deutschem Untemehmensforscher seit geraumer Zeit gewünscht hatten: Die 
DGU schloss sich mit der Gesellschaft für Operations Research (früher 
AKOR) zu einer einheitlichen deutschen Gesellschaftler Deutschen Gesell- 
schaft für Operations Research (DGOR) zusammen. So ist das Ende der DGU 
gleichzeitig der Beginn der DGOR, einer Gesellschaftler man nur wünschen 
kann, dass sie in Fortsetzung der Arbeit der Deutschen Gesellschaft für Un- 
ternehmensforschung und der Gesellschaft für Operations Research der Fort- 
schritt auf dem Gebiete des Operations Research fördern möge. 

Die zehn Jahre des Bestehens der DGU bezeichnen in Deutschland die Ent- 
wicklungsphase des Operations Research oder der Unternehmensforschung, 
wie sie weitgehend in Deutschland genannt wurde, von den absoluten Anfängen 
bis zu dem Stadium in dem in weiten Kreisen der Wirtschaft und der Lehre 
die Nützlichkeit der Anwendung dieser Methoden nicht mehr bezweifelt wird. 
Dies kommt auch im Programm dieser Jubiläums tagung zum Ausdruck: 

Zum einen waren die Themenkreise noch nie so weit gespannt wie in diesem 
Jahr. Sie reichen von praktischen Anwendungen auf den verschiedensten Ge- 
bieten einschliesslich der Politik, des Umweltschutzes und der Familien- 
planung auf der einen Seite bis zu bedeutenden theoretischen Ergebnissen, 
von denen noch nicht vorauszusehen ist, auf welchen Wissensgebieten sie 
einmal angewandt werden. 

Zum anderen konnte bisher noch nie eine so hohe internationale Beteiligung 
verzeichnet werden wie bei dieser Tagung. Neben Untemehmensfors ehern 
aus Nachbarländern wie Belgien, den Niederlanden, Österreich und der 
Schweiz, konnten wir auch Kollegen aus Grossbritannien, den USA und Ka- 
nada begrüssen. 

Um sowohl den Teilnehmern an der Tagung als auch den Interessierten, die 
nicht an der Tagung teilnehmen konnten, Gelegenheit zu geben, sich mit den 
Ergebnissen der Tagung vertraut zu machen, wurden die Tagungsbeiträge 
in den vorliegenden "Proceedings” veröffentlicht. Dieser Brauch, von der 
DGU auf ihrer letzten Jahrestagung begonnen, soll ebenfalls von der DGOR 
weitergeführt werden. Möge er Praktikern und Forschem gleichennassen 
nützen. 



Prof. Dr. H. - J. Zimmermann 
(1. Vorsitzender der DGZ) 
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Unternehmensforschung in Wissenschaft, Wirtschaft und Politik 
von H. Künzi, Zürich 



1. Einleitung 



Die Unternehmensforschung, auch Operations Research genannt, 
hat ihre Wurzeln weitgehend in der angewandten Mathematik und in 
der Statistik. Man geht nicht fehl, wenn man die Theorie des neuen 
Forschungszweiges als Teilgebiet der angewandten Mathematik be- 
zeichnet. Aus dieser Ueberlegung heraus kann man die Unternehmens- 
forschung auch als diejenige Richtung der angewandten Mathematik be- 
trachten, die sich dem Anwendungsgebiet der Wirtschaft und der Po- 
litik zuwendet. Sicher gehen theoretische und praktische Ansätze aus 
dem Gebiete der Unternehmensforschung schon sehr weit zurück; 
aber das, was wir heute als eigentliche Unternehmensforschung be- 
zeichnen, weist auf eine kaum dreissigjährige Geschichte hin, deren 
Geburtsort eindeutig in den Vereinigten Staaten zu suchen ist. In Eu- 
ropa hat man sich dann rasch für die wertvollen Gedankengänge in- 
teressiert, vor allem die Engländer und Franzosen. In den ersten Jah- 
ren dieser Entwicklungsgeschichte beschäftigte man sich vorwiegend 
mit theoretischen Betrachtungen. Im Vordergrund stand vor allem die 
Theorie der linearen und nichtlinearen Optimierung. Rasch entstand 
auch die Spieltheorie und die umfassende Entscheidungstheorie. Das 
machtvolle Einsetzen des Computers und der Computerwissenschaften 
gaben der Unternehmensforschung gewaltige Impulse, so zum Beispiel 
in der Simulationstheorie. Was Deutschland und die Schweiz anbetrifft, 
muss heute festgestellt werden, dass dieser neue Zweig der ange- 
wandten Mathematik relativ spät Fuss fassen konnte. Man beschränkte 
sich sporadisch auf die Uebernahme einiger amerikanischer Modelle 
und fand dann sehr oft, dass sich diese nicht ohne weiteres für unse- 
re Verhältnisse eigneten. 
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Die eigentliche Geburtsstunde in Deutschland wie in der Schweiz 
kann im Jahre 1961 fixiert werden, als nämlich in beiden Ländern 
entsprechende Fachverbände ins Leben gerufen wurden. 

In Deutschland war es die Deutsche Gesellschaft für Unternehmens- 
forschung (DGU) und auf der südlichen Seite des Rheines die Schwei- 
zerische Vereinigung für Operations Research oder SVOR. Die bei- 
den Gesellschaften, die somit auf eine zehnjährige Tätigkeit zurück- 
blicken können, hatten zum Ziele, die Unternehmensforschung theore- 
tisch und praktisch zu fördern. Es muss als Glücksfall bezeichnet 
werden, dass in beiden Gesellschaften massgebende Persönlichkeiten 
aus der Wissenschaft und der Wirtschaft aktiv am Aufbau mitarbeite- 
ten. Obschon Skeptiker und Kritiker in beiden Ländern nicht fehlten, 
entwickelte sich das neue Forschungs- und Anwendungsgebiet rasch. 

Es entstanden Lehrstühle und Institute, und die fortschrittlichen Wirt- 
schaftsunternehmungen sowie politische Instanzen suchten Kontakte und 
Anknüpfungspunkte im noch bescheidenen Garten der Unternehmens- 
forschung. 

Heute darf mit berechtigtem Stolz festgestellt werden, dass 
sich aus den zarten Pflänzchen n Unternehmensforschung", die vor 
zehn Jahren in Deutschland und in der Schweiz gepflanzt wurden, an- 
sehnliche Bäume entwickelt haben, die nicht mehr aus dem Landschafts- 
bild wegzudenken sind. Diese Bäume weisen jetzt schon zahlreiche 
Aeste und Verzweigungen auf; trotzdem benötigen sie aber auch in 
Zukunft noch weitere Pflege. Man muss sie vor allem stärken, da- 
mit sie den verschiedenen Sturmwehen, die nicht selten im politischen 
und wirtschaftlichen Landschaftsgarten ihr Unwesen treiben, standhal- 
ten können. 
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2. Operations Research in Wirtschaft und Politik 

Nach zehn Jahren intensiver Arbeit in Theorie und Praxis ist 
festzuhalten, dass in beiden Ländern, wenn auch nicht in allen Spar- 
ten, so aber doch in einigen wichtigen Teilgebieten, der Anschluss 
an die weltweite Entwicklung gefunden werden konnte. Das trifft 

vor allem in den zentralen Gebieten der Optimierung, der Ent- 
scheidungstheorie und teilweise beim Gebiet der Simulation zu. 
Natürlich verfügen wir noch bei weitem nicht über den Forschungs- 
stand wie das z. B. in den USA der Fall ist; aber an internationalen 
Kongressen ist die Präsenz aus den beiden Ländern jeweils beacht- 
lich. Wirtschaftsbetriebe verfügen teilweise schon über eigene Opera- 
tions Research-Gruppen, und verschiedene Unternehmungen lassen 
sich durch private oder Universitätsinstitute beraten. Die frühere 
Skepsis ist weitgehend zusammengeschrumpft, und praktisch alle 
fortschrittlichen Betriebe, die sich mit entsprechenden Problemen 
au seinander zu setzen haben, erkennen den Wert der neuzeitlichen 
Methoden. Etwas zaghafter ist die Politik oder die öffentliche Ver- 
waltung in Richtung des Einsatzes der Unternehmensforschung vor- 
gegangen. Neben Aemtern, die sich gerne beraten lassen, gibt es 
noch immer recht konservative Stellen, die einer wissenschaftlichen 
Durchleuchtung politischer Aufgaben ablehnend gegenüberstehen. Wir 
alle wissen aber, dass zwischen einem umfangreichen wirtschaftli- 
chen Unternehmen und einer grossen öffentlichen Verwaltung weit- 
gehende Parallelen bestehen. Je länger je mehr wird die öffentli- 
che Verwaltung zum eigentlichen Unternehmen, ja zum Grossunter- 
nehmen, das nur noch funktionieren kann, wenn es von einem mo- 
dernen Management und von einer neuzeitlichen Betriebsführung ge- 
lenkt wird. Vielleicht bedarf es an solchen Stellen noch in vermehr- 
tem Masse der Orientierung und der Information. 

Unternehmensforschung gehört nicht nur in die wirtschaftliche 
Unternehmung, um diese zum Blühen zu bringen, sondern auch in 
die öffentliche Verwaltung, um diese zum besseren Funktionieren 
zu bringen. 
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3. Ausbildungsprobleme 



Wenn heute noch nicht überall von einem durchschlagenden 
Erfolg der Unter nehmensf or schung, vor allem im wirtschaftlichen 
Sektor, gesprochen werden kann, so ist dies sicher in hohem Masse 
darauf zurückzuführen, dass leider in Deutschland und in der Schweiz' 
noch recht wenig erstklassige Spezialisten vorhanden sind. Das be- 
zieht sich sowohl auf die Wissenschafter wie auf die entsprechenden 
Praktiker. Zahlreiche Projekte, die sich hervorragend für Opera- 
tions Research- Studien geeignet hätten, konnten aus diesem Grunde 
nicht behandelt werden. Natürlich brauchen wir in der Unternehmens- 
forschung nicht ausschliesslich Theoretiker, die die entsprechende 
mathematische oder statistische Seite beherrschen, sondern man ist 
in hohem Masse auf solche Spezialisten angewiesen, die einerseits 
die Theorie kennen, aber daneben in der Lage sind, auch das prak- 
tische wirtschaftliche oder volkswirtschaftliche Problem zu über- 
blicken, um darauf den theoretischen Apparat anzuwenden. Der eigent- 
liche Spezialist in der Unternehmensforschung muss daher vielseitig 
sein, theoretisch gründlich ausgebildet und in den verschiedenen prak- 
tischen Sparten bewandert. Leider fehlen uns in Deutschland und in 
der Schweiz noch weitgehend typische Ausbildungsstätten der Unter- 
nehmensforschung. Der Einzelne ist auf sich allein angewiesen; er 
muss sein Studium aus mehreren Fakultäten selber Zusammentragen, 
was sicher keinen allzugrossen Anreiz für diese Fachrichtung auszu- 
üben vermag. Meines Erachtens sollte in der Neugestaltung der Lehr- 
pläne an unseren Universitäten dieser offensichtlichen Lücke Rechnung 
getragen werden. Es muss untersucht werden, ob ein Studium der Un- 
ternehmensforschung, das sich auf Mathematik, Wirtschaftswissen- 
schaft und Technik abstützt, sich als sinnvoll erweist, oder ob die 
Ausbildung - ich persönlich weise eher in diese Richtung - in Form 
von ,! post graduate 11 Studien anzustreben sei. So oder so benötigen 




7 



wir sicher in Deutschland wie in der Schweiz wissenschaftliche 
Zentren, die im Gebiete der Unternehmensforschung eng miteinan- 
der verbunden sind. 



4. Erreichtes in der Unternehmensforschung 



Es würde den Rahmen dieses Referates bei weitem sprengen, 
wenn ich das bisher Erreichte im Bereiche der Unternehmensfor- 
schung in Deutschland und der Schweiz auch nur kurz skizzieren 
wollte. Vielmehr liegt mir daran, einige Schwerpunkte aufzuzeich- 
nen, mit denen ich durch die Aktivität meines früheren Institutes 
für Operations Research und elektronische Datenverarbeitung der 
Universität Zürich konfrontiert wurde. Ich beschränke mich auch , 
und dies wiederum im Zusammenhang mit meiner jetzigen Tätigkeit 
im öffentlichen Dienst, auf vorwiegend staatliche Probleme. 

A) Volkswirtschaftliche Probleme 

In Zusammenarbeit mit dem eidgenössischen Volkswirtschafts- 
departement haben wir in den letzten Jahren die modernsten mathe- 
matischen und technischen Hilfsmittel eingesetzt, um zentrale volks- 
wirtschaftliche Projekte, vorwiegend auf dem Gebiet der Kriegsvor- 
sorge, lösen zu helfen: an erster Stelle erwähne ich den sogenannten 
schweizerischen Anbauplan. Seine Aufgabe besteht darin, die Bebauung 
des knappen schweizerischen Bodens So zu gestalten, dass unsere 
Bevölkerung in Notzeiten aus der eigenen Produktion ernährt werden 
kann. In Analogie zum bekannten "Plan Wahlen" aus dem zweiten 
Weltkrieg wurde ein grosses Modell mit mehreren hundert Relationen 
und Variablen aufgestellt. Dieses umfangreiche Modell, das sich exak- 
ter mathematischer Methoden bedient, ist in der Lage, das mathe- 
matische Optimum mit Sicherheit anzugeben. Die Berechnung kann 
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nur auf einer sehr leistungsfähigen elektronischen Rechenanlage 
erfolgen. Die Rechenzeit beläuft sich auf ungefähr eine Stunde. 

Dabei muss natürlich erwähnt werden, dass die Ausfertigung des 
Computer Programms eine umfangreiche Arbeit darstellte. Doch 
dieses Programm kann, wenn es einmal vorhanden ist, immer 
wieder benutzt werden. 

Neben diesem Anbauplan werden im Rahmen der neugeschaf- 
fenen Sektion KOR (Kriegswirtschaftliches Operations Research) 
verschiedene weitere Projekte studiert, die hier nur kurz erwähnt 
werden sollen: 

a) Vorbereitung einer Lebensmittelrationierung im Kriegsfall unter 
Einsatz von Computern 

b) Probleme der Mehl- und Brotversorgung. Dabei gilt es zu er- 
mitteln, auf welche Weise die Brot Versorgung der Schweiz im 
Kriegsfall kostenoptimal gesichert werden kann. Fragen treten 
auf in der Form: Welche Ortsgetreidestelle hat welche Mühle 
zu beliefern, und zu welchem Verbraucher Zentrum wird das 
Mehl dann weitergeleitet, so dass der Mehlbedarf gedeckt wird 
und die Kosten minimal werden. 

c) Zu unseren interessanten Problemen gehört der Operations Re- 
search-Teil des zur Zeit im Aufträge stehenden Agrargutachtens 
des Bundesrates über die langfristige Struktur der schweizeri- 
schen Landwirtschaft. Hier steht die optimale Produktions Struk- 
tur für zwei Varianten im Vordergrund, nämlich 

1) bei einem allfälligen EWG- Anschluss unter spezieller Berück- 
sichtigung der EWG- Preise, 

2) bei keinem EWG- Anschluss. 

Auch bei dieser Aufgabe ist die Kriegsvorsorge zu berücksichti- 
gen, im Zusammenhang mit der minimalen Anbaubereitschaft. 

d) Die Nachfragefunktionen für die ca. 20 wichtigsten Agrarprodukte 
wurden mittels ökonometrischer Methoden geschätzt. Die Kenntnis 
solcher Nachfragefunktionen der wichtigsten Agrarprodukte kann 
der schweizerischen Agrarpolitik bestimmt wichtige Dienste- leisten. 
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e) Besonderes Interesse möchten wir der Aufstellung volkswirt- 
schaftlicher Gesamtmodelle widmen. Solche Modelle haben 
sich vor allem mit dem Problem der allgemeinen Interdepen- 
denz, d. h. mit der Tatsache, dass alle wirtschaftlichen Vor- 
gänge miteinander verkettet sind, au seinander zu setzen. 



B) Militärische Probleme 



a) Evaluation von Kampfflugzeugen 

Bei diesen Operations Research- Studien geht es im Wesent- 
lichen darum, aus den dem Eidgenössischen Militärdepartement 
angebotenen Flugzeugtypen denjenigen Typ zu ermitteln, der die 
verlangten militärischen Wirkungen mit den relativ geringsten Ko- 
sten erzielt. Zur Untersuchung eines sehr wichtigen Teilproblems, 
nämlich der Ueberlebenswahrscheinlichkeit eines Kampfflugzeuges 
im Kriegseinsatz, wurde ein Operations Research- Modell zur Si- 
mulierung von Luftkämpfen auf dem Elektronenrechner entwickelt. 
Für dieses Modell wurden die Einflussfaktoren auf die Bewegung, 
die Sichtung und auf den Luftkampf zweier Flugzeuge oder Flug- 
verbände untersucht. Durch das möglichst realistische Nachbilden 
(Simulieren) von Luftkämpfen mit dem Computer soll insbesondere 
die Bedeutung wichtiger technischer Eigenschaften, wie zum Bei- 
spiel das Beschleunigungsvermögen, für die Ueberlebenswahrschein- 
lichkeit eines Kampfflugzeuges abgeschätzt werden. Je grösser sei- 
ne Lebenserwartung im Kriegseinsatz ist, umso mehr Einsätze 
können mit demselben Flugzeug geflogen werden, bevor es durch 
Abschuss zerstört wird. Um zuverlässige Aussagen über die Ue- 
berlebenswahrscheinlichkeit en erhalten zu können, mussten auf ei- 
nem Computer gegen 100 000 Luftkämpfe simuliert werden. 

b) Weitere militärische Operations Research- Aufgaben 

In Kürze sei noch auf einige Studien von militärischen 
Operations Research- Arbeiten stichwortartig hingewiesen, die in 
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den letzten Jahren in der Schweiz behandelt wurden: Probleme 
der Munitionslagerung - Simulation von Bombardierungen eines 
Flugplatzes - optimaler Erneuerungs zyklus bei Motorfahrzeugen - 
Transportmodelle - Pertstudien - Panzerabwehr simulation - Mo- 
delle zur Tiefflieger er f as sung u. a. m. 



C) Verkehrsplanungen 

Die zukünftig^ Bewältigung des Strassenverkehrs stellt die 
kommunalen, kantonalen und eidgenössischen Behörden vor mannig- 
faltige Aufgaben, die nur in Zusammenarbeit von Vertretern aus 
verschiedenen Fachgebieten mit Aussicht auf Erfolg in Angriff ge- 
nommen werden können. Für solche umfangreiche Arbeiten stellt 
das Operations Research wiederum geeignete Verfahren und Modell- 
strukturen zur Bearbeitung zur Verfügung. Ebenso erweisen sich 
die modernen Hochlei stung s- Rechenautomaten als ein äusserst wert- 
volles Hilfsmittel, welches zur Bearbeitung eines sehr grossen Auf- 
gabenspektrums verwendet werden kann. 

In den grossen Städten erfordert der öffentliche Verkehr eine 
zunehmend gründlichere Verkehrsplanung. Es ist zum Beispiel bei 
der Entwicklung einer Verkehrs Steuerung das Ziel, den Verkehrs- 
ablauf durch geeignete Massnahmen nach bestimmten Gesichtspunk- 
ten zu optimieren (minimale Wartezeit, grösstes Verkehrs volumen, 
etc. ). Analytische Verfahren eignen sich nicht immer zur Behand- 
lung solch komplexer Probleme; deshalb benützt man oft grosse 
Elektronenrechner und simuliert mit Operations Research- Methoden 
den zu untersuchenden Prozess. Beim Aufstellen von Simulations- 
modellen zeigen sich zudem oft weitere, bis dahin noch nicht er- 
fasste Zusammenhänge. Aus den Simulationsergebnissen können 
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Schlüsse auf eine günstige Fahrpiangestaltung der Strassenbahnen 
gezogen werden; ebenso vermitteln die Simulationsergebnisse An- 
gaben über die Auswirkungen der untersuchten Ampel st euer Strate- 
gien auf die übrigen Verkehrsteilnehmer, sowie Hinweise auf all- 
fällige notwendige bauliche Veränderungen. 



D) Bauplanung 

Die Investitionen im Bauwesen sind in den letzten Jahren 
enorm angestiegen. Durch eine sorgfältige Auswahl von Standor- 
ten sowie durch zweckmässige Raumdimensionierungen und Raum- 
anordnungen kann der Nutzwert von Gebäuden wesentlich erhöht wer- 
den. Bezüglich der Baukosten ist zu bemerken, dass durch die Wahl 
geeigneter Bauverfahren, durch die Bestimmung optimaler Serien- 
grössen von Bauteilen und durch die Wahl einer optimalen Bauge- 
schwindigkeit wesentliche Kosteneinsparungen erreicht werden 
können. 

Bei den Berechnungen wird ein gegebenes Zahlenmaterial un- 
ter vielen Gesichtspunkten ausgewertet und mittels Verfahren des 
Operations Research optimiert. Die Auswahl der Rechenoperationen 
erreicht in der Regel einen Umfang, der nur durch den Einsatz 
elektronischer Grossrechenanlagen in wirtschaftlicher Weise bewältigt 
werden kann. 



Es lag mir daran, an diesen wenigen Beispielen den Einsatz 
der Unternehmensforschung in der Schweiz zu erläutern. 

Oft mutet es sonderbar an, dass diese Wirkung s voll en Metho- 
den nicht schon früher verwendet wurden. Ich glaube, dass unsere 
grossen Mathematiker des letzten Jahrhunderts ohne weiteres in dei 
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Lage gewesen wären, den erforderlichen theoretischen Apparat 
zu liefern. Der Grund dieses späten Einsatzes des Operations Re- 
search lässt sich aber ganz einfach erläutern: Ohne Computerein- 
satz lohnt es sich nicht. Operations Research zu betreiben. Es 
lässt sich eindeutig feststellen, dass die Geburtsstunde des Opera- 
tions Research ziemlich genau mit der Geburtsstunde des Com- 
puters zusammenfällt. 

Es soll unser Bestreben sein, diese beiden Komponenten, 
nämlich das Operations Research und den Computer, am richtigen 
Ort einzusetzen und auf die richtige Weise zu fördern. Die obigen 
Beispiele, die sich durch zahlreiche andere noch erweitern liessen, 
haben uns eindrücklich vor Augen geführt, wie nützlich diese beiden 
Komponenten sind. 

5. Der Ausbau der Unternehmensforschung auf nationaler Ebene 



In den vorangegangenen Abschnitten glaube ich deutlich ge- 
zeigt zu haben, wie sehr weite Kreise heutzutage auf die Unter- 
nehmensforschung angewiesen sind. Die Konsequenz daraus lautet: 
Förderung dieses Zweiges sowohl in wissenschaftlicher wie auch in 
praktischer Richtung auf breiter Ebene. 

Die beiden Vereinigungen tun das ihrige dazu. Neben der DGU 
und der SVOR bemühen sich noch weitere Institutionen in ver dankens- 
werter Weise für einen gezielten Ausbau der Unternehmensforschung. 
In hohem Masse sehe ich darin aber eine staatliche Aufgabe. Was 
dem Staate direkt und indirekt nützt, hat er entsprechend zu fördern. 

Unter Mitwirkung von Kollege Professor W. Kr eile habe ich 
ein Programm ausgearbeitet, das der* staatlichen Förderung der 
Wirtschaftswissenschaften im allgemeinen und der Unternehmensfor- 
schung im speziellen in Deutschland und der Schweiz Rechnung trägt. 
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Es kann wie folgt zusammengefasst werden: 

1. Sowohl in Deutschland wie auch in der Schweiz sollten je ein 
oder zwei Zentren für Unternehmensforschung an Universitäten 
oder Technischen Hochschulen errichtet werden. An diesen 
Zentren sollten diejenigen Unternehmensforscher vereinigt wer- 
den, die Theorie und Praxis von hoher Warte aus betreiben. Die 
Berufung an diese Zentren sollte durch ein unabhängiges Gremium, 
das auch die hervorragendsten ausländischen Fachvertreter um- 
fasst, auf Zeit erfolgen (ca. fünf Jahre). Nach dieser Zeit tritt 
der Berufene wieder in seine alte Stellung zurück. 

Eine genügende Zahl von Stellen für Assistenten und Schreib- 
kräfte müssen für das Zentrum zur Verfügung stehen. 

2. Es sollten alle möglichen Massnahmen ergriffen werden, um 
eine Verbindung von reiner Forschung und Praxis herzustellen. 
Hierzu wird vorgeschlagen: 

a) Regelmässige Abordnung von leitenden Beamten in der Ver- 
waltung und Persönlichkeiten in Führungspositionen der Wirt- 
schaft zu Kursen an der Universität, in denen sie mit dem 
Fortgang der Wissenschaft vertraut gemacht werden und 
ihrerseits die Forschung auf dringende Probleme der Praxis 
hinweisen. 

b) Zeitweise Abordnung von Wissenschaftern in die Verwaltung, 
um bestimmte Planung sauf gaben zu leiten, für die sie ihrer 
Forschungsrichtung nach kompetent sind. 

3. Der Staat sollte die Unternehmensforschung in erheblichem Masse 
fördern. Besonderes Gewicht ist auf die Grundlagenforschung zu 
legen. Die angewandte Forschung könnte durch direkte Regierungs- 
aufträge für bestimmte Projekte, ähnlich wie in den USA, ge- 
fördert werden. 
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Die wirtschaftliche Zukunft beider Länder hängt vom Stand 
der Wissenschaft und der Technik ab. Die wissenschaftliche 
Vorherrschaft der USA ist weitgehend durch die grössere 
Förderung der Forschung durch die Regierung bedingt. Unsere 
Regierungen sollten hier entsprechend gleichziehen, damit wir 
nicht weiter zurückfallen. 

4. Sowohl in Deutschland wie in der Schweiz sollte die Studien- 
reform für das Gebiet der Unternehmensforschung vorange- 
trieben werden. 

5. Die Zentren in Deutschland und in der Schweiz sollten eng 
Zusammenarbeiten, um dadurch eine noch grössere Aus- 
strahlungskraft zu gewinnen. 

Wiederum erachte ich die DGU und die SVOR als Institutionen, 
die wesentlich zur Realisation dieses Programms beitragen können, 
das sicher dafür Gewähr bieten würde, dass die Unternehmens- 
forschung in unseren beiden Ländern die ihr zustehende Bedeutung 
erlangen wird. 
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Beziehungen zwischen Mathematik und Wirklichkeit 
in mathematischen Modellen von empirischen Strukturen 
von A. Jaeger, Bochum 



1. EINLEITUNG 

Der amerikanische Mathematiker G. B. DA.NTZ IG bemerkte einmal in 
einem Interview mit einem Journalisten, man könne mit einem Modell 
aus dem von ihm entwickelten Gebiet der linearen Programmi erung 
streng beweisen, daß von allen möglichen Arten des Zusammenlebens 
der Geschlechter in einer sozialen Gruppe (Monogamie, Bigamie, Po- 
lygamie usw.) die Monogamie optimal sei. Darauf soll der Journalist 
nur kopfschüttelnd mit dem Kommentar reagiert haben: Dann müsse 
DANTZIG sicherlich mit der falschen Art von Modellen gearbeitet 
haben ([2], S. 367). 

Diese köstliche Geschichte ereignete sich schon vor mehr als fünf- 
zehn Jahren. Heutzutage sind weitaus vielfältigere Mißverständnisse 
bei Verwendung der Vokabel Modell denkbar, da sich diese zu einem 
allseits beliebten Modewort mit wissenschaftlichem Nimbus entwik- 
kelt hat. Die umgangssprachliche Abgegri f fenhei t dieser Vokabel 
kann nun leicht zu dem Fehlschluß führen, daß sich auch der Begriff 
des abstrakten Modells im Rahmen einer Wissenschaft eigentlich nur 

in eine Metapher oder in ein vages Schlagwort auflöst. Mit Hilfe 
weniger elementarer mathematischer Grundbegriffe läßt sich jedoch 

der wissenschaftliche Mod^l 1 begri f f präzisieren. Der wesentliche 
Hinderungsgrund für das Verstehen der exakten Definition des mathe- 
matischen Modelles besteht für den Nicht-Spezialisten lediglich da- 
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rin, daß das Stichwort mathematisch bei ihm eine sofortige abwei- 
sende Reaktion auslösen kann. Denn in der Schule mag das Fach, wel- 
ches unter der Bezeichnung Mathematik gelehrt wird, meist nur 
Schrecken und ein Gefühl der Ohnmacht vermittelt haben, und viele 
sogenannte "Gebildete" wollen mathematische Überlegungen gar nicht 
anstellen, ja geben sogar ein Fehlen von mathematischen Kenntnissen 
ohne Scham und fast mit Stolz zu, während sie wohl gleichzeitig 
Lücken auf schöngeistigen Gebieten sorgfältig zu verbergen suchen. 

Es kann hier nicht der Ort sein, im einzelnen zu analysieren, wo- 
rauf diese bedauerliche Entwicklung gerade in der heutigen, so sehr 
durch die exakten Naturwissenschaften geprägten Zeit zurückzuführen 
ist; aber vielleicht muß man ungenügender Relevanz des Mathematik- 
unterrichtes zur realen Welt einen erheblichen Anteil der Verant- 
wortung zuschreiben. Dabei liegt, wie in diesem Referat behutsam 
entwickelt werden soll, in dem Verstehen des Zusammenhanges zwi- 
schen Wirklichkeit und Mathematik gerade der Schlüssel zum Verständ- 
nis des Model 1 begri ff s der Wissenschaft. 

2. ABSTRAKTIONEN UND KLASSEN 

Wie die Entwicklung der menschlichen Sprache, mit Ihren akustischen 
und später auch optischen Symbolen, durch eine Folge von Verfeine- 
rungen von Abstraktionen aus Gegebenheiten der realen Welt erfolgt, 
so geschieht dies auch mit der Entwicklung jeder Wissenschaft, ein- 
schließlich der Mathematik. Abstraktion bedeutet hierbei das Weg- 
lassen von Eigenschaften, die Trennungsmerkmale beinhalten, so daß 
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die dann verbleibenden Eigenschaften Gemeinsamkeiten signalisieren. 
Unter Verwendung des Begriffes der Menge , dem heutzutage ja 
schon immer mehr die Schulkinder (und deren Eltern) ausgesetzt 
sind, wenn sie die sogenannte Neue Mathematik lernen müssen, läßt 
sich jeder Abstraktionsvorgang dahingehend interpretieren, daß je- 
dem Element einer Menge von Gegebenheiten der realen Welt oder von 
Begriffen genau ein Element einer neuen Menge von abstrakteren Be- 
griffen zugeordnet wird. Durch diese Abbildung, wie der 
Mathematiker gewöhnlich eine derartige Zuordnung nennt, wird im 
allgemeinen verschiedenen Elementen der Ausgangsmenge ein und das- 
selbe Bild der neuen Menge zugeteilt. So kann man etwa, aus- 
gehend von der Menge der verschiedenen Baumarten, so z.B. von den 
Begriffen Tanne, Fichte, Kiefer, zu dem Oberbegriff Nadelbaum, aber 
z.B. von den Begriffen Eiche, Buche, Birke, zu dem Oberbegriff 
Laubbaum und danach von den Begriffen Nadelbaum und Laubbaum zu dem 
Oberbegriff Baum übergehen (siehe Fig. 1). 




Fig. 1 .Abstraktionsschritte in der Botanik 



Jeder solche Abstraktionsschritt führt zu einer Vergröberung der 
Beschreibung und Auffassung der betrachteten Objekte, z.B. Gegeben- 
heiten der realen Welt, bei welcher das ursprüngliche Gleichsein 
dieser Objekte zu einem Xhnlichsein - ln bezug auf die Ubriggeblie- 
benen Gemeinsamkeiten, also hier im Beispiel die Gattungseigenschaf- 
ten von Bäumen - vergröbert wird. Oie Elemente der Ausgangsmenge 
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werden hierdurch in Klassen aufgespal ten , wobei jede 
Klasse genau aus all denjenigen Elementen besteht, denen dasselbe 
Bild, also dasselbe abstraktere Objekt, zugeordnet ist. 

3. AUSSAGEFORMEN UND RELATIONEN 

Der menschliche Geist konstruiert nun aber Abstraktionen nicht - 
oder wenigstens vorwiegend nicht nur - um ihrer selbst willen, er 
bedient sich ihrer auch, um Eigenschaften, Beziehungen, Gliederun- 
gen, Gesetzmäßigkeiten, Zusammenhänge, kurz Organisationsmuster 
oder -merkmale, zu finden oder zu präzisieren. Eine solche Bezie- 
hung trat soeben schon auf, die Beziehung der Zugehörigkeit von 
zwei Elementen einer Menge zu einer Klasse im obigen Sinne. Derar- 
tige Beziehungen lassen sich mit Hilfe der Sprache (der Umgangs- 
sprache oder einer wissenschaftlichen Fachsprache mit ihren Fach- 
symbolen, einschließlich der mathematischen Sprache, je nach den 
zur Betrachtung stehenden Objekten) folgendermaßen beschreiben: 

Man stelle eine Behauptung auf, die grammatikalisch wie ein Satz 
aussieht, aber etwa ln derselben Welse variable Größen enthält wie 
z.B. die sogenannten “Unbekannten" x,y, ... bei mathematischen Glei- 
chungen, oder wie die "großen Unbekannten" X,Y bei den Steckbriefen 
von Verbrechern. Genauer gesagt, derartige Behauptungen sollen ein 
Symbol oder mehrere Symbole als grammatikalischen Satzteil oder 
als Satzteile enthalten, von denen jedes jeweils stellvertretend 
fUr ein beliebiges Element einer gleichzeitig spezifizierten Menge 
steht. Unabhängig davon, für welchen Wahrheitsbegriff man sich ent- 
schieden hat (ob, grob ausgedrückt: wahr z.B. logisch beweisbar 
oder experimentell verlf izlerbar bedeuten soll), wird die Frage 
nach der Wahrheit einer solchen mit Variablen behafteten Aussage, 
die A u s s a g e f o r m ([5], S. 40) oder auch Satz- 
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form, Aussagefunktion oder S a t z - 
funktion genannt wird, im allgemeinen erst dann einen Sinn 
haben, wenn man die Variablen durch fest herausgegriffene Elemente 
der zugrundegelegten Mengen ersetzt, "durch Konstan- 
te spezialisiert" (siehe [4], S. 12). Steht z.B. 
x stel 1 vertretend für eine Zahl (genauer: für ein Element der Men- 
ge aller Zahlen), so ist die Aussageform "x ungleich 3" falsch für 
die Spezialisierung x = 3 (d.h., wenn die Variable x durch die 
Konstante 3 ersetzt wird), aber richtig für alle anderen Speziali- 
sierungen. Oder steht X stellvertretend für eine Person (genauer: 
für ein Element aus der Menge aller Menschen) und Y für eine Be- 
schäftigung (genauer: für ein Element aus der Menge aller Beschäf- 
tigungen), so ist "X übt zur Zeit die Beschäftigung Y aus" eine 
Aussageform mit zwei Variablen. Beispiele anderer Aussageformen, 
zu denen sich jeweilige Bezugsmengen leicht konstruieren lassen: 

" x + y * z 4 " , 

" X hat an der Universität U am Tage T die Prüfung P 
bestanden " , 

" X ist mit Y verheiratet " , 

" X ist Vater von Y " , 

H Der Beobachtungswert x ist gleich dem Beobachtungs- 
wert y gemessen worden "(Meßgleichheit, 
d.h. Gleichheit im Rahmen der Meßgenauigkeit). 

" Der Computer hat ausgerechnet, daß x denselben Wert 
hat wie y" (Rechengleichheit, d.h. 
Gleichheit im Rahmen der Rechengenauigkeit). 

" X ist Oberbegriff von Y " . 

Durch jede derartige Aussageform mit einer Bezugsmenge oder bei meh- 
reren Variablen auch mit mehreren (nicht notwendigerweise verschie- 
denen) Bezugsmengen 1st nunmehr eine Beziehung zwischen den Elemen- 
ten dieser Mengen gewonnen worden, die n - s t e 1 1 1 g e 




20 



Relation in oder auf der Bezugsmen- 

g e bzw. zwischen den Bezugsmengen genannt wird, wenn 

n die Anzahl der darin auftretenden verschiedenen Variablen ist 
( [1], S. 256, [4], S. 101 u. S. 109). Für den Spezialfall einer 
einstelligen Relation ist auch die Bezeichnungsweise Eigen- 
schaft üblich. Man sagt, daß (bei n = 1) ein Element s aus 
der Bezugsmenge bzw. (bei n > 1) eine Folge (z.B. ein Paar, ein Tri- 
pel, ein Quadrupel usw.) Sj, s 2> ..., s p von Elementen aus den je- 
weiligen Bezugsmengen die Relation erfüllt oder 
zulässig für die Relation ist, wenn bei 
der Spezialisierung: x durch s ersetzt ("x = s"), bzw. den Spezia- 
lisierungen: Xj durch s 1 ersetzt, x 2 durch s 2 ersetzt, ..., x n durch 
s n ersetzt ("Xj * Sj, x 2 « s 2 , ...» x p * s n ") die Aussageform in 
eine Aussage umgewandelt wird, welche wahr ist. Wegen dieses Zusam- 
menhanges kann man eine Aussageform auch eine Relations- 
behauptung für die Bezugsmenge(n) nennen ( [4], S. 146). 
(Genauer gesagt: der Mathematiker definiert die Relation einfach 
als Zusammenfassung derjenigen Elemente, bzw. Elementepaare bzw. 
Elementefolgen, für welche die zugrunde liegende Behauptung wahr ist, 
also als die Menge der zulässigen Elemente, bzw. Elementepaare bzw. 
Elementefolgen . ) 

4. STRUKTUREN UND STRUKTURKOMPLEXE 

Betrachtet man eine Menge M zusammen mit einer Relation R oder mit 
mehreren Relationen R 1# R 2 , ..., die durch Behauptungen mit Varia- 
blen aus dieser Menge erklärt sind, so spricht man von einer struk- 
turierten Menge, kürzer von einer Struktur, und schreibt 
abkürzend (M,R) bzw. (M; Rj, R 2 , ...). (Man findet auch andere Be- 
zeichungsweisen, wie z.B. Relationengebilde, 
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relationales System, Relativ, Sy- 
stem, Ganzheit, Gestalt oder in speziellen 
Fällen Graph, Graphoid. Siehe z.B. [1] , S. 209 
und 284, [4], S. 119, [5], S. 60) Manchmal muß man sogar mehrere 
Mengen Mj, » ...» zusammen mit mehreren Relationen, die sich je- 
weils auf einzelne oder alle dieser Mengen beziehen, zugleich be- 
trachten, z.B. wenn man Wertordnungen, Bewertungen (siehe [4], S. 17) 
für Dinge, Phänomene, Ergebnisse, Alternativen usw. einführen will; 
in einem solchen Falle schreibt man natürlich abkürzend (Mj, M 2 » 

Rj> R 2 * • ••)• Auch dann soll die Zusammenfassung dieser Mengen und 
der mit ihnen betrachteten Relationen Struktur oder viel- 
leicht besser Strukturkomplex (wenn man ausdrück- 
lich betonen will, daß mehrere Mengen betrachtet werden) genannt 
werden. Für die zugrundeliegenden Mengen benutzt man Vokabeln wie 
Grundmenge(n) , Basismenge(n) oder 
Strukturträger. Das klingt vielleicht bisher alles 
nicht übermäßig mathematisch, denn es wurde bewußt der mathemati- 
sche Charakter dieser Überlegungen selbst auf die Gefahr verbaler 
Unschärfe hin heruntergespielt; aber es sollte mit Nachdruck darauf 
hingewiesen werden, daß die moderne Mathematik nicht nur solche 
Strukturen (in diesem Sinne), die auf "Zahlen und Figuren" aufge- 
baut oder von "Zahlen und Figuren" abgeleitet sind, untersucht, son- 
dern deradezu als die Theorie abstrakter Strukturen an sich angese- 
hen werden kann. Vielleicht sollte man sogar die oft emotional -ne- 
gative Vokabel Mathematik durch eine neue Wortbildung wie etwa Struk- 
turik ersetzen. 
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5. MORPH I SMEN VON STRUKTUREN 

Häufig beobachtet man nun zwei Strukturen gleichzeitig, vor allem 
wenn sie eine gewisse Art von struktureller Ähnlichkeit aufweisen, 
die zunächst für den Fall einer einzigen Grundmenge und einer ein- 
zigen Relation bei beiden Strukturen vorgestellt werden soll: Es 
läßt sich eine Abbildung von der zugrundeliegenden Menge M der ei- 
nen Struktur, Ausgangs- oder Original Struk- 
tur genannt, in die zugrundeliegende Menge N der anderen Struk- 
tur, Bildstruktur genannt, finden, welche die beson- 
dere Eigenschaft aufweist, daß aus der Gültigkeit der Relation R 
der ersten Struktur, der Ausgangs- oder Origi- 
nalrelation, für gewisse Elemente von M automatisch die 
Gültigkeit der Relation S der zweiten Struktur, der B i 1 d re- 
lation, für die entsprechenden Bildelemente in N folgt. Dazu 
ist es allerdings unbedingt erforderlich, daß die zwei zu verglei- 
chenden Relationen dieselbe Stellenzahl haben, d.h. 
daß S auch n-stellig sein muß, wenn R n-stellig ist. Wenn in diesem 
Sinne die Bildrelation gewissermaßen die Ori gi nal rel ation widerspie- 
gelt, spricht man von einer kompatiblen Struk- 
turabbildung oder von einem (Struktur-) 
Morphismus oder Homomorphi smus von (M,R) 
ln (N,S) (siehe z.B. [l], S. 288, [4], S. 123). 

Bedeutet z.B. R die Relation der Gleichheit, definiert durch die 
Aussageforn "x 1st gleich y", ln der Menge M, und bedeutet S die Re- 
lation der Gleichheit, aber jetzt für die Menge N definiert, so 1st 
ganz offensichtlich jede beliebig* Abbildung von M ln N ein Morphls- 
■us von (M,R) ln (N,S). Bedeutet dagegen S die Relation der Ungleich- 
heit für N, definiert durch die Aussageforn "x 1st nicht gleich y". 
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so ist nun plötzlich keine Abbildung von (M,R) in (N,S) ein Mor- 
phismus . 

6. ISOMORPHISMEN VON STRUKTUREN 

Bei Abbildungen von Mengen kann man in vielen Fällen durch Änderung 
der Zuordnungsrichtung wieder zu einer Abbildung, der inver- 
sen oder umgekehrten Abbildung, gelangen, nämlich 
genau dann, wenn in der zweiten Menge N jedes Element ein Bild ei- 
nes Elementes der ersten Menge, aber niemals Bild mehrerer Elemente 
der ersten Menge M ist. Bei dieser Umkehrung vertauschen sich dann 
gewissermaßen die Rollen von Bild und Original. Eine solche spezi- 
elle Abbildung nennt der Mathematiker eindeutig, bi- 
j e k t i v oder einfach umkehrbar. Typisches Beispiel 
ist die Partnerzuordnung in der Tanzschule, wenn kein Mauerblümchen 
übrig bleibt. Ist nun eine Abbildung M in N, die einen Morphismus 
von (M,R) in ( N #S ) ergibt, umkehrbar, so ist die Frage interessant, 
ob auch die umgekehrte Abbildung einen Morphismus bildet, und zwar 
natürlich dann von (N,S) in (M,R). In einem solchen speziellen Falle 
hat man es gewissermaßen mit struktureller Gleichheit zu tun und 
spricht von einem Isomorphismus der Ausgangsstruktur 
auf die Bildstruktur (siehe z.B. [1], S. 290, [4], S. 123, [5], 
S. 129). 

Ein typisches, wohlbekanntes Beispiel für einen Isomorphismus läßt 
sich mit Hilfe des Oberganges von den positiven Zahlen zu ihren Log- 
arithmen konstruieren; er bildet den Hintergrund für die Benutzung 
des Rechenschiebers oder der Logarithmentafel. Genauer gesagt, geht 
es hierbei um folgendes (siehe Fig. 2): Man nehme für M die Menge al 
ler positiven Zahlen, für R die dreistellige Relation, welche durch 
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die Relationsbehauptung "x . y * z" definiert ist, für N die Menge 
aller (reellen) Zahlen, für S die dreistellige Relation, welche 
durch die Relationsbehauptung "u + v * w" definiert ist, sowie die- 
jenige Zuordnung, welche jeder Zahl a ihren Logarithmus: log a zu- 
ordnet. Wenn nun für drei positive Zahlen a^ • a 2 * a 3 gilt, so 
gilt ;ür ihre Logarithmen: log a 1 + log a 2 ■ log a 3 , und umgekehrt. 
Also ist diese Zuordnung ein Strukturisomorphismus von der Struk- 
tur (M,R) auf die Struktur (N,S). (Dieses Beispiel zeigt übrigens 
auch gleichzeitig, auf welche Weise sich die elementaren Rechen- 
operationen der Schule als dreistellige Relationen auffassen lassen.) 




Fig. 2 Das Rechnen mit Logarithmen 



7. MORPHISMEN VON STRUKTURKOMPLEXEN 

Diese Begriffsbildungen lassen sich nun in naheliegender Weise auf 
Strukturkomplexe übertragen, wobei die verbale bzw. formelmäßige 
Beschreibung eigentlich komplizierter 1st als der zugrundeliegende 
Gedanke. Um die Unübersichtlichkeit nicht allzu groß werden zu las- 
sen, soll allerdings zunächst angenommen werden, daß die zu betrach- 
tenden Strukturkomplexe jeweils nur eine zugrundeliegende Menjo ha- 
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ben. Man gehe aus von einem Strukturkomplex (A; P^, ••• P p)» 

welcher sich aus der Menge A und einer Anzahl von Relationen P^, 

. .., Pp in dieser Menge zusammensetzt, und bringe diesen Struktur- 
komplex mit einem zweiten (i.a. anderen) Strukturkomplex (B; Q^, 

..., Q„) in Verbindung, welcher sich seinerseits aus der Menge B* 
und den Relationen Q^, in dieser Menge zusammensetzt. Dabei 

braucht keinesfalls verlangt zu werden, daß die Anzahl der Relatio- 
nen des ersten Strukturkomplexes gleich der Anzahl der Relationen 
des zweiten Strukturkomplexes ist, also unbedingt p * q erfüllt sein 
muß. Es soll nun von einer Abbildung des Ausgangskomplexes in den 
Bildkomplex gesprochen werden, wenn erstens eine Abbildung der Men- 
ge A in die Menge B festgelegt ist und zweitens jeder Relation der 
Ausgangsstruktur eine Relation der Bildstruktur mit gleicher Stel- 
lenzahl zugeordnet ist. Dabei ist keineswegs vorausgesetzt, daß zwei 
verschiedenen Ori di nal rel ati onen unbedingt zwei verschiedene Bild- 
relationen zugeordnet sein müssen. Nun besteht wiederum der Wunsch, 
die verschiedenen Zuordnungen so vorzunehmen, daß sich die Bild- 
struktur irgendwie ähnlich wie die Ori gi nal Struktur verhält; genauer 
gesagt, es ist nun das Folgende erwünscht: Erfüllen Ori gi nal el emente 
eine Ori gi nal rel ati on , so sollten die entsprechenden Bildelemente 
die entsprechende Bildrelation erfüllen. Trifft das allgemein zu, so 
wollen wir in wichtiger Verallgemeinerung der üblichen Definition 
weiterhin von einem Morphismus oder von Kompatibilität sprechen. 

BEISPIEL: A sei eine Menge verschiedener (physikalischer, chemischer, 
ökonomischer, soziologischer oder anderer) Zustände. Durch Meß- oder 
Beobachtungsanordnungen werden numerische Daten bestimmt, wobei an- 
genommen werden möge, daß verschiedene Zustände verschiedene Meßda- 
ten ergeben. Es seien die Relation P^ durch die Relationsbehauptung 
"x und y ergeben denselben Meßwert", die Relation P 2 durch die Re- 




26 



1 ati onsbehauptung "x und y sind gleich“ definiert. Die Menge B be- 
stehe aus den obigen Meßdaten, Q sei die Relation der Gleichheit 
in B. Unter diesen Voraussetzungen ergibt sich ein Strukturmorphi s- 
mus von (A; P^, P 2 ) in (B;Q) in dem von uns definierten Sinne, wenn 
jedem Zustand aus A sein Meßwert in B und beiden Relationen Pj, ? z 
die Relation Q zugeordnet wird (siehe Fig. 3). 




Fig. 3 




Ein Strukturmorphi smus 



Das folgende Beispiel ist sehr ähnlich. A sei die Menge aller reel- 
len Zahlen. Ein Computer soll Rechenoperationen ausführen, wozu er 
Rundungsregeln benutzt. Es seien nunmehr die Relation P 1 durch die 
Relationsbehauptung “x und y ergeben dieselbe Zahl bei der Abrundung 
durch den Computer", P 2 durch die Relationsbehauptung "x und y sind 
gleich" definiert, frie Menge B bestehe aus allen für die Zwecke die- 
ses Computers abgerundeten reellen Zahlen. Q sei die Relation der 
Gleichheit in dieser Menge B. Wiederum ergibt sich ein Strukturmor- 
phismus von (A; Pj, P 2 ) in (B;Q) in unserem Sinne auf die folgende 
Weise: Jeder reellen Zahl aus A wird die zugehörige abgerundete Zahl 
aus B und beiden Relationen P^, P 2 die Relation Q zugeordnet. 

Wir können noch allgemeinere Typen von Strukturmorphlsmen definieren. 
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wenn man zwei Strukturkomplexe (...» M, . ...» R, ...) und 
(..., N, ...; .... S, ...) mit jeweils mehreren zugrundegelegten 
Trägern und jeweils mehreren betrachteten Relationen auf die fol- 
gende Weise in Beziehung setzt: Als erstes ordne man jeder Menge M 
des ersten Komplexes, also jeder Originalmenge, als Ganzes eine der 
Mengen N des zweiten Komplexes, also eine Bildmenge, zu, wobei eben- 
falls nicht gefordert wird, daß auf diese Weise zwei verschiedene 
Originalmengen unbedingt mit zwei verschiedenen Bildmengen in Ver- 
bindung gebracht werden. Zweitens führe man eine Abbildung von je- 
der Originalmenge in die ihr im ersten Schritt zugeteilte Bildmenge 
ein. Drittens ordne man jeder Relation R der Ori gi nal Struktur eine 
Relation Q der Bildstruktur derart zu, daß zwei Bedingungen erfüllt 
sind: Einmal muß wieder die Stellenzahl beider Relationen gleich 
sein; und zum anderen muß, wenn sich eine Variable in der Relations- 
behauptung für eine Relation R auf eine bestimmte der Originalmen- 
gen M bezieht, die entsprechende Variable in der Relationsbehauptung 
für die diesem R* zugeordnete Relation S sich auf die diesem speziel- 
len M zugeordnete spezielle Bildmenge N beziehen. Die Begriffe Ori- 
ginal- und Bildelemente sowie Original- und Bildrelationen werden 
entsprechend wie früher definiert. Folgt nunmehr für alle Relationen, 
die zwischen den Elementen des ersten Strukturkomplexes definiert 
sind, aus ihrer Gültigkeit für Ori gi nal el emente stets zwangsläufig 
die Gültigkeit der entsprechenden Bildrelation für die entsprechen- 
den Bildelemente, so wollen wir auch in dieser sehr allgemeinen Si- 
tuation noch von einem Morphismus sprechen. 

8. DIE STRUKTUR EINES AUSSCHNITTES DER WELT 

Auch wer diese letzten Gedankengänge nicht ganz nachvollziehen mag, 
well in Ihnen zu viele Zuordnungen gleichzeitig betrachtet werden. 
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wird sicherlich erahnen, daß Strukturmorphi smen geeignet sind, 
Denkvorgänge zu beschreiben oder zu approximieren, welche vorgenom- 
men werden, wenn man Organisationsmuster in einem Ausschnitt der 
realen Welt aufsucht. 

Zunächst möge unter der Struktur eines Ausschnittes der Welt oder 
auch eines gedanklichen Gebildes (also z.B. der Struktur des Gehirns, 
der Struktur eines Unternehmens, der Struktur einer Sprache, der 
Struktur einer Wissenschaft usw.) die Gesamtheit aller Gesetzmäs- 
sigkeiten dieses Ausschnittes der Welt oder dieses abstrakten Ge- 
bildes verstanden werden, d.h. etwa eine Zusammenfassung der Form 
(Mj, M 2 , ...; R^, R 2 > •••)» welche sich aus den Gesamtheiten der 
Mengen , die zu diesem Ausschnitt gehören, und der Gesamtheit al- 
ler Beziehungen R. zwischen Elementen dieser Mengen zusammensetzt. 

J 

Allerdings muß man hierbei erwarten, daß - außer in ganz leicht zu 
übersehenden Fällen - noch erheblich allgemeinere Arten von Bezie- 
hungen als die hier bisher skizzierten Relationen benötigt werden, 
z.B. Begriffsbildungen, wie sie die Mathematik in topologischen oder 
wahrscheinlichkeitstheoretischen Strukturen definiert, doch vermut- 
lich auch noch neu zu erfindende Begriffsbildungen, die der heutige 
Mathematiker noch nicht kennt und die speziellen Strukturuntersu- 
chungen bei Forschungen in irgendwelchen Wissensgebieten entsprin- 
gen können. (Man kann tatsächlich vermuten, daß durch Entdeckung 
neuartiger komplizierter Gesetmäßigkeiten , z.B. In den biologischen 
Wissenschaften, sehr fruchtbare Impulse fUr die Untersuchung der 
dadurch angeregten mathematischen Begriffsbildungen ausgehen, d.h. 
z.B. neuartige mathematische Gebiete entstehen können.) 
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9. MATHEMATISCHE MODELLE YON EMPIRISCHEN STRUKTUREN 

Daß man nun l.a. nicht sämtliche Gesetzmäßigkeiten eines Ausschnit- 
tes der Welt auffinden und sprachlich beschreiben kann, dürfte recht 
einleuchtend sein: Allein die exakte Beschreibung der zugrundelie- 
genden Mengen - ganz zu schweigen von den Beziehungen zwischen die- 
sen Mengen - kann faktisch und semantisch schwierig und unmöglich 
sein. Ernste Streitfragen mögen entstehen, ob eine solche Beschrei- 
bung vollständig ist oder sogar, ob diese erstrebte Vollständigkeit 
überhaupt prinzipiell erreicht werden kann. Außerdem sind Struktur- 
komplexe häufig dynamischer Natur, d.h. sie verändern sich, und so 
1st es oft zwingend erforderlich, den Zeitpunkt oder Zeitraum genau 
festzulegen, auf den man die Beschreibung festlegen will. Man wird 
sich oft auf nicht ganz klar abgegrenzte, offene Bezugsmengen sowie 
auf einige der von dem Sachverhalt her möglichen Relationen (oder 
allgemeineren Arten von Beziehungen), auf eine offene 
Tei 1 Struktur , beschränken müssen. Konstruiert man nun 
unter diesen erschwerten Bedingungen eine Strukturabbildung eines 
Ausschnittes der empirischen Welt in eine mathematische Struktur, so 
wird man auch dann nach einem Morphismus in der hier definierten 
Form bezüglich sämtlicher ausgewählten und in Verbindung gebrachten 
Relationen (oder allgemeineren Beziehungen) streben. Falls man ei- 
nen solchen Morphismus inständig erhofft oder dieser mehr oder weni- 
ger genau vorliegt, nennt man die Strukturabbildung ein m a t h e 
matisches Modell der empirischen Struktur ([4], 

S. 125). ( Häufig wird auch nur die mathematieche Bildstruktur als 
Modell bezeichnet , wobei dann die Originalmengen und Ori gi nal rel a- 
tionen und die verbindenden Zuordnungen stillschweigend als bekannt 
vorausgesetzt werden.) Derartige Modelle können zu deskriptiven oder 
präskriptiven Zwecken aufgestellt werden, d.h. um die empirische 
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Struktur zu beschreiben und analysieren zu können oder um sie be- 
einflussen und regulieren zu können. 

Es kann nun nicht genug betont werden, daß die moderne Mathematik 
weit mehr Hilfsmittel als das Rechnen mit Zahlen und Zahlenfolgen 
und das Hantieren mit geometrischen Figuren zur Hand hat, um Pro- 
bleme in solchen mathematischen Bildstrukturen zu lösen, die Fra- 
gestellungen in der empirischen Struktur entsprechen. Z.B. lassen 
sich die eng miteinander verknüpften Methoden der Relationentheorie, 
mathematischen Logik, Graphentheorie und Booleschen Matrizentheorie 
in immer vielfältigerer Weise auf Probleme anwenden, welche Frage- 
stellungen aus solchen empirischen Gebieten widerspiegeln, bei de- 
nen quantitative Begriffe noch relativ wenig entwickelt und unter- 
sucht sind (z.B. bei einigen Gesellschaftswissenschaften, siehe [4]). 

Die Anwendungsmöglichkeiten von mathematischen Modellen sind äußerst 
vielseitig und breiten sich auf immer mehr Gebiete aus. Neben recht 
allgemein bekannten Modellen aus den Natur- und Ingenieurwissen- 
schaften haben sich ja inzwischen z.B. auch Modelle als außerordent- 
lich wertvoll erwiesen, mit denen man Vorhaben plant, Ketten logi- 
scher Schlußfolgerungen nachvollzieht und ihre logische Stichhal- 
tigkeit prüft, komplexe Folgerungen aus einfachen Grundsätzen ab- 
leitet, Sprachen analysiert, gegen die Luftverschmutzung und gegen 
Krebs kämpft, gegen oder für Zigarettenkonsum wirbt, Weltraumflüge 
simuliert usw. 

10. DER PROZESS DER MODELLBILDUNG 

Der Prozeß der Model 1 bi 1 dung geht nun, etwas ausführlicher beschrie- 
ben, folgendermaßen vor sich: Ein Ausschnitt der realen Welt werde 
zunächst verbal beschrieben. Bereits dieser erste Prozefr 1st eigent- 
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lieh ein Abbildungsvorgang der oben beschriebenen Art. Ausgangs- 
punkt ist i.a. ein außerordentlich verwickelter Strukturkomplex, 
den wir empirisch erkennen. Der Versuch der Beschreibung ergibt zu- 
nächst eine verbal formulierte Struktur, die auf unseren Abstrak- 
tionen basiert, wobei wir erstreben, daß man durch diese Abbildung 
einen Strukturmorphi smus gewinnt. Daß jedoch bereits dieser Über- 
gang i.a. keinen Strukturisomorphi smus ergibt, ist in Kreisen der 
Semantiker wohlbekannt. Diese sprechen von dem Grundsatz 
der Nicht-Allgemeinheit, welcher den wesent- 
lichen Unterschied zwischen Wirklichkeit und Beschreibung der Wirk- 
lichkeit betont. (Gleichgültig wie viel man über eine "Sache", ei- 
nen "Vorgang", eine "Eigenschaft oder irgendetwas anderes" aussagt, 
man kann i.a. nicht alles darüber aussagen. Selbst ganz einfache 
Aussagen sind nur innerhalb gewisser Grenzen wahr. Siehe [3], S.16.) 

Ausgehend von dieser ersten, semantisch bestimmten Struktur suchen 
wir dann durch Abstraktion oder durch eine Folge von Abstraktionen, 
die wir gedanklich zu einer Gesamt-Abstraktion zusammenschalten kön- 
nen, zu einer weiteren Struktur derart überzugehen, daß möglichst 
ein Strukturmorphi smus bei dem Übergang gegeben ist. Dabei erhofft 
man immer noch, daß über die Kompatibilität hinaus ein Strukturiso- 
morphismus besteht, d.h. daß man von den Bildern zu den Originalen 
und von den Relationsbehauptungen in den Bildstrukturen zu den Re- 
lationsbehauptungen in den Ori gi nal Strukturen auf eindeutige Weise 
zurückkehren kann. 

Aber leider ist eine solche Hoffnung auch in ganz einfachen Fällen 
recht fragwürdig. Im Grunde genommen bewegen sich nämlich unsere 
Aussagen nur in einem abgeschlossenen System der Abstraktionen, die 
wir geprägt haben, und nicht in dem offenen System der Wirklichkeit. 
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Dieser Erkenntnis sind sich die Menschen oft nicht bewußt, so daß 
die Abstraktionen mit der Wirklichkeit verwechselt werden. Wenn 
plötzlich die Erleuchtung auftritt, daß Aussagen sich nicht auf 
die Wirklichkeit, sondern nur auf Abstraktionen beziehen, so kann 
leicht als eine Folgerung aus dieser plötzlichen Erkenntnis der Vor 
wurf der "Manipulation" entspringen. Ein typisches Beispiel hierfür 
ist die oft gehörte Klage: "Statistiken lügen!" (obwohl natürlich 
auch hier, ebenso wie bei allen anderen Typen von Aussagen, die 
Möglichkeit einer Manipulation vorliegen kann, nämlich durch be- 
wußt einseitige, tendenziöse Auswahl der Informationen). 

11. DIE UNVOLLSTÄNDIGKEIT VON MATHEMATISCHEN MODELLEN DER WIRK- 
LICHKEIT 

Zum Schluß seien noch einmal die prinzipiellen Schwierigkeiten ei- 
ner Idealisierung eines Ausschnittes der Wirklichkeit, wie sie bei 
einer Model 1 bi 1 dung vorgenommen wird, kurz zusammengestellt. 

1. Die empirische Ori gi nal Struktur ist möglicherweise 
nicht klar abgrenzbar oder nicht intersubjektiv 
verständlich beschreibbar. 



2. Die einzelnen Abbildungen sind möglicherweise nur 
scheinbar, nur bei oberflächlicher Betrachtung oder 
nur beim heutigen Stand unserer Erkenntnis Morphismen. 



3. Ob bei einer solchen Abbildung ein Isomorphismus 
vorliegt, ist oft sehr fraglich. 

Bei einer Rückübersetzung der Erkenntnisse, welche über die mathe- 
matische Bildstruktur gewonnen wurden, in Informationen über die 
Wirklichkeit können daher drei wesentliche Fehlerquellen entstehen: 

1. Die Wirklichkeit umfaßt mehr Ori gi nal el emente oder 
andere Ori gi nal rel ati onen , als aus den mathematischen 
Bildern ersichtlich ist. 
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2. Oie Ähnlichkeit zwischen Wirklichkeit und mathe- 
matischem Bild ist unzureichend. 



3. Die abstrakte Gleichheit zwischen Wirklichkeit 
und mathematischem Bild 1st eine Fiktion. 

Daß man trotz dieser gravierenden Einschränkungen vom Modell recht 
fruchtbare Rückschlüsse auf die Wirklichkeit ziehen kann, zeigen 
die Erfolge der Praxis. Der erfolgreiche Praktiker kennt nämlich 
die Grenzen seiner Modelle und hütet sich, diese zu überschreiten. 



LITERATUR 



[1] H. BEHNKE, R. REMMERT, H.-G. STEINER, H. TIETZ (Herausgeber) 

Das Fischer Lexikon der Mathematik 1. 

Fischer-Bücherei, Frankfurt 1966 



[2] G.B. DANTZIG: 

Lineare Programmierung und Erweiterungen (Deutsche Be- 
arbeitung von A. Jaeger). 

Springer-Verlag, Berl in-Heidelberg-New York 1966 

[3] S.I. HAYAKAWA (Herausgeber): 

Wort und Wirklichkeit. 

Verlag Darmstädter Blätter, Darmstadt 1970 

[4] A. JAEGER, K. WENKE: 

Lineare Wirtschaftsalgebra. 

Teubner, Stuttgart 1969 

[5] W. LEINFELLNER: 

Einführung in die Erkenntnis- und Wissenschaftstheorie. 
Bibliographisches Institut, Mannheim 1965 




Numerical Aspects of Pollution Abatement Problems: 
Optimal Control Strategies for Air Quality Standards 
by K. O. Kortanek and W. L. Gorr, Pittsburgh 



ABSTRACT 

The Environmental Protection Agency, U.S.A., has proposed National 
Air Quality Standards for classes of air pollution, designed to protect 
public health and welfare by setting limits on levels of pollution in the 
air. They apply to all areas of the United States. Air quality standards 
give rise to an infinite number of restrictions, called the "air quality 
constraint set," on the ambient air quality of a region; that is to say, 
the concentration at every ground level point must be less than or equal 
to the standards on the average. 

Diffusion modeling bridges the gap between emissions and ambient 
concentrations and permits specification of the air quality constraint 
set. We use one particular form of diffusion model, by drawing on its 
algebraic properties, which it also shares with many diffusion models. 

We present an approach to determine regional emission regulations 
that are (1) feasible in the sense that they lead to compliance with air 
quality standards and (2) optimal in the sense that the total economic 
impact due to implementation is minimized. Under goal (1), we set forth 
a procedure more sophisticated than heuristic procedures such as those 
guaranteeing compliance on a fixed rectangular set of grid receptor points 
in a region. 

Our procedure does not presuppose a finite set of constraints 
arbitrarily set in advance, but the correct finite set is determined by 
the entire model and the data. The model will determine the location of a 
set of maximum receptor points in the region, which then will be the optimal 
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location of sampling stations. If compliance is achieved at these strate- 
gic points in the region during the designated time interval, then 
compliance is achieved at all points in the region with respect to an 
optimal air quality standard. 
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1 . Introduction 

The approach to pollution abatement by Imposition of air quality 
standards appears to be one way in which the Environmental Protection 
Agency (EPA) is following. A standard specifies a level of concentrations 
of pollutant which is not to be exceeded on the average within a prescribed 
period of time. It appears that a prescription of standards must hold at 
all ground level (two-dimensional) coordinates. Current diffusion models 
appear to fall into two main categories- -short -term or episodal, and long- 
term. Most of the diffusion models in the literature are based on a solution 

(statistical in nature) to an initial-value, boundary-value problem which 
yields a spatial distribution of a pollutant within the class of named 
multivariate probability distributions. They are differentiated through 
assumptions as to the type of emission source, functional complexity of the 
parameters, and the nature of the iteration method employed. See Herzog [33] 
and also Sutton [51], Pasquill [43], [44], Gifford [25], Turner [54], 
Martin-Tikvart [40], Fortak [23], Roberts, et al [46], and others. 

Normally, pollution policies are based on a few dominating sampling 

points in a region, which may first be divided into a grid of areas. The 

assumption then is that if pollution is below a standard at the sampling 
stations, then it is below the standard at all. points in the region. Thus 
it appears that the region requires that emission rates of pollutant sources 
be below certain levels, demanding a fortiori that the sources take what- 
ever measures necessary that appear economically feasible. 

It is the case that the 1967 Clean Air Act and the new 1970 Clean 
Air Amendments [12] are stated in terms of air quality standards. There- 
fore, part of pollution abatement efforts must rely heavily on diffusion 
modeling. 




One of the most important uses of input data for a mathematical 
Multiple- source urban air pollution model is a real time short-term pre- 
diction of concentration for the entire urban area (see Fortak [23] and 
others) . 

However, as Professor Fortak emphasizes, the problems of city 
planning are connected with the^ problem of simulation and prediction of 
long-term climatological ground-level concentration fields. "Here, time 
series of concentrations are of minor interest; instead statistics of 
observed or calculated concentration fields for given long periods of time 
are important (Fortak [23], page 9-2). 

In this paper we develop an optimization model involving estima- 
tions for the annual ground level pollution due to area or point sources 
of pollution involving approximately unimodel pollution surfaces of the 
kind referenced above. Using as decision variables the weight fraction 
of SO 2 and particulate removed at the sources, we illustrate this class 
of models by a single objective function of minimizing total cost from 
all sources while keeping total emissions below standards at all coordinates 
of the region. Since both short and long-term have been related implicitly 
by Legislation (see [12]), we also develop a selective abatement procedure 
model which, for the first time, considers short and long-term standards 
simultaneously. 

We therefore achieve a model formulation which is consistent with 
any functional relationship whatever, and in particular, well suited when 
such relationships are based on a single parameter such as time, t. No 
longer is the model dependent on the grid size, or the space between grids. 
In fact, we prove that only a finite number of points of concentration in 
the region need enter an optimal solution, as ground level coordinates for 
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the basic pollution surfaces. Thus, by appealing to recent results on 
generalized moment problems [29] , [30] , we render a continuous parameter 
problem into an equivalent one with only a finite number of parameter values 
needed. 

As already stated by several authors, the functional complexity 
of the parameters and the iteration method employed are intrinsic parts of 
diffusion modeling. It appears, however, that the literature has not yet 
recognized that optimization problems stemming from models such as the 
Martin-Tikvart , Fortak or Roberts, et al. models are instances of semi- 
infinite programming problems [6], [7], [8]. Once this identification is 
made, it is possible to give a generalized moment type problem whose 
solution is equivalent to solving a (semi-infinite) Martin-Tikvart problem. 
In fact, following Gustafson [29] and Gustafson-Kortanek-Rom [31], it is 
seen that a generalized moment problem is "dual” to the Martin-Tikvart 
problem. We thus develop models which display the functional complexity 
of paramters, and incorporate efficient iteration methods of solution, 
namely algorithms for solving generalized moment problems as developed in 
[29] and [31] , based on Newton-Raphson subroutines and implementable via 
computer codes. 

2. Pollution Abatement Problem (PAP) Formulations 

2.1 A Two-Dimensional Minimum Cost Diffusion Model 

We cite an optimization problem frequently found in the literature, 
as presented for example by Mr, Wilpen Gorr [26] , stemming from the work 
of Martin-Tikvart [40] . 



39 



The annual average ground level pollution due to N stationary 
sources is calculated as: 



( 2 . 1 ) 



X^x.y) 



N 

: £ Q, 'Mx.y) 
j=i J J 



where 



(2.2) X,j, - annual average concentration at ground level coordinate 

3 

(x,y) in units of ^gm/m , 



(2.3) 



= present average emission rate of pollutant in mgm/sec 
at source j , and 



(2.4) 



= basic pollution surface for source j, approximately 
unimodal . 



Decision variables are defined by: 

(2.5) Eg = weight fraction of SO^ removed at source j 

J 

(2.6) Ep = weight fraction of particulate removed at source j. 

j 



The semi-infinite optimization problem (finite number of variables, 
finite number of constants) becomes in this notation: 



N 

.?‘s. ( E s.) 

J-l J ^ J' 
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(2.7) subject to 
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(x,y) < S , for all (x,y)€A 
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- Eg > - E S U 1< j<N. 

j j 

In this problem we use: 

(2.8> Q = S0 o emission rate at source j 

s j 

(2.9) Cg ^E \ = cost of attaining E 

S j V S j/ S j 

(2.10) S = annual SO^ air quality standard 

(2.11) A = air pollution control region (county, etc.), 

(2.12) E^ , E^ = lower and upper bounds on the E - variables. 

J J J 

Observe that (2.7) may be further refined by functional relation- 
ships involving a single parameter such as time, t. For example, x and y 
coordinates may be converted to functions of time x(t), y(t) where distances 
are measured along trajectories. 

Note also that a similar optimization problem may be given for 
particulates involving decision variables E r , either by itself or mixed 

j 

with problem (2.7). 

Observe that problem (2.7) has as many linear constraints as there 
are coordinate points in the control region, namely an infinite number of 
them. However, the number of variables is finite, and therefore problem 
(£.7) it a semi-infinite programming problem of the type developed and 
fttudied by Charne8 -Cooper Kortanek [7], [8], and [9]. 
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2.2 A One -Dimensional Analog of Problem (2.7) 

In order to relate recent results on generalized moment theory 
of Gustafson, Kortanek and Rom to PAP, we consider a one -dimensional analog 
of optimization problem (2.7), where for convenience we use a linear 
objective function. ^ 



N 

min E Cg E g 

j = l j j 



( 2 . 2 . 1 ) 



s . t . 



T' ( 1_E s .) Q S. 

J=1 ' V J J 



for all t € I 



E s. ^ E S. 

J J 



- E s.^ - E s. 

J J 



1< j<N. 



Here I is a closed interval, taken for example to be [0,1]. The N sources 

. ^ 2 

are at distinct points in I. The functions i|fj(t) 55 e ' where 

source j would be at the point 6^ in I, are examples of unimodal surfaces 
of normal distribution type. 



2.2.2 A Restatement of Problem (2.2.1) in Convenient Notation 



Continuing our goal of relating problem (2.2.1) to moment theory, 
we write problem (2.2.1) in the following ^ equivalent form: 



inf 



N 



E y 

r=l 



r 






r 



Formulation (2.2.1) may be equivalent to a two-dimensional one via certain 
general functional relationships. In some cases» therefore, (2.2.1) may 
not in effect be a one- dimensional problem initially* 
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for all t € [0,1] 



1< r < N 

Now the y 's are the variables, p- 's the cost, and u (t) are continuous 
functions on I. Here "inf" replaces "min" because in the general theory, 
an infinum may not be a minimum. Problem (2. 2. 2.1) is analogous to problem 
D in [29] and [31], but with bounded variables. These extra bounding 
conditions give a more general problem than D, but one which can be handled 
by extensions of their numerical methods. 

2.2.3 The Moment Mathematical Mpdel - See [29] 

For numerical solution of problem type (2. 2. 2.1), we actually 
propose to numerically solve its "dual" problem. To this end we turn now 
to [29] and [311 for development of the dual moment problem. 

Let E be the linear space of functions a which are of bounded 
variation on [a,b] , continuous to the right in (a,b) , and normalized by 
the condition Qf(a) = 0 . We introduce the norm ||or|| defined by 

II “11 = Ja 1 d “( t )l 



(2.2. 2.1) subject to Sy u (t) > ^>(t) 



r=l 



. L 
y > y 

r- 'r 



- y > - y 
7 r — r 



Thus E is a Banach space. 
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Let F be the linear space of functions, continuous on [a,b] 
denoted by C[a,b] and normed by 

\\ f\\ = max \ f (t) | 
t e [ a ,b] 



Thus F is a Banach space. 

If L is a bounded linear functional on C[a,b] then there exists 
an a in E which is independent of f and such that 

L(f) = Jg f (t) do'(t) 

We therefore can introduce the pairing < , > between E and F 

< f , c*> =f b f(t)dcy(t) f e F , aeE . 

J a 

Denote by R n the Euclidean n-dimensional vector space. Choose n functions 
u^u^,...,^ from F and define the operator U:E*R n by U(Qf) = y where 

r b u (t)dff(t) = y , r = l,2,...,n 
Ja r r 

Now let $ be a fixed function in F and \i a fixed vector in R n . Then we 
can formulate the problem 

(2. 2. 3.1) Compute sup 

subject to 
and 



<$ , ot> 

U(o0 = p. 

S' 

oi 
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Let C denote the set of functions o' which are not decreasing and 
of bounded variation over [a,b] . It is easily verified that C is a closed 
cone* W.e now use the procedure outlined in [29] and formulate the pair of 
dual problems P and D 

P: max<$,o , > D: min (p,,y) 

T 

subject to U(o0 = \x subject to U (y) - $ e C* 

a ec 

Here ( , ) denotes scalar product in R n . The adjoint operator U T 
of U, : R n -*F is defined by the requirement 

(2. 2. 3. 2) (U(or),y) = <or, U T y> 

for all o' in C and y in R n . It* is easily proved that 




and that this is the only way of satisfying (2. 2. 3. 2). Next we define 
C* by the relation 

C* = [f : f s F, «ec, < f, a>> 0} 

We Immediately find 

C* = [f : f e F, f(t)>0, V t e[a,b]} 
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Now we can formulate P and D in the following way 



(2. 2. 3. 3) P: sup J^(t)dcv(t) 



subj. to u r (t)da(t) = p* r , 

s 

a 

r= 1,2,. ..,n 



D: inf S y p. 

r=l 



subj. to Hy u (t)>$(t), t e[a»b] 
r=l r 



Problem (2.2.2 .1) is a slight extension of problem D in Gustafson 
[29], having the additional bounding conditions on the y^-variables. 

2.3.4 Some Properties Leading to an Algorithm . See [29] 

Gustafson [29] and Gustaf son-Kortanek-Rom [31] consider regularised 
versions of problem D (denoted by D p ) where additional redundant constraints 
| yj < F are adjoined. These are introduced for purposes of solution 
rather than being part of the problem structure* 

Hence the moment problem associated with (2 » 2 .2.1) becomes: 

(2.3.4. 1) sup J* ^(t)dof(t) +Z)(y^ v^ - y^ v^ 

r=l 

s.t. J** u r (t) d of(t) + v^ - v r ” r * 1>2,. •• » n 

^ L u . n 

a v r , v r > 0 
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The following Theorems are related to analogous ones in [29] and 
[31] . Their proofs follow analogously also, providing one replaces 
| yj < F by the pair u^ > y^ , y f < y^ for each variable r. 



Theorem 1^ (See [29]) Assume some u.(t) > 0 ^ for all t e [0,1] , and 



that max |^(t)l/u.(t) < yV . 
t e [0,1] J J 



Then only two cases for problems 
(2.2.2. 1) and (2.3.4. 1) are possible: 



1. Both are consistent and have bounded solutions 

2. (2.3.4. D is inconsistent and (2.2.2. 1) has unbounded solutions. 

Theorem 2 The objective value of (2.3.4. 1) can be taken on for an 
integrator with at most n points of increase, i-.e. a step function. 



Remark It is this theorem which permits replacement of continuously 

parameterized problem by a discrete one. However, there is no apriori 
reason for the mass points to be ’’evenly spaced" such as the case for 
a grid system. 

The following theorem yields a system of nonlinear equations for 
the n points of increase which leads to an algorithm (see [29], [31]) 
presented in the appendix. 



1 ) 



This will be true for u^(t) 
is a 'ceby^ev system, see [28]. 



if 



the system 
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Theorem 3 (See [29]) Assume an optimal solution of (2. 3. 4.1) is given 
by: 



(i) 


mass 


at point 


t ± c [0,1] . 


, for 


i = l,2, 


(ii) 


U 

q 


of the v U > 0 , 
r 


namely 


U 

V 


U 

•*-’ V r U 

q 


(iii) 


L 

q 


of the v L >0, 
r 


namely 


L 

V s, ’■ 


L 

" ,V s L 



L ö - ö u 

i q 

Let be an optimal solution of (2. 2. 2.1). Then the 

following equations are satisfied: 

q t u 

(T3.1) ?^P i u r ( t i ) + v r " v r = ; r==1 » 2 »-**» n 

i-1 

n 

(T3.2) T< y^r-C^) = ^ (t ± ) ; i = l,2,...,q 

r=l 

(T3.3) y r = y V j = l,2,...,q U 

j J 

(T3.4) y = y^ j = 1 ,2 , . . . ,q L 

S j J 

(T3.5) Tjy u'(t ) =^> (t.) i = l,2,...,q . 

r.l r 1 1 

These necessary conditions lead to an algorithm essentially 
identical to [29] and [31] . 

As is apparent from [31] , several different states can occur, 
which are not encountered in finite linear programming. Actually, a 
variety of analytical problems arise depending on assumptions imposed on 
the functions u^,U 2 >...,u . For example, if they form a complete. Ceby^ev 

system, then the results are simpler but still permit the treatment of a 
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wide range of problems. An algorithm for numerical solution of the general- 
ized moment problem for one dimension has been developed in Gustafson- 
Kortanek-Rom [31], Gustafson- Rom [32] and Gustafson [29]. 

3 Conditional-Adaptive Planning 

In industrial engineering contexts such as multi-product produc- 
tion scheduling [36], [49], the sequential nature of production planning 
has been explicitly taken into account. Long- and short-term horizons have 
been designated, with the short-term decisions being implemented and the 
long-term providing overall guidance. By implementing the short-term plan, 
we increase the chance of a quicker response to change while implicitly 
incorporating the intertemporal character of the entire horizon (both short 
and long) to shape the current period decision. In fact, by using con- 
straint oriented devices such as ’’horizon-postures" recursively, one may 
indirectly minimize changeover and start-up effects across all horizons. 

In pollution abatement, it appears that both short and long-term 
horizons have been related via a standards approach. Short-term standards 
impose higher concentrations vis-a-vis long term ones. On the other hand, 
however, very high short-term stagnation pollution potentials have resulted 
in temporary and economically expensive control measures such as fuel 
switching, cutting back on operations, etc. 
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The problem of designing emission limitations for both short and 
long-term air quality standards can be a difficult one. Short-term 
diffusion modeling appears to be a much less precise tool than long- 
term diffusion modeling, the reason being statistical (law of averages) 
in nature. In some cases, such as with SC^» it appears that the long 
term air quality standards are dominating in the sense that PAP- 
designed regulations would also guarantee compliance with 24-hour 
standards. Here one can demonstrate with simple calculations that the 
effect of "episode’ 1 curtailments of emissions has negligible effect on 
annual average concentrations. Thus the PAP formulation is applicable. 

On the other hand, fine particulate standards appear to have the 
short term requiring the more restrictive regulations. In this case 
one might redefine the annual standard using one of the many empirical 
short term/ long term concentration relationships so that the redefined 
standard dominates the short term standard. Again the PAP formulation 
would be relevant in a practical way since we can more readily model 
annual concentrations. 

It appears at this point that the issues of changeover, horizons 
and conditional-adaptive planning are ones which cannot be avoided if there 
is interest in numerical results including implementation. We are interested 
in developing "laws" or "functional relationships" wherever possible in 
order to simplify complex relationships. We cannot \/ait^ hoover. 
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for the complete master plan which delineates the function task of each 
entity or "grid cell" of the city. We intend to develop analytical tools 
which assist in prediction and control of a system without complete dependency 
on perfect understanding, as detailed by Professors Cooper, Eastman, 

Johnson and Kortanek, e.g. see the section "A. The Modeling Effort," 
pp. B742-B745 of [48]. 

4. Optimum Emission Control Strategies Achieving Air Quality Standards 

The current best known methods for the formulation and evaluation 
of emission control strategies (see [22], [52] and [53] below) involve 
searching through a prespecified list of strategies to see which of these 
meet air quality standards at all grid points in a region, and which minimize 
total regional cost. These methods are trial and error procedures and 
involve the computer through combinatorial simulations. Such heuristic proqedur 
(see [53]) appear to be of real value in getting a start on the problem of 

establishing implementation plans. The essential ingredients of the 
method are: (1) have the user specify various emission standards for 

particular pollutants in a region, (2) for each of these compute the 
changes in pollution surfaces and test for compliance with air quality 
standards at a finite set of receptor points in the region, and (30 of 
those strategies which do comply at these points, select one which minimizes 
total regional cost and is technically feasible with respect to available 



devices . 
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We view our development on optimal control strategies to be 
undertaken along with these. We have a more sophisticated procedure to 
be sure, but one which could be used to measure the effectiveness of the 
heuristic procedure, including justifying those instances where the 
heuristics alone are enough. The optimal procedure could single out 
instances, however, where additional refinements are necessary, thereby 
reducing the chances of embarrassing consequences from a straightforward 
application of the heuristic approach to all cases. 

5. Some Methods for Further Restrictions on the E^ - Reductions : Towards 

Regulatory Policies [39] ^ 

5 . 1 Min-Max Percent Reduction, per Source 

Since we have already indicated elsewhere preliminary responses 

to potential problems that least cost strategies may impose for regulatory 

policy, [39] , we indicate some ways in which our technology could embrace 

modifications in order to obtain a more accurate picture of the total impact 

of implementation of control strategies. The E^-variables introduced in 

section 2 are the weight fraction of S0 o removed at source j , E c , or the 

2 Sj 

weight fraction of particulate removed at source j , E p . These variables 

j 

yield selective abatement procedures, and lead to regional emission regula- 
tions incorporating possibly different reductions for different sources. 

^We would like to acknowledge helpful discussions with Professors Blumstein, 
Cassidy and Walters of the School of Urban and Public Affairs. See also 
reference [1]. 
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The numerical values of these emission reductions are determined that are 
(1) feasible in the sense that they lead to compliance with air quality 
standards throughout the region and (2) optimal in the sense that the 
"total" economic impact due to implementation is minimized. 

In section 2 we introduced lower limits and upper limits for each 
source reduction E ^ , which would be tuned in practice according to a host 
of attributes of particular sources depending on the region in which they 
are located, local politics, etc. Thus, we require that each be bigger 
than or equal to a certain lower limit E ^ , and smaller than or equal to 
an upper limit, E 1 ?. In terms of a scenario at this time, an agency could 
specify the min-max reduction percentages within which we would seek 
emission reductions, E ^ , which would be now the least cost, but from among 
those meeting goal (1) . 

It is probably the case that this modification is still too fine 
for practice. However, it does lead into the following modification based 
on classifications of source types. 



5.2 - Competitive n-Person Game Relations Between the E^ 's 



In order to not rely exclusively on a "price" system to distri- 
bute total regional impact among sources, we can introduce various 
competitive game theoretic concepts which involve dominance and stability 
(including political) . Included in this approach would be a weighting 
procedure among the E^ -classes which might involve a regulator's estimate 
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of a mutually stable relationship among sources in a region. Another con- 
cept could involve the notion of "undominated" emission regulations to be 
determined so that one source cannot develop a "justified objection" 
against any of the other sources. 

In any of the modifications (1) through (4) , our available 
technology will yield systems measures for the cost-benefit consequences 
of their imposition. For example, we can obtain the change in "total 
regional cost" as measured one particular way before the introduction of 
a regulatory modification and afterwards. In this way, we can get a 

measure on possible inequities and increased control costs above and beyond 
the initial way chosen to measure "total regional cost," as for example, 
set forth in our models in sections 2 and 3. 

In summary , any of these modifications suggested by legal 
regulatory considerations would be goals required to be met in addition 
to (1) , while attempting to seek (2) , an optimal least cost solution 
under the altered circumstances. It is because of additional constraints 
and restrictions that we have dubbed our numerical solution methods. 



constrained generalized moment techniques 
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Entscheidungskriterien — Hilfe oder Hindernis beim Handeln? 
von M. Rutsch, Karlsruhe 



Man spricht heute viel von den sogenannten Entschei dungs- 
wi ssenschaf ten (Decision Sciences). Zu diesem Komplex gehören 

- Theorien und Techniken des Optimierens, 

seine verschiedenen Varianten (z.B. dynamische 
Optimierung) und die einschlägigen kombina- 
torischen Hilfsmittel ; 

- Kon trol 1 theori e und Lerntheorie, sowie die 

- eigentliche (statistische) Entscheidungstheorie. 

Wenn es um eine Sache schlecht bestellt ist, entsteht 
gewöhnlich rasch eine Wissenschaft, man denke etwa 
an den expandierenden Komplex der Erziehungswissenschaften. 
Der wissenschaftliche und vor allem der literarische Aus- 
stoß reagiert sehr ra'Sch und häufig überproportional auf 
die neue Nachfrage. Ist daraus zu schließen, daß es wie 
um das Erziehen, so auch um das Entscheiden schlecht be- 
stellt ist? Einige Tatsachen sprechen dafür: 

Wie beim Eftziehen, so reichen auch beim Entscheiden die her- 
kömmlichen Methoden, die sich unter früheren Verhältnissen 
durchaus bewährt hatten, nicht mehr aus. Hinzu kommt, daß 
man sich nicht mehr damit begnügen will, noch funktionie- 
rende Verfahren weiter zu verwenden, ohne nachgeprüft 
zu haben, warum sie funktionieren, wieweit man sich da- 
rauf verlassen kann, und ob es nicht doch etwas Besseres 
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gibt: man ist progressiv, man experimentiert ... das ge- 
hört heute zum guten Ton. Aber auch von dieser Mode- 
strömung abgesehen, macht sich ein tiefes und weitgehend 
berechtigtes Unbehagen über Entscheidungen bemerkbar, 
von denen wir hören. Jede Maßnahme, die sehr viele Men- 
schen betrifft und für sie meistens eine Belastung oder 
Einschränkung darstellt (z.B. Steuer- oder Fahrpreiser- 
höhung), wird kritisiert. Die Journalisten der verschie- 
denen Medien fragen: War das nötig? War das richtig? Wir 
hören von eklatanten Fehlentscheidungen, von sinnlosen 
oder schädlichen Entscheidungen; von Entscheidungen, die 
unausgereift und unüberlegt realisiert wurden, d.h. ohne 
reifliche Überlegung aller Konsequenzen. Man berichtet uns 
von Entscheidungen, die nicht auf einwandfreie Weise zu- 
standegekommen sind, sondern z.B. unter dem Druck von 
Interessengruppen. Die widersprüchlichsten Entscheidungen 
lösen einander in planloser Folge ab. Wir hören von Miß- 
ständen, die doch auch nichts anderes sind als das 
Resultat unterlassener oder hinausgeschobener Aktionen, 
unkoordi ni erter und verzögerter Entscheidungen. 

Bei diesem trostlosen Stand der Dinge wäre es einfach 
wunderbar, wenn die Entscheidungswissenschaften den ver- 
antwortlichen Individuen oder Gremien die Arbeit des Ent- 
scheidens abnehmen und im richtigen Augenblick die 
richtige Entscheidung aufspüren und anzeigen würden. 

Wenn das Entscheidungsproblem einen gewissen Reifegrad 
erreicht hat (d.h. wenn genügend viele Vorentscheidungen 
schon gefallen sind), besteht die Arbeit der Entschei dungs 
findung darin, ein optimales Aktionsprogramm (eine optima- 
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le "Strategie") zu berechnen. Wir werden auf dieser Tagung 
Gelegenheit haben, von neuen Algorithmen und Modellen zu 
hören, die mit großem Scharfsinn zur Analyse, Prognose 
und optimalen Planung demographischer, kommunaler, 
ökologischer, industrieller und gesamtwi rtschaftlicher 
Prozesse ausgedacht wurden und sicher nicht verfehlen 
werden, uns zu beeindrucken. Es mag sogar den Anschein 
haben, als gäbe es Lösungen, wo garkeine Probleme sind, 
oder wie Jacques Prevert so schön absurd sagt, als gäbe 
es "reponses ä tout, questions ä rien Und ich brauche 

kaum darauf hinzuweisen, daß Optimierung, immer und über- 
all, ein Fetisch zu werden droht und zum Teil schon gewor- 
den ist; so sagt der Chef in einer der Zeichnungen des 
NEW YORKER zu seinen Angestellten: "Don't just sit there - 
optimize something!" Trotz dieser kritischen Vorbehalte 
wird der Bereich, in dem OR-Techniken vernünftigerweise 
und mit Aussicht auf Erfolg zur Entwirrung und Lösung von 
Entscheidungsproblemen eingesetzt werden können und hoffent- 
lich auch eingesetzt werden, in Zukunft immer größer 
werden. Darüber dürfte es nach den programmati sehen Aus- 
führungen von Herrn KÖNZI keinen Zweifel geben. 

Prob! emati sch ist dagegen die Hilfe der Entscheidungs- 
wissenschaften bei Aufgabenstellungen, die noch nicht 
soweit aufbereitet sind, daß alle wesentlichen Eingabe- 
daten vorliegen und die optimale Aktion "nur" noch ausge- 
rechnet zu werden braucht, und erst recht in Problem- 
situationen, wo eine derartige Aufbereitung einer Dena- 
turierung gleichkäme. Kann auch da der Entscheidungstheo- 
retiker den eigentlichen Entscheidungsträger, den für die Ent- 
scheidung Verantwortlichen, ablösen? Kann er in die Rolle des 
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Entscheidenden schlüpfen wie Siegfried in die Gestalt 
Gunthers, als dieser um Brunhild warb - oder ist auch das 
nur ein Mythos? 

Wir wollen gleich sagen, für welche Art von Entscheidungen 
es nicht geht: Für die “inspirierten" Entscheidungen, die 
eine glückliche Eingebung des Augenblicks darstellen, für 
künstlerisches und allgemeiner: kreatives Schaffen. "Je ne 
pouvais rien faire d'apres une decision prise d'avance", 
sagte Franjoi s ‘ Mauri ac in einem seiner letzten Gespräche, 

Ein Kunstwerk frappiert durch das Unerwartete und Einmalige; 
seine Entstehung kann nicht als "Auswahl unter mehreren 
Formen" gedacht werden - obwohl wir zuweilen seine Genese 
(in Vorstudien) und das tastende Suchen nach ihm verfolgen 
können , 

Das instinktive Entsche i dungsverhal ten der Tiere ist stammes- 
geschichtlich erworben . Instinktives Entsche i dungs verhal- 
ten gibt es auch beim Menschen; in manchen Situationen 
ist es unerläßlich, fn anderen schädlich - immer aber ist 
es von Reflexion und Analyse bedroht, die solche in- 
stinktiven oder reaktiven Entscheidungsschemata in uns 
"löscht". Der Mensch muß mit seinen Entscheidungen immer 
wieder ganz vorn, am Nullpunkt anfangen; das ist ein Privi- 
leg und zugleich eine große Belastung. Daher ist der Mensch 
im allgemeinen froh, wenn er Entscheidungen los wird. 
Entscheidungen machen, wie wir alle aus eigener Erfahrung 
wissen, müde und reizbar. Wir können sehr viel mehr Wahr- 
nehmungen aufnehmen als Entscheidungen treffen - ich glaube 
200 bits/min gegenüber 10 - 50 bits/min - und wir ermüden 
bei 50 bits Entscheidungen pro Minute sehr schnell. Das 
an sich schon anstrengende Geschäft des Entscheidens wird 
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häufig noch dadurch erschwert, daß der Entschei dungsträger 
sich in seiner Rolle nicht wohlfühlt. Psychologen sehen es 
geradezu als Indiz für schlechte Anpassung in einer beruf- 
lichen Stellung an, wenn jemand Angst hat, selbst kleinere 
Entscheidungen zu treffen. Häufig sind pathologische Stö- 
rungen des Entschei dungsverhal tens zu beobachten, z.B. 
die Unfähigkeit, die der Position des Betreffenden ange- 
messenen Entscheidungen zu treffen, oder die Neigung, 
solche Entscheidungen durch eine übergewissenhafte Analyse 
endlos hinauszuschieben und damit zu umgehen. In dem Büch- 
lein "The Peter Principle" von L. J. PETER und R. HULL findet 
man etliche fingierte, aber zutiefst realistische Krank- 
heitsbilder zur Illustration gehemmten oder blockierten 
Entschei dungs Vermögens . 

Aus all diesen Gründen sucht der Mensch einen Freund, 

Bei chtvate v ' , Advokaten, Astrologen, Mathematiker oder sonst 
eine "Vatergestalt", die ihm seine Entscheidungen abnimmt. 

Ein Kollektiv delegiert die Entscheidungsbefugnis an ein 
Führungsgremium oder an eine Einzelperson, z.B. an eine 
Persönlichkeit, die durch die Kühnheit ihrer politischen 
Konzeption imponiert, oder an jemanden, der "Glück" hat 
oder im Rufe steht, Glück zu haben - ein wichtiger Ge- 
sichtspunkt im Altertum bei der Wahl eines Feldherrn für 
eine kriegerische Expedition. Durch solche Übertragung 
der Entscheidungsbefugnis hofft man, sich gegen Fehlent- 
scheidungen abzusichern, und manchmal hat der Berater oder 
neue Entscheidungsträger tatsächlich bessere Kenntnisse, 
größere Erfahrung oder größeres Geschick in der Lenkung 
komplexer und langfristiger Unternehmungen. 

Der Stellvertreter, dem wir das Entscheiden auftragen, ist 
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uns natürlich "verantwortlich" und in gewissem Maße haft- 
bar für die Schäden, die durch seine Amtsführung evtl, ent- 
stehen. Dadurch versucht man, ihn psychologisch und mit sei- 
ner sozialen und materiellen Existenz an den Erfolg des Un- 
ternehmens zu binden. Das ist ein sehr komplexer, um nicht 
zu sagen zwielichtiger Vertrag, eine Art Seelenhandel oder 
Pakt mit dem Teufel (Wer ist der Teufel?). Es kommt vor, 
daß die Beteiligten ein unverschuldetes Fiasko mit ihrer 
Existenz bezahlen - beispielsweise der Verantwortliche, auch 
wenn er weder (absichtlich oder fahrlässig) den Zielvorstel- 
lungen des Kollektivs zuwider entscheiden hat, noch riskanter 
vorging als dem Kollektiv akzeptabel erschien. Andererseits 
kann ein Einzelner durch seine persönliche, selbst unbegrenzte 
Haftung den durch seine Entscheidung entstandenen Schaden 
eines großen Kollektivs nicht reparieren. Zu dieser Proble- 
matik der Entscheidungsdelegation kommt die ebenfalls grund- 
sätzliche Schwierigkeit hinzu, dem neuen Entscheidungsträger 
die richtigen Maßstäbe für seine Entschei dungs|ät i gkei t mitzu- 
geben: Welche Wert- und Zle 1 vorstel 1 ungen sollen zugrundege- 
legt werden? Welche Haupt- und Nebenwirkungen sind (getrennt 
oder vereint) zu veranschlagen? Welches Risiko darf man äußer- 
stenfalls eingehen? 

Das sind, bei Entscheidungen, die ein Kollektiv betreffen, 
außerordentlich brisante wirtschaftliche oder gesellschaftli- 
che Vor-Entschei düngen . Sie müssen von den Betroffenen im 
wesentlichen selbst entschieden werden, und kein Entschei dungs 
Wissenschaftler kann ihnen diese Arbeit abnehmen. Die 
Aktionäre einer Versicherungsgesellschaft vertrauen die 
Führung der Geschäfte dem Vorstand und dem Aufsichtsrat 
an, deren Mitglieder auf Grund ihrer Fähigkeiten, Erfah- 
rungen ,Zuverlässi gkei t und wegen ihres erwiesenen .oder 
erwarteten maßvollen Risikoverhaltens ausgewählt wurden. 

Dieser unvermeidbaren, durchaus üblichen und in den meisten 
Fällen gerechtfertigten Entscheidungsdelegation schließt 




65 



sich möglicherweise eine zwei te an: Die genannten Gremien 
versuchen ihre Entscheidungsbefugnis wei terzudelegieren 
an einen Versicherungsmathematiker, wenn sie z. B. vor der 
Frage stehen: In welcher Höhe soll sich die Gesellschaft 
rückversichern, wenn bei ihr ein 10 Millionen % teures 
Schiff gegen eine jährliche Prämie von 300000 % versichert 
ist? KARL BORCH, bei dem man dieses Beispiel mit weiteren, 
amüsanten Einzelheiten und sarkastischen Bemerkungen nach- 
lesen kann (The Economics of Uncertainty, 2.13 ff.), be- 
zweifelt, daß ein noch so gelehrtes statistisches Gutachten 
mehr als numerische Details zur Beantwortung der Frage bei- 
steuern kann; das Entscheidungsproblem ist schon auf seine 
einfachste, nicht weiter reduzierbare Form gebracht und muß 
von den Leitungsgremien der Gesellschaft in dieser Form 
bewältigt werden. 

Wann immer wir eine Entscheidung delegieren, vertrauen wir 
einem z.T. unbekannten Mechanismus die Entscheidungsfindung 
und damit unser Schicksal an. und treffen damit die fol- 
genschwerste Entscheidung. Das gilt auch für die Hilfe, die 
wir von Entschei dungswi ssenschaftfta , speziell Entscheidungs- 
theorie und Entscheidungskriterien erwarten: ohne Kenntnis 
der dort verwendeten Prinzipien, ohne Prüfung, ob wir die vom 
Entscheidungsmodell benötigten Daten mit der notwendigen 
Präzision bereitstellen können, wird die erhoffte Hilfe 
u.U. eher eine Gefahr - und jedenfalls ein unkal kul iertes 
Risiko - darstellen. 

Nehmen wir an, es sei uns gelungen, alle diese vielfältigen 
und vielschichtigen Bedenken auszuräumen! Dem Entschei- 
dungstheoretiker sei es also gestattet, sich an die Stelle 
des ursprünglichen Entscheidungsträgers zu versetzen. Dann 
bleibt noch eine grundsätzliche und beunruhigende Schwierig- 
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keit: die Ungewi ßhei t . Gemeint ist die Ungewißheit darüber, 
was bei jeder der ins Auge gefaßten Handlungsalternativen 
herauskommt. Welches sind die erkenntni s 1 ogi sehen Wurzeln 
dieser Ungewißheit? Ihr Entstehungsgrund ist in erster Linie 
darin zu sehen, daß die Handlungsalternativen durch Aussagen 

a l* a 2 a n 

beschrieben werden und die Umweltzustände ebenfalls durch 
Aussagen 



b l’ b 2 * • • ‘ > b m 

beschrieben werden. Im Rahmen einer "Theorie", die unser 
einschlägiges Wissen verkörpert, können wir Implikationen, 
etwa 



ableiten. Die Aussage e.. beschreibt die gemeinsame Konse- 

' J 

quenz der Handlung a^ und des Umwel tzustandes bj nach un- 
serem besten Wissen. Das heißt: Wenn ich das freigewählte a^ 
realisiere (= die Aussage a.. "wahr mache") und weiß, daß 
b. zutrifft, dann wird die Konsequenz e.. Wirklichkeit. 

w ' J 

Wenn ich aber z.B. nur weiß, daß bj oder b k zutrifft, dann 
ist auch die Handlungskonsequenz nur ungenau bekannt: 

a i a (b. v b k ) =► e.j . v e. k . 

Natürlich sind auch die ursprünglichen Konsequenzaussagen 
e.. nicht absolut präzis und erschöpfend; aber selbst wenn wii 

' J 
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annehmen, daß si e alle Angaben in hinreichender Präzision ent- 
halten, die den von der Entscheidung Betroffenen interessieren 
dann kann doch bereits die ungenauere Aussage 



eine beängstigende Ungewißheit über den Ausgang der Aktion 
a i darstellen. Daran knüpfen sich eine Reihe von Aufgaben 
und Problemstellungen. Wir »ollen einige von ihnen kurz 
skizzieren. 



Das Ordnen der Konsequenzen gemäß unseren Präferenzen berei- 
tet Schwierigkeiten, da es sich bei ihnen prinzipiell um 
ungefähre Beschreibungen 

e^j oder e^j v e^ usw. 

von Sachverhalten handelt, die unser Wohlergehen berühren. 

Man kann sich diese ungefähren Ergebnisse etwa als Teil- 
mengen (Purtkthaufen oder Intervalle) auf einer Skala vor- 
stellen, nicht als präzise Punkte darauf. Die Schwierig- 
keiten werden verschärft, wenn wir die Aktion a^ nicht genau 
festlegen können oder bj nicht genau kennen (s.o.). Die Schwie- 
rigkeiten werden gemildert, wenn wir die Umweltaussage b^ 
präzisieren können, oder wenn uns eine bessere Theorie 
präzisere Impli kate e.. zu den Prämissen a. und b. lie'- 

' J I J 

fert. 



Die - nach Möglichkeit totale - Präferenzordnung der e^j 
dient dem Zwecke, die Aktionen zu ordnen, die eine oder 
mehrere dieser Konsequenzen hervorbringen können, und damit 
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letzten Endes zur Auswahl einer optimalen Aktion. Selbst 

wenn die e. . total geordnet sind (d.h. je zwei von ihnen 

■ u 

sind präferenzmäßig vergleichbar), wir aber nicht wissen, 
welches b. zutrifft, können wir deswegen trotzdem die 

sJ 

Aktionen nicht ohne weiteres auch total ordnen, weil ja 
jeder Aktion ein ganzes Bündel von Konsequenzen entspricht. 
Die Entschei düngskri terien haben gerade den Zweck, auch 
unter solchen Umständen aus der(Präferenz-)Ordnung der Kon- 
sequenzen eine ( Präferenz- )0rdnung zwischen den Aktionen her 
zuleiten. Wir müssen also Entscheidungskriterien ent- 
wickeln für den Fall, daß die Konsequenzen der einzelnen 
Aktionen nicht eindeutig bekannt sind. Beispielsweise 
sei der wahre Umweltzustand nur in Form der ungenauen 
Aussage F bekannt und von den zwei Handlungen a^, a 2 wisse 
man , 



bei F kann a^ zu, e^ oder auch zu ej führen, 
zu e 2 oder auch zu e£. 

Man kann F formal in vier feinere Zustände aufspalten und 
erhäl t 




Dieser schöne Trick nützt ga-rnichts für die Entscheidung 
zwischen a^ und a 2 , außer in speziellen Situationen wie 
den folgenden: 
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1. Jede- bei F mögliche Konsequenz von a^ ist besser als je- 
de Konsequenz von a 2 : dann ziehen wir a^ der Aktion a 2 vor. 
Dieses selbstverständliche Vorgehen läuft darauf hinaus, 
daß wir zwischen den Teilmengen der Menge E aller möglichen 
Konsequenzen eine Ordnung Xj einführen durch: 

{e j , 6 ^ , ... } ^ ® 2 ** • • • ^ 

falls jedes e^, e^, ... besser ist als jedes e 2 , e 2> ... 

Diese partiel le Ordnung auf £(E) übefcträgt sich in offen- 
sichtlicher Weise auf die Menge A der Aktionen und kann 
eventuell sogar, wie im obigen speziellen Fall, auf der Teil- 
menge der zur Entscheidung stehenden Aktionen (oben: a^ und 
a 2 ) eine totale Ordnung induzieren. 

2. Wir sind bereit, uns nach dem Maximin-Prinzip zu entschei- 
den. Das bedeutet: eine Aktion ist uns um so lieber, je weniger 
schlechte Konsequenzen sie hervorbri ngtn kann. Dieser Ent- 
stellung entspricht folgende Ordnung >2 au ^ $(£■)'• 

min {ej, ej, ... > > min {e 2 , e 2 , ...} , 

wo > die Präferenzordnung auf E bezeichnet. Ist X total , so 
gilt dasselbe von X 2 . Die Ordnungsrelation X 2 enthält Xj. 

3. Wir sind in der Lage, mehr über die Möglichkeit bzw. 
Unmöglichkeit der feineren Zustände bj, b 2> ... heraus- 
zufinden. Wir mögen beispielsweise feststellen, daß 

e^ und e 2 unvereinbar 



sind und. daß 

ej'sfc. a^ und e^ a^ . 
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Dann kann es keinen Feinzustand geben (man beachte, 

daß die Unvereinbarkeit von e^ und dazu nicht ausreicht). 

Nach Streichung der zu b^ gehörigen Spalte in der obigen 

Ergebni smatri x mag sich vielleicht Dominanz (im Sinne 

von , siehe Nr. 1 ) zwischen a^ und herausstel 1 en , etwa 

wenn 



> e 0 >* e n e 



2 ’ 



Die Streichung von Feinzuständen nach Nr. 3 zerstört nichts 
an der Relation ^ der Nr. 1 , sondern vergrößert sie 
höchstens. Dagegen kann die Relation >2» die schon maximal 
ist, u.U. in ihr Gegenteil verkehrt werden. Beides ist aber 
nur möglich, wenn eine Konsequenz, z.B.e2> als Ergebnis 
von a2 völlig ausfällt, wenn also im obigen Beispiel b^ und 
b2 gestrichen werden - das entspricht einer Einengung des 
Ergebnisses von a2 und F von {e2> e£} auf { e 



Die letzte Überlegung ist wichtig, wenn wir es nicht mit 
ei nem Grobzustand F zu tun haben, sondern mit mehreren 
Grobzuständen F 1 , F2 , . . . . Es liegt dann eine "ungenaue" 
Ergebni smatri x 



a i 

a 2 



E 11 E 12 
E 2 1 E 22 



vor mit Ergebnissen E.-^» 

' J 



die Teilmengen von E sind. Als Prä- 
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ferenzordnung zwischen den E. . (d.h. auf £(E)) können wir 

I J 

>1 oder >*2 zugrundelegen. Es werde angenommen, daß über 
der Menge der Grobzustände völlige Ungewißheit herrsche. 
Ähnlich wie vorhin können wir dann wieder in speziellen 
Situationen zu einer Entscheidung finden: 

1. Wir stellen möglicherweise Dominanz von a^ Uber alle 
ander Aktionen a^ (i>l) fest, d.h. 

VI» V j : E u E 1r 

Dann werden wir selbstverständlich a^ zu realisieren suchen. 
Wenn wir uns jedes Ej als eine Menge von Feinzuständen vor- 
stellen, dann ist die eben angesprochene Dominanz gleichbe- 
deutend mit Dominanz bzw. der Menge (JF. von Feinzuständen. 

J 

2. Wir entscheiden uns nach dem Maximin-Prinzip. Dann kann 
die obige Ergebnismatrix ersetzt werden durch folgende ein- 
fachere Matrix: 





*1 * 


*2 


a l 


min E ^ 


min Ej 2 


a 2 


min E£i 


min E 2 ^ 






• 



und wir können auf sie das Maxi mainkri teri um zur Entschei dungs 
findung anwenden. Ein interessantes Mischkriterium ergibt 
sich übrigens, wenn man stattdessen das Dominanzprinzip 
auf die letzte Ergebnismatrix anwendet. Dieselbe Ergebnis- 
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matrix kann ferner auch im Risikofall (s.u.) zugrundegelegt 
werden . 

Die eben besprochene Entscheidungslage - völlige Ungewiß- 
heit Uber einer Menge von Grobzuständen ¥j, ¥ 2 , ... - kann 
man sich so entstanden denken, daß zunächst ein ganz grober 
Zustand ¥ vorlag. ¥ wurde nach dem früher erwähnten Trick 
in Feinzustände zerlegt, von denen einige nach Nr. 3 ge- 
strichen werden konnten. Mehrere Feinzustände, deren Ergeb- 
nisse bei a^, a 2 »... in beliebigen Kombinationen Vorkommen 
können (d.h. keine solche Ergebniskombination wurde ge- 
strichen), werden schließlich zu einem ¥. zusammengefaßt. 

V 

Bei derartigen Umformungen und überhaupt beim Vergleich 
von Entscheidungsproblemen bemerken wir die allgemeine 
Wirksamkeit folgender Komplementär! tät : Je feiner wir die 
Fallunt.erschei düngen 



und b. 



machen, um so präziser werden natürlich die e.j, aber lfm so 
schwerer ist es auch, das wahre bj herauszufinden oder 
präzisere Wahrscheinlichkeitsaussagen darüber zu beschaffen. 
Betrachten wir etwa, als Gegenstück zur völligen Ungewißheit, 
den Risikofall für Grobzustände! Gegeben seien also 
präzise numerische Wahrschei nl i chkei teh p j , Pj,... für die 
Grobzustände ¥j , ¥ 2 Spalten wir jedes 

■ b ] * b j * ••• 



in feinere Zustände auf, so tritt an die Stelle der ein- 
deutigen Wahrscheinlichkeitsverteilung über ¥j, ¥ 2> ... 
eine Menge von Wahrscheinlichkeitsverteilungen 
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über der' wieder völlige Ungewißheit herrscht. Es 
übrigens zwei Auffassungen über diese Menge von Wahr- 
scheinlichkeitsverteilungen: man kann alle Verteilungen 
mit T P 1 ^ = P- betrachten, oder alle Verteilungen mit 
k 3 3 

= p". für jeweils ein k = k(j), j = 1,2,.... 

Die völlige Ungewißheit kommt also durch die Hintertür 
wieder herein und muß durch ein Kriterium wie das Maxi- 
minkriterium bewältigt werden. Gleichbedeutend ist die An- 
wendung des Bayesschen Kriteriums auf die Entschei dungs- 
matri x 

Pi P2 

h *2 



min Ejj min 

min E g j min ^2 



a l 

a 2 



unter Verwendung einer erwartungstreuen Nutzenfunktion auf 
E. Das ist der Sinn des von G. MENGES so genannten Adap- 
tationskriteriums oder des SCHNEEWEISS-Kri teri ums . 

Wie wir jetzt gesehen haben, kann sich eine zunächst vor- 
handene präzise Wahrscheinlichkeitsverteilung zu einer 
(begrenzten) Ungewißheitssituation "verwischen" oder 
"verwackeln", wenn man nämlich die Zustände verfeinert. 
Setzt man aber in der statistischen Entscheidungstheorie 
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nicht voraus, daß auf dem Raum auch der feineren Zustände 
eine wohl besti mmte Wahrschei nl i chkei ts vertei 1 ung vorliegt, 
oder verspricht man nicht gar, eine solche zu beschaffen? 

In der Tat erwecken einige Entschei dungstheori en diesen 
Anschein, beispielsweise diejenige von LEONARD J. SAVAGE, 
die in verfeinerter Form auch auf dieser Tagung zur Sprache 
kommen wird. Der Entschei dungs theoreti ker legt dem (fik- 
tiven) Individuum, dem er bei der Entscheidung helfen 
will, Lotterien zur Auswahl vor und destilliert aus 
dem (hoffentlich) konsistenten Wahl verhal ten des 
Individuums dessen persönliche Wahrscheinlichkeitsver- 
teilung. Schauen wir uns diesen Prozeß einmal etwas näher 
an, um herauszufinden, ob er tatsächlich zu einer präzisen 
Wahrscheinlichkeitsverteilung führt! 

Wir benutzen zu diesem Zweck ein sehr einfaches Modell, 
dessen zentrale Figur ein Zufallsexperiment mit zwei 
möglichen Ausgängen ist. An jeden Ausgang können wir 
ein Ergebnis e^, e^ usw. knüpfen. Mit 

e l ° e 2 

bezeichnen wir die Lotterie, bei der das betreffende 
Experiment einmal durchgeführt wird und beim Eintreffen 
des ersten bzw. zweiten Ausgangs e^ bzw. e^ als Preis 
ausgehändigt wird. Wir verwenden nur zwei Ergebnisse, ein 
"gutes", e, und ein "schlechtes", 0: 

0 e . 



Daraus können wir aber mit Hilfe der Verknüpfung ® unendlich 
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viele zusammengesetzte Lotterien bilden, z.B. 

0 ° 0 , e « e , 0 © e , e © 0 , (0 ° e) © e usw. 

Auf diese Weise entsteht eine freie Wortalgebra mit der Ver- 
knüpfung © und dem Alphabet {0,e} , auf der eine Präferenz- 
ordnung definiert sei. Letztere entspricht dem Wahl verhal ten 
des Individuums gegenüber den Lotterien. Die geforderte 
Konsistenz dieses Wahl verhal tens drückt sich in einigen 
Relationen aus wie 

0 * 0 ~ 0 



oder 



0 • e -< e © 0 -* ( 0 © e ) © e < (e © 0) © e 

oder der PFANZAGLschen *Bi symmetri e und noch anderen Relationen. 
Wir illustrieren diese Konsistenzforderungen an den ersten 
Schri tten ties Auswahl prozesses und zeigen, wie sie eine 
sukzessive numeri sehe, Impl ementati on der Wahrscheinlichkeit 
ermögl i chen . 

Zunächst hat sich unser Individuum zu entscheiden für 

0 o e y e 0 0 oder 0 © e K e © 0. 

Das ist gleichwertig zur Frage: Dünkt ihm der zweite Aus- 
gang wahrschei nl i eher als der erste, oder ist es umgekehrt? 

Und für die zunächst noch fraglichen numerischen Wahrschein- 
lichkeiten p des ersten Ausgangs und q = 1-p des zweiten 
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Ausgangs bedeutet es die Festlegung auf 

p<q, d.h. p< ^ oder p>q, d.h. p> ^ . 

Wir nehmen an, das Individuum habe sich für die zweite Mög- 
lichkeit entschieden. Dadurch hat es für die dreigliedrigen 
Lotterien auch schon entschieden, daß 

0 ® 0 ® 0<0 ® ( 0 ® e ) ■< 0 ® (e • 0) < ( 0 ® e ) © o 

■< (e ® 0) ® 0. 

Frei entscheiden kann sich das Individuum aber noch für (oder 
gegen ) 

(e®0)® 0 )* ( 0 ® 0) ® e . 

Hat es sich dafür entschieden, dann verfeinert bzw. erweitert 
sich die obige Präferenzkette zu 

... (0®e)® 0-<0®e~(0®0)® e-<(e®0)® 0. 

Jetzt ist nur noch offen, ob man die Lotterie (0 ® e) » e 
zwischen die letzten beiden Glieder der Präferenzkette steiler 
will oder rechts davon. Wenn auch das geklärt ist, sind die 
dreigliedrigen Lotterien alle in eine totale Ordnung gebracht., 
Was bedeuten die letzten beiden Entschei düngen . für die nume- 
rischen Wahrscheinlichkeiten p und q? 

(e o 0) ® 0 >- (0 * e) * e bedeutet p >pq+q und somit 

® 0) ® 0 > (0 ® 0) « e folgt notwendig aus der vorige 
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Prä ferenzent sehe i dung. 

(0 • e) o e > (e ® 0) » 0 > (0 • 0) « e bedeutet 
pq + q > p > q , d.h. y ( /5"- 1 ) • 

Schon diese allerersten Schritte des Präferenzentscheidungs- 
prozesses lassen erkennen, daß die unbekannten Wahrschein- 
lichkeiten (wir hatten es oben nur mit einer zu tun, nämlich 
mit p) mit wachsender Genauigkeit in immer kürzere Intervalle 
eingeschlossen werden - daß sie aber nie mit endgültiger 
Präzision festliegen. Wahrschei nl i chkei tsf nterval 1 e statt 
exakter, punktförmiger Wahrscheinlichkeiten wurden auch von 
B. 0. K00PMAN und von I.J. GOOD untersucht, die mit "oberen" 
und "unteren" Wahrscheinlichkeiten arbeiten. Ebenso wie mit 
ungenauen Ergebni s angaben muß die Entscheidungstheorie auch 
mit solchen unscharfen Wahrscheinlichkeitswerten umzugehen 
lernen. Solange sie das nicht kann, besteht die Gefahr, daß 
sie die groben empirischen Informationen mit vorgetäuschte?? 
Präzision zu einer "optimalen" Entscheidung verarbeitet. 




Kommunale Anwendungen 
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Digitale Simulation der räumlichen Stadtentwicklung 
von J. Meise und M. Wegener, Frankfurt a. M. 



1 . Simulation sozioökonomi scher Systeme 

Mathematische Simulation ist eine Form des wissenschaftlichen 
Experiments , bei der 9 anstatt mit der Wirklichkeit selbst 9 mit 
einem numerisch- logischen Modell der Wirklichkeit experimen- 
tiert wird. 



Durch ihren experimentellen Ansatz unterscheidet sich mathe- 
matische Simulation von analytischen Lösungsverfahren mathe- 
matischer Systeme, mit denen eine Variable des Systems aus 
den übrigen Variablen allgemein errechnet wird. Simulation 
bedeutet dagegen die Ermittlung des Systemzustands für eine 
spezielle Konstellation von Werten der Systemvariablen. Eben- 
so unterscheidet sioh mathematische Simulation von Optimie- 
rungsverfahren, mit denen auf analytischem Wege diejenigen 
Werte der Systemvariablen festgestellt werden, die eine Ziel- 
funktion unter Berücksichtigung der einschränkenden Bedingun- 
gen des Systems am besten erfüllen. Die Ermittlung einer op- 
timalen Lösung ist in der Regel nicht das Hauptziel der Simu- 
lation. Selbstverständlich dienen auch Simulationsverfahren 
letzten Endes der "Optimierung” des untersuchten Systems. We- 
sentliches Ziel der Simulation ist jedoch die Gewinnung von 
Erkenntnissen über das Verhalten des beobachteten Systems un- 
ter unterschiedlichen Bedingungen. 

Mathematische Simulationsverfahren gewinnen wachsende Bedeu- 
tung für die Erforschung und Planung sozioökonomisoher Syste- 
me. Die Organisation der menschlichen Gesellschaft als Ganzes 
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und ihrer Teilbereiohe besitzt Eigenschaften eines komplexen 
dynamischen Systems« Viele ihrer Wirkungs Zusammenhänge sind 
zirkulär und können als positive oder negative Rückkopplungs- 
schleifen interpretiert werden« Es ist möglich, Teilbereiche 
des sozialen und ökonomischen Gesamtsystems in Form von Dif- 
ferential- oder Differenzengleichungen darzustellen« 



Allerdings entziehen sich solche Gleichungs Systeme häufig ei- 
ner analytischen Lösungt die meisten beobachteten Abhängig- 
keiten sind prinzipiell nichtlinear (Wachs turns Vorgänge, Ver- 
haltenswahrscheinlichkeiten, Verteilungs- und Kostenfunktio- 
nen) ; viele von ihnen sind sogar nicht glatt oder nloht stetig, 
haben dichotomischen oder Stufencharakter« Zudem kann bei nur 
geringer sachlicher, räumlioher und zeitlicher Differenzierung 
des Modells die Zahl der Gleichungen des Systems sehr groß 
werden« 



Simulations verfahren dagegen erlauben es, aueh umfangreiche 
Systeme mit beliebigen numerischen und loglsohen Verknüpfungen 
mathematisch einfach und anschaulich in ihrer seitlichen Ent- 
wicklung darzustellen« 



Aber nicht nur deshalb finden mathematische Simulations verfah- 
ren wachsende Anwendung auch im Bereich der sozioökonomisohen 
Planung« Der experimentelle Charakter der Simulation ent- 
spricht in besonderer Weise dem iterativ ablaufenden Ent Schei- 
dungsprozeß bei der Gestaltung der menschlichen Umwelt« Dieser 
unterscheidet sich, trotz prinzipieller Ähnlichkeit der Zlel- 
Mittel-Relatlon, nicht unwesentlich vom Ablauf überwiegend 
teohnisoh-wirtsohaftlioher Planungent Während dort in der Regel 
ein einziges, eindeutig und operational formuliertes Planungs- 
siel vorliegt, muß bei Planungen im sozioökonomisohen Bereich 
typischerweise von mehreren, nicht miteinander vereinbaren und 
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nloht einmal in sieh selbst widerspruchsfreien, komplexen Ziel- 
bündelungen ausgegangen werden, die won verschiedenen Gruppen 
der Gesellschaft mit unterschiedlichem Nachdruck und häufig 
auf unterschiedlichen Ebenen des Problemverständnisses vertre- 
ten werden. Der Weg von Erkennen eines Planungsproblens bis 
zur sohließliohen Entscheidung für eine konkrete Maßnahme zu 
seiner Lösung ist hier ein mühsamer und langwieriger Prozeß, 
in dem nicht nur die Lösungen, sondern auch die Ziele mehr- 
fach korrigiert werden müssen, bis ein befriedigender Aus- 
gleich zwischen den divergierenden Zielen und den zur Verfü- 
gung stehenden Mitteln gefunden worden ist. 

Experimente mit Simulationsverfahren können diesen Prozeß der 
schrittweisen Annäherung an eine Problemlösung in hervorragen- 
der Weise unterstützen. Als besonderer Vorteil erweist sich 
dabei, daß die Arbeit mit einem Simulationsmodell relativ vor- 
aussetzungslos, d.h. mit geringem Informat ions stand über das 
Planungsproblem selbst, über die Konstellation der Ziele und 
ihre möglichen Konflikte begonnen werden kann. Optimierungs- 
verfahren erfordern dagegen als ersten Schritt die Formulie- 
rung einer hier notwendigerweise mehrdimensionalen Zielfunk- 
tion mit all ihrer grundsätzlichen Problematik. Bei der Simu- 
lation erfolgen Bewertung und Auswahl mit gutem Grund außer- 
halb des eigentlichen Modells. Stattdessen soll die Arbeit 
mit dem Modell einen Lernprozeß über die Zusammenhänge des 
untersuchten Systems, über die Folgen und Wechselwirkungen 
von Planungseingriffen Ingangsetzen, der es ermöglicht, itera- 
tiv zu immer "besseren" Lösungen vorzustoßen« 

In jüngster Zelt mehren sich die Bemühungen, das Konzept der 
mathematischen Simulation auch auf die Probleme der langfri- 
stigen Stadtentwicklung zu übertragen. Die Stadt wird dabei 
verstanden als ein komplexes, dynamisches System sozialer, 
ökonomischer und technischer Beziehungen, das sich ln seiner 
räumlichen und zeitlichen Dimension verändert. Die Veränderun- 
gen des Systems können in begrenztem Maße durch Planungsein- 
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griff« beeinflußt oder gesteuert werden. Die langfristige Pla- 
nung dieser Eingriffe ist Stadtentwioklungsplanung* 

2. Bas POLIS-Modsll 

2.1 Zielsetzung 

Seit 19^9 wird an Battelle« Institut an einen Modell als Ent- 
sohe idungshilfe für die Stadtentwioklungsplanung gearbeitet* 

Bas Modell setzt sieh aus drei Teilnodellen zusanmen 9 die 
stufenweise aufeinander auf bauen (Abb* l)s 

- Bas erste Teilnodell (Prognosemodell) soll die Möglichkei- 
ten und Grenzen des Bevölkerung«- und Wirtschaftswachstums 
für die Stadt als Ganzes prognostizieren* 

- Bas zweite Teilmodell (Simulationsmodell) soll die Auswir- 
kungen alternativer Planungs Varianten der räumlichen Stadt- 
entwioklung erkennbar machen* 

- Bas dritte Teilmodell (Bewertungsmodell) soll durch die 
Bereitstellung eines überschaubaren Bewertungsverfahrens 
die Rationalität der Entscheidung zwischen den Planungsva- 
rianten erhöhen helfen* 

Von den drei Teilmodellen liegt jetzt das Zweite Teilmodell, 
das Simulationsmodell, in anwendbarer Form vor* Bis erste Er- 
probung des Modells mit Baten einer konkreten Stadt ist in- 
zwischen abgeschlossen* Hierzu stellte das Stadtplanungsamt 
der Stadt Köln Baten über die Bevölkerung«- 9 Wirtschaft«- 9 Be- 
bauung«- und Verkehrsstruktur Kölns zur Verfügung* Bie «eitere 
Erprobung und Anwendung des Modells für Planungsaufgaben Kölns 
erfolgt, in enger Abstimmung mit dem Stadtplanungsamt, im Rah- 
men eines Förschungsauftrags des Bundesministeriums für Städte- 
bau und Wohnungswesen* 

Bas Simulationsmodell POLIS ist eine Entscheidungshilfe bei 
der langfristigen Planung der räumlichen Stadtentwioklung* Es 
macht die vielschichtigen Auswirkungen alternativer Kombina- 




t Ionen Ton Flanungsmaßnahmen auf den Gebieten Flächennutzung 
und Verkehr in ihrer gegenseitigen Abhängigkeit unter Berück- 
sichtigung des Zeitablaufs erkennbar und liefert Bas is Infor- 
mationen für die mittel- bis langfristige Investitionsplanung 
der Stadt. 

Jede Stadt steht laufend vor grundlegenden Entscheidungen zur 
räumlichen Stadtentwicklung. Vo soll ein neues Wohngebiet an- 
gelegt verden? Soll man in der City weitere Hochhäuser zulas- 
sen? Sollte man mehr für den Straßenbau ausgeben? Oder den 
U-Bahn-Bau beschleunigen? Oder Tiefgaragen in der City bauen? 
Die Diskussion über solche Fragen wird nicht nur durch wider- 
streitende Interessen verdunkelt 9 sondern auch durch den Um- 
stand, daß niemand verläßliche Angaben über die wahrscheinli- 
chen Folgen der möglichen Handlungsalternativen machen kann. 

Der Vergleich zwischen ihnen, eine Bewertung ihrer Vor- und 
Haehteile nach rationalen Kriterien ist unmöglich. Die Folge 
ist bestenfalls die Lähmung aller Entscheidungen, im schlim- 
meren Fälle die Entscheidung für einen schlecht geprüften Plan. 

Hier setzt das POLIS-Modell ein. Es ermöglicht, schnell und 
mit geringem Aufwand eine Vielzahl von Handlungsalteraativen 
der Stadtverwaltung experimentell durchzuspielen. Es trägt 
durch eine differenzierte Darlegung der wichtigsten technischen 
ökonomischen und sozialen Konsequenzen der einzelnen Entsohei- 
dungsaltemativen dazu bei, diejerigen Kombinationen von Pla- 
nungsmaßnahmen auszuwfihlen, welche den Zielvorstellungen und 
finanziellen Möglichkeiten der Gemeinde am besten entspricht. 

2.2 Die Stadt als System 

Die Teilsysteme Flächennutzung und Verkehr und die Beziehungen 
zwischen ihnen bilden zusammen das r&umlioh- zeitliche System 
Stadt. Es erhält seine räumliche Dimension durch die Eintei- 
lung des Stadtgebiets und seines Umlands in eine Anzahl von 
Teilzonen, seine zeitliche dadurch, daß nicht nur ein einziger 
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Systemzustand untersucht wird, sondern die Veränderung des 
Systems über mehrere Zeitschritte (Perioden) hinweg. 

Die Beziehungen zwischen Flächennutzung und Verkehrsaufkommen 
sind seit langem bekannt und bilden die Grundlage jeder Ver- 
kehrsprognose. Aue prognostizierten Flächennutzungsdaten der 
Teilzonen des Untersuchungsgebiets wird auf die zukünftigen 
Verkehrsströme zwischen ihnen geschlossen. Der Ausbau des Ver- 
kehrssystems erfolgt entsprechend der prognostizierten Ver- 
kehrsnachfrage. Weit weniger erforscht sind die Einflüsse, die 
vom Verkehrssystem, also dem Verkehrsangebot her, auf die Flä- 
chennutzung einwirken. Dabei gehören diese Einflüsse zu den 
Hauptantriebskräften der räumlichen Stadtent Wicklung. Im 
POLIS-Modell werden sie mit Hilfe der Variablen "Zugänglich- 
keit” erfaßt. "Zugänglichkeit" ist ein Maß für die Lagegunst 
einer Zone in bezug zu den Aktivitäten in allen Zonen unter 
Berücksichtigung der vorhandenen Verkehrs Verbindungen. Eine 
Verbesserung der Erschließung einer Teilzone erhöht ihre Zu- 
gänglichkeit und damit ihre Attraktivität für die Errichtung 
von Wohnungen und Betrieben, und, sofern sie über Baulandre- 
serven verfügt, wird sie, ceteris paribus, mehr Neubautätig- 
keit auf sich ziehen als andere, weniger gut erschlossene Zo- 
nen. Als Folge erhöht sich das Verkehrsaufkommen von und zu 
dieser Zone, das Verkehrssystem muß weiter ausgebaut werden, 
und so weiter ... (Abb. 2). 

Dieser Kreislauf der "Stad tent Wicklung" findet sein Ende an 
einer Grenze, die durch die technischen Möglichkeiten bzw. die 
Bereitschaft der Gesellschaft, diese Möglichkeiten zu erproben, 
gesetzt wird. Die meisten großen Städte haben diese Grenze 
heute bereits erreicht! Das Verkehrsaufkommen übersteigt die 
Kapazität der Verkehrswege, die Zugänglichkeit der am stärk- 
sten belasteten Teilsone sinkt trotz verstärkter Verkehrsin- 
vestitionen. Gleichzeitig beansprucht gerade der Ausbau der 
Verkehrsanlagen immer größere Anteile des knappen Badens, der 
damit anderen Nutzungen entzogen wird. 




2.3 Logik des Modell« 



Mo genannten Teilsysteme und Zusammenhänge bilden das logi- 
sche Grundgerüst des POLIS-Modells. Der Zustand der einzelnen 
Teilzonen des Untersuchungsgebiets (interne Zonen) wird dureh 
Datengruppen wie "Einwohner nach Altersgruppen", "Beschäftig- 
te nach Wirtschaftsbereichen”, "Wohnungen nach Baualter und 
Bauzustand 11 und "Flächen nach Hut zungsarten " , der Einfluß der 
Umlandzonen (externe Zonen) durch ihre Einwohner, Arbeits- 
plätze und Erholungsgebiete im Modell repräsentiert. Das Ver- 
kehrssystem wird durch Streckendaten der wichtigsten Teil- 
strecken der beiden konkurrierenden Verkehrsarten, des öffent- 
lichen Personennahverkehrs (OPNV) und des Individualverkehrs 
auf der Straße (PKW) im Modell dargestellt. Die Streckendaten 
enthalten u.a. Angaben über Streckenart, Streckenlänge, Reise- 
zeit, Kapazität, verkehrende Linien und Zugfolge. 

Diese Inventardaten werden mit ihren Ausgangswerten in das 
Modell eingegeben. In der Folge unterliegen sie verschiedenen 
Arten von Veränderungen (Abb. 3 ): Einwohner bekommen Kinder, 
ziehen um, zu oder fort, altern und sterben, die Anzahl der 
Beschäftigten nimmt zu oder ab, ihre Verteilung auf die ein- 
zelnen Wirtsohaftsbereiche ändert sich mit dem Strukturwandel 
der Wirtschaft, Wohnungen werden neu- oder umbebaut, abgeris- 
sen oder in Büros umgewandelt, Flächen werden anderen Nutzung 
gen zugeführt, Straßen oder öffentliche Verkehrslinien werden 
neu- oder ausgebaut, auf anderen Linien wird der Verkehr ver- 
stärkt, verringert oder eingestellt. 

Abfolge und Ausmaß der Veränderungen werden durch Vorgaben ge- 
steuert, die für jeden Simulationslauf frei gewählt werden 
können. Die Vorgaben enthalten Annahmen über die Entwicklung 
technischer Planungsgrößen und Richtwerte, Kosten- und Finam- 
zierungsdaten sowie Parameter, die, Ergebnis modellertemer £ 
empirischer Untersuchungen, Lage und Verlauf einzelner Modell- 
funktionen festlegen. Zu den Vorgaben gehören auch die Ergeb- 
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niese des ersten Teilmodells, des (Makro-) Prognosemodells, 
bzw. bis zu dessen Fertigstellung Annahmen über den Gesamt- 
einwohner- und Besohäftigtenzuwachs des üntersuohungsgebiets 
sowie den Einwohner- und Besohäftigtenzuwachs der Umlandzonen. 
Die wichtigsten Vorgaben jedoch sind die Planungsmaßnahmen • 

Planungsmaßnahmen sind alle bewußten Eingriffe der Stadt oder 
eines anderen Planungs trägere in die Stadtentwioklu ng. Gegen- 
wärtig unterscheidet das Modell drei Gruppen von Planungsmaß- 
nahmen. Die erste Gruppe bilden institutionelle Maßnahmen, d.h. 
in erster Linie Maßnahmen der Bauleitplanung wie Festlegung 
von Art und Maß. der zulässigen Nutzung in jeder Zone. Diese 
Flächenausweisungen bedeuten allein noch keine eigentliche 
Investition oder Bautätigkeit, sie setzen lediglich den Rah- 
men für die im Modell simulierten Entscheidungen der privaten 
Investitoren und Bauherrn. Die zweite Gruppe von Planungsmaß- 
nahmen sind direkte Eingriffe der Stadt oder eines anderen 
Planungsträgers in die Flächennutzung der Zonen. Das Modell 
verarbeitet ca. 60 verschiedene Arten von Planungsmaßnahmen 
auf den Gebieten 

- Wohnungsbau, 

- Gewerbeansiedlung, 

- Gerne inbedarf Seinrichtungen, 

- Erholungseinriohtungen, 

- Einkaufseinriohtungen, 

- Ruhender Verkehr, 

- Naherschließung, 

- Ver- und Entsorgung, 

- Städtische Baulandreserve, 

- Stadterneuerung. 

Bei allen Maßnahmen kann für jede Zone ein eigenes, zeitlioh 
abgestuftes Bau- und Dur ohführungs Programm vorgegeben werden. 
Die dritte Gruppe von Planungsmaßnahmen schließlich sind die 
Verkehrsmaßnahmen. Das Modell erlaubt die Eingabe beliebiger, 
zeitlioh abgestufter Programme zum Ausbau der öffentlichen 
Verkehrsmittel und des Straßennetzes, gegebenenfalls auch 
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die experimentell« Einführung neuartiger, bisher nicht vor- 
handener Verkehrssysteme (Abb. 4, 5). 

Di« Simulation einer Periode beginnt beim Verkehrsangebot am 
Anfang der Periode* Aus den Ausgangsnetzen und den Verkehrs- 
bauten der zurückliegenden Perioden ergeben sioh die Verkehrs- 
netze der Periode* Aufgrund des Verkehrsangebots und der vor- 
handenen Plächennutzungsverteilung werden Zugänglichkeitswerte 
aller Teilzonen berechnet* Aus Zugänglichkeit und anderen 
Attraktivitätsmerkmalen der Zonen wird die wahrscheinliche Ver- 
teilung der für die betreffende Periode erwarteten Bautätig- 
keit auf die zur Verfügung stehende Baulandreserve ermittelt* 

Es folgt die Fbrt schrei bung des Einwohner- und Beschäftigten- 
bestandes der Zonen unter Berücksichtigung der Bautätigkeit» 
der natürlichen Bevölkerungsentwicklung, der Veränderung der 
Fläohenanf orderungen und Be legungs Ziffern sowie die Berechnung 
sämtlicher Folgefläohen* Dann kann aus den mittleren Einwohner- 
und Beschäftigtenzahlen der Periode das Verkehrsaufkommen so- 
wie die wahrscheinliche Verteilung der Verkehrsströme auf die 
Teilstrecken beider Verkehrsarten berechnet werden* Damit ist 
die Simulation einer Periode beendet* Alle Inventardaten ha- 
ben neue Zustandswerte erhalten, die nächste Periode kann be- 
ginnen» 

2*4 P*ftfp»*«»« ufbau 

Das Programm des Polis-Modells besteht aus zwei selbständigen 
Hauptprogrammen mit je mehreren Unterprogrammen (Abb« 6)« Das 
erste Programm POLIS 1 enthält die Programmtelle zur Analyse 
einer Verkehrs Variante« Es erzeugt ein Übergabeband mit allen 
notwendigen Daten sämtlicher Ausbaustufen beider Netze einer 
Variante« Mit diesem Band als Eingabe bandiwird im zweiten Pro- 
gramm POLIS 2 die eigentliche Simulation der Stadtentwieklung 
durchgeführt« 

Das Netzanalyseprogramm folgt im wesentlichen gebräuchlichen 
Verfahren der Verkehrsplanung* Zunächst wird aus den Strecken- 
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d&ten des Grundnetze« und etwaiger Hetzveränderungen das zu 
untersuchende Hetz zusammengestellt. Diese« Hetz wird in 
Form der sogenannten Hetzbeschreibung für die Weiterverarbei- 
tung aufbereitet. Die Hetzbeschreibung ist die Grundlage für 
die Ermittlung der "Bäume" der kürzesten Hege aller internen 
zu allen internen und externen Zonen sowie der für die Vege 
benötigten Reisezeiten. Das Verfahren wird für beide Hetze in 
gleicher Weise angewendet. Ein Unterschied besteht lediglloh 
darin, daß im Hetz der öffentlichen Verkehrsmittel bei der 
Wahl des kürzesten Weges die je nach Linienführung und Zugfol- 
ge unterschiedlichen Ums teige Zeiten berücksichtigt werden müs- 
sen. Bei den Fkw-Fährten dagegen werden die erforderlichen 
Zeiten für das Aufsuchen eines Parkplatzes in der Eielzone zu 
den eigentlichen Reisezeiten hinzugerechnet. 

Die Reisezeiten sind die Grundlage für die Berechnung der Zu- 
gängliohkeitswerte der Zonen. Reisezeitfaktoren als Funktion 
der Reisezeit drücken den Widerstand aus, der auf dem jewei- 
ligen Verkehrsnetz zwischen zwei Zonen zu überwinden ist. Je 
kleiner dieser Widerstand ist, desto mehr trägt eine Zone zur 
Zugänglichkeit der anderen Zone bei. Andererseits hängt der 
Beitrag einer Zone zur Zugänglichkeit einer anderen von der 
Größe der Aktivitäten in ihr ab. 

Als nächstes erfolgt die Simulation der Fläohennutzungsent- 
wicklung. Der Einwohnerzuwachs einer Zone setzt sich zusammen 
aus natürlicher Bevölkerungsentwicklung und Wanderungssaldo. 
Hit Hilfe der für die Gesamtstadt prognostizierten Wohnraum- 
belegung läßt sich der Gesamtzuwachs an Wohnraum in einer Pe- 
riode berechnen, der erforderlieh ist, wenn der vorgegebene 
Gesamteinwohner Zuwachs realisiert werden soll. Die Verteilung 
des Gesamtwohnraumzuwachses auf die Zonen erfolgt, wie alle 
Verteilungsoperationen des Modells, proportional zu dem je- 
weiligen "Entwicklungspotential" der Zone. Das Entwioklungs- 
potential ist allgemein die mit der Attraktivität der Zone ge- 
wichtete Kapazität, in Falle des Wohnraumzuwaohses die Flä- 
ohenreserve der Zone, gewichtet mit einer Funktion der Zu- 
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gänglichkeit und anderer Attraktivitätsmerknale» Ist der Zu- 
wachs an Wohnräunen bekannt 9 kann die Einwohnerkapazität der 
Zone am Ende der Periode berechnet werden» Dabei ist nicht nur 
das Absinken der Wohnraumbelegung im Stadtdurchsohnitt zu be- 
rüoksichtigen 9 sondern auch die für ein Wohngebiet mit einer 
bestimmten Altersstruktur und einem bestimmten Baualter und 
Zustand der Wohnungen charakteristischen Abweichung vom Stadt - 
durchsohnitt. 

Auf ähnliche Weise wird der Zuwachs an Arbeitsplätzen im pro- 
duzierenden Gewerbe und im Dienstleistungsgewerbe auf die Zo- 
nen verteilt. Anschließend werden mit Hilfe spezifischer Flä- 
ohenbeiwerte die Bruttogeschoßflächen neu zu errichtender 
Wohn- und Nichtwohngebäude sowie das erforderliche Nettobau- 
land, die Flächen für Gerne inbedarf seinrichtungen 9 die Grün- 
flächen und Erschließungsflächen berechnet» Für den Fäll f daß 
die Nachfrage nach Bauland das Angebot in einer Zone über- 
steigt 9 können besondere Maßnahmen zum Ausgleich der Öbemach- 
frage vorgesehen werden» Es folgen die Fortschreibung der Ein- 
wohneraltersgruppen, der Beschäftigten in den einzelnen Wirt- 
schaftsbereichen 9 der Wohnungen und Flächen sowie die Aufstok- 
kung der Flächenreserve durch Abriß vorhandener Bausubstanz» 

Dann kann das Verkehrsaufkommen der Periode berechnet werden» 
Die Fahrten ln beiden Richtungen zwischen zwei Zonen werden 
dabei direkt aus den mittleren Einwohner- und Beschäftigten- 
zahlen der Periode unter Benutzung der Reisezeitfaktoren be- 
rechnet» Der Anteil der Fährten mit öffentlichen Verkehrsmit- 
teln zwischen zwei Zonen wird nach dem Reisezeltverhältnis ge- 
schätzt, wobei die Modal-Split-Funktion von Periode zu Periode 
entsprechend dem Ansteigen des Motorisierungsgrads verschoben 
wird» 

Die Fahrten werden abschließend auf die Teilstrecken der er- 
mittelten kürzesten Wege umgelegt» Durch Zusammenführen von 
Strepkendaten und Teilstreckenbelastungen können Leistungsdaten 
der beiden Verkehrsnetze wie Ausnutzung des Platzangebots bzw» 
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Kapazitätsauslastung der Teilstreoken, Personenkilometer und 
Reisezeiten nach Strecken- bzw« Straßenarten ermittelt werden« 
Kapazitätsübersehreitungen im Straßenverkehr haben Stauungen 
zur Folge, denen duroh Zuschlag zu den Reisezeiten Rechnung 
getragen wird« 

Die Ausgabe der Simulationsergebnisse jeder Periode erfolgt 
auf dem Zeilendrucker in Form von Tabellen Diagrammen« Au- 
ßerdem kann in jedem Stadium der Simulation die räumliche Ver- 
teilung beliebiger Modellvariablen in Karten ausgegeben werden 
(Abb« 7)« Hach Durchlaufen aller Perioden wird die Entwicklung 
der Flächennut zung und des Verkehrsaufkommens der Ge samt stadt 
sowie einer Reihe von Kennziffern zur Bewertung der Flächen- 
nutzungs- und Verkehrs Struktur der Stadt in Tabellen und Dia- 
grammen dargestellt« 

Daneben werden die wichtigsten durch die untersuchte Variante 
verursachten einmaligen und wiederkehrenden Einnahmen und Aus- 
gaben der Stadt und anderer beteiligter Gruppen (Bund/Land/ 
Kreis, Verkehrsbetriebe, Verkehrsbenutzer, Bauherrn, Mieter) 
errechnet und in einer Matrix der Zahlungsströme der einzelnen 
Perioden dargestellt« Abschließend werden die Zahlungsstrom- 
salden der einzelnen Gruppen, auf das Ausgangs jahr bezogen, 
abgezinst und akkumuliert und mit den entsprechenden Salden 
einer hypothetischen "Nullvariante" verglichen« Die Null Vari- 
ante ist definiert als "Entwicklung ohne Planungseingriffe" 
und dient als Vergleichs bas is für alle Planungs Varianten« 

3« Bewertung der Ergebnisse 

Es ist die Aufgabe des POLIS-Mödells , Bas is Informationen für 
die vergleichende Bewertung der Planungsvarianten zu liefern« 
Implizite Wertungsvorgänge im Planungsprozeß sollen duroh lo- 
gisch konsistente und nachvollziehbare Verfahren überlagert 
werden« 

Die Beziehungen zwischen den untersuchten Maßnahmenkombinatio- 
nen und ihren Konsequenzen wird durch das Simulationsmodell 




bergest eilt» Die quantitativ® Erfassung der Auswirkungen einer 
Variante in Simulationsmodell soll ein differenziertes Bild 
nicht nur der teohnisohen v sondern auoh der ökonomischen und 
sozialen Aspekte dieser Auswirkungen ergeben. In der Nutzen- 
bilanz zum Beispiel sind die Auswirkungen des Verkehrssystems 
der untersuchten Variante für verschiedene Gruppen dargestellt 
(Abb. 8)s die Benutzer des Verkehrssystems , die Einwohner ins- 
gesamt sowie Minderheiten wie alte und arme Einwohner, die Be- 
triebe und Institutionen und schließlich die Umwelt. Die Aus- 
wirkungen werden durch Kennziffern aus gedrückt, deren Entwick- 
lung über die Simulationsperioden verfolgt werden kann. So 
wird die Qualität des Verkehrssystems aus der Sicht der Be- 
nutzer durch Merkmale wie Geschwindigkeit, Stauzeit, Warte- 
zeit, Fußweg und Sitzplatzangebot bestimmt. Für die Einwohner, 
Betriebe und Institutionen äußert sie sich in der Zugänglich- 
keit zu ausgewählten Zielaktivitäten. Zugängliohkeitsmaße wer- 
den außerdem gesondert für die Wohnstandorte der alten und ar- 
men Einwohner unter ausschließlicher Berücksichtigung der Ver- 
sorgung mit öffentlichen Nahverkehrsmitteln berechnet. Die ex- 
ternen Auswirkungen des Verkehrssystems auf die Umwelt und die 
damit verbundenen sozialen Kosten sind durch Meßgrößen wie 
Luftverschmutzung, Lärmbelästigung, Behinderung, Unfälle und 
Sachs eliäden erfaßt. 

Die finanziellen Implikationen der Planungsvarianten dienen 
als eine weitere Bewertungsggundlage • In ihnen kommen die Re- 
alisierbarkeit der Variante sowie die Relation zwischen Auf- 
wand und Ergebnis der Planung zum Ausdruck» Durch die diffe- 
renzierte Angabe der Verteilung der Gesamtkosten auf die' ein- 
zelnen Gruppen der Stadt wird die Frage, welchen Interessen 
die Planungsvariante nützt, ln das Bewußtsein gerückt. 

Zwekk des eigentlichen Bewertungsmodells wird es zunächst 
sein, die vielfältigen, verschieden dimensionierten und ska- 
lierten Auswirkungen der Planungsvarlanten vergleichbar dar- 
zustellen. Dem Sntscheldungsträger sollen die unverfälschten 
Implikationen der verschiedenen Varianten, geordnet und ver- 
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gleiohbar gegenübergestellt, snr Interpretation überlassen 
werden« Um jedoch aus der Fülle der dabei anfallenden Baten 
eine eindeutige und nachvollziehbare Rangfolge der Varianten 
hinsichtlich ihres Nutzens ableiten zu können, sollen darü- 
berhinaus die verschiedenen Nutzeneffekte einer Variante in 
einen Gesamtnut zennaß zusammengefaBt werden. Dies bedingt eine 
Gewichtung der Planungsziele sowie die Zuordnung von Nutzen- 
Funktionen zu den einzelnen Kerknalen der Varianten« Die Be- 
wertung soll in den verschiedenen Stadien der Planung wieder- 
holt vorgenommen werden. 
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Abb. 2 Fl&ohttimut sung und Yerkehrssjstea 






Periode (1) 
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Abb* 3 Inrentardaten und ihr# Veränderung 
im Ablauf einer Simul&tionaperiede 




1 WOHHUHGSBAU 2 GETOEBEAHSIEDLDHG 3 GEMEIHBEDAB7SEIHBICHTOHGH I 
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Abb. 4 Planungsmaßnahmen zur Fläohennutzung: Eingabe foraular (Sohena) 











UKUUSTHICKOr (OHTY) UMBitrSäBCEER ( KT) 
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Abb, 5 Verkehrs baum&ßn&hmen : Eingabeformular (Schema) 



Periode E - Länge FQ - Zugfolge 

A-Kneten T/G - Zeit/Geschwindigkeit Z - Zone 

B-Knoten L - Linien BJ - Baujahr 
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Abb. 6 Programmaufbau des POLIS- 
Hodells : 

POLIS 1 (oben) 

POLIS 2 (unten) 



Eingabe Programm Ausgabe 
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Abb. 7 Beispiele für Ausgabe der Sinulationsergebnisse (Köln)s 
Zugänglichkeit, Baulandreserre , Einwohnerdichte, Reise- 
zeiten zu Beginn (links) und an Ende (reohts) der Simu- 
lation 
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Ein Modell zur optimalen Wärmeversorgung in Städten 
durch leitungsgebundene Energieträger 
von K. P. Liesenfeld, Aachen 

1 . Einführung 

Im Rahmen der allgemeinen Bemühungen, die Verunreinigungen 
unserer Umwelt zu verringern, sind Maßnahmen zur Reinhaltung 
der Luft von besonderer Bedeutung. Im Raumwärmesektor bietet 
sich dazu (bei genügend hoher Wärmedichte des betrachteten 
Gebietes) die Verwendung leitungsgebundener Energieträger an, 
da diese einen weitgehend imissions- und emissionsfreien Be- 
trieb von Wärmeanlagen ermöglichen. 

Man nennt Energieträger leitungsgebunden, wenn sie zum Trans- 
port an den Endverbraucher fest installierte Systeme - Netze - 
benötigen: also Gas, Heizwasser, Dampf und elektrischer Strom. 

Ein Energieversorgungsunternehmen (EVU), dem die Versorgung 
eines räumlich abgegrenzten Gebietes - etwa einer Stadt - mit 
den Energieträgern Gas, Heizwasser und elektrischer Strom ob- 
liegen möge, kann nun vor der folgenden Problemstellung stehen 

Das EVU will innerhalb seines Versorgungsgebietes in den der- 
zeitig nicht mit leitungsgebundenen Energieträgern versorgten 
Wärmemarkt expandieren, ohne jedoch einen Teil seiner bishe- 
rigen Abgabe an leitungsgebundenen Energieträgern, der mit 
Grundlast bezeichnet sei, zu beeinträchtigen und ohne sein 
ausgenblicklich bestehendes Energieverteilungs - und Beschaf- 
fungssystem ausbauen zu müssen. Die einzigen Investitionen, 
die es bezüglich seiner Expansion tätigen will, seien die, 
welche zur Erstellung von Hausanschlüssen notwendig sind. 



Das Modell ist Bestandteil eines Forschungsauftrages, der 
am Lehrstuhl für Unternehmensforschung an der RWTH Aachen 
mit Unterstützung der WIBERA A.G. durchgeführt wurde. 
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Der augenblicklich nicht mit leitungsgebundenen Energieträgern 
bediente Teil des Wärmemarktes sei mit Wärmepotential bezeich- 
net. Dieses Wärmepotential ist örtlich qualitativ und quanti- 
tativ bestimmt durch: 

a) eine mögliche Substitution von nicht leitungsgebundenen Ener- 
gieträgern, die einen Teil des bereits bestehenden Wärmebe- 
darfs decken, durch leitungsgebundene, 

b) eine mögliche Bedienung des noch latenten Wärmebedarfs mit 
leitungsgebundenen Energieträgern . 

Das so definierte Wärmepotential soll sich in der Hauptsache 
auf Raumheizungs zwecke beschränken. 

Die vom EVU angestrebte Absatzexpansion soll so erfolgen, daß 
sie unter Einhaltung gewissener Nebenbedingungen , die später auf 
geführt werden, einen maximalen Periodengewinn bei konstanten 
Preisforderungen erbringt. 

Zur Bedienung des Wärmepotentials mit leitungsgebundenen Ener- 
gieträgern sollen Gas, Heizwasser und elektrischer Strom zum 
Schwachlasttarif (Schwachlaststrom) zur Verfügung stehen. 



2 . Skizze des mathematischen Modells 

2.1 Darstellung der räumlichen Verteilung des Wärmepotentials 
im Modell 



Zur Darstellung der räumlichen Verteilung des Wärmepotentials 
wird das Versorgungsgebiet in i(i=l,...,n) Zellen unterteilt. 
Diesen Zellen werden dann die entsprechenden Einzelabnehmer 
und deren Wärmepotentiale zugeordnet. Damit bestimmt sich das 
Wärmepotential a^ b (i=l,...,n; b=l,...,4) einer Zelle i aus der 
Aggregation der Wärmepotentiale aller einer Zelle zugeordneten 
Einzelabnehmer, wobei der Index b(b=l,...,4) zwischen den Ab- 
nehmergruppen, in die die Einzelabnehmer eingeteilt werden kir- 
nen, differenziert. Der Index b kennzeichnet mit: 




b*l: Abnehmer, die Wärme zu Zwecken der Einzelheizung 

b=2 : Abnehmer, die Wärme zu Zwecken der Etagenzentralheizung 

b = 3: Abnehmer, die Wärme zu Zwecken der Kellerzentralheizung 

b*4: Abnehmer, die Wärme zu sonstigen Zwecken 

benötigen* 

Aus Abb. 1 ‘können die Beschränkungen und die Bedienungsmöglich- 
keiten des Wärmepotentials einer Zelle, nach Abnehmergruppen 
getrennt, für die einzelnen Energieträger 1(1 = 1, 2, 3) entnommen 
werden, wenn man dem Schema in Abb. 1 folgende Interpretation 
zugrunde legt : 

M 1 M : Bedarfsdeckung mit Energieträger 1 bei der Verbraucher- 

gruppe b möglich 

”0 M : Bedarfsdeckung nicht möglich. 



Für Wärmebedarfsdeckung 
möglicher leitungsge- 
bundfener Energieträger 


b = l 


l 

1 


Abnehmergruppe 

b=2 1 b = 3 ! 

1 1 


b = 4 


Gas 1=1 


1 


"T 

i 


1 ' 


! 

1 1 


1 


Heizwasser 1=2 


0 


1 


0 1 


1 ' 


1 


Schwachlast- i - q 

ström i-o 


1 


i 

1 


1 

o | 

1 


1 

0 , 
I 


1 


Wärmepotentiale a^ b 


a ü 


I 

l 


a i2 


* . 

a . ~ 
i3 | 


a i4 


in der Zelle i 




i 

— 1— 


1 


1 

L- 





Abb. 1 

Beschränkungen und Bedienungsmöglichkeiten des Wärmepotentials 

einer Zelle 
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Mit Hilfe von Abb. 1 sei nun noch die Matrix R * (r^ 1 ^ 
definiert, mit 



1, falls Bedarfsdeckung mit Energieträger 1 bei 
der Verbrauchergruppe b möglich ^ 

0, sonst. 



2 . 2 Die Zielfunktion des Modells 



Um die Zielfunktion darstellen zu können, werden zunächst die 
Variablen 



( i = l , . . . ,n) 

y ib » y ib ^ 0 (1 = 1, 2, 3) (2) 

(b = l , . . • ,4 ) 



eingeführt. Sie beschreiben die Wärmeleistung, die der Abneh- 
mergruppe b zur Bedienung des in der Zelle i bestehenden Wär- 
mepotentials a^ b mit der Lieferung des Energieträgers 1 zur 

Verfügung gestellt wird. Führt man jetzt noch die Preiskoeffi- 
1 ) 

zienten 

( i = l , . . . ,n) 

(1 = 1, 2, 3) (3) 

(b = l , . . .>) 



ein, so läßt sich die zu maximierende Zielfunktion schreiben 
zu: 



Z <yi£> « £ £ 



i=l b=l 



1 = 1 



1 r l 

ib r b 



'ä 



(4) 



Die Preiskoeffizienten c^ setzen sich aus folgenden Bestand- 
teilen zusammen: 



Bei Gas muß wegen unterschiedlich hoher Erlöse eigentlich 
zwischen den Kundenarten "bisheriger Gaskunde" und "bis- 
heriger Nichtgaskunde" unterschieden werden. Darauf wird 
hier jedoch aus Gründen einer vereinfachenden Darstellung 
verzichtet . 
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a) urbeits- und leistungsabhängigen Erlösen , 

b) arbeitsabhängigen Beschaffungskosten, 

c) erfolgswirksamen Auflösungen von Baukostenzuschüssen, 
die der Kunde bei Errichtung eines Hausanschlusses evtl, 
an das EVU zu zahlen hat , 

d) jährlichen Abschreibungen des EVU's auf die für einen 
Hausanschluß erstellten Anlagen, 

e) Transportkosten der einzelnen Energieträger. 

Auf eine nähere Beschreibung der Art und Weise wie diese Be- 
standteile in den Preiskoeffizienten berücksichtigt werden, 
soll hier verzichtet werden. 



2 . 3 Die Restriktionen des Modells 

2 . 3.1 Bedienung der den Zellen zugeordneten Wärmepotentiale 

Die einzelnen Energieträger 1 (1 = 1, 2, 3) können einander bei der 
Bedienung des in der Zelle i (l=l,...,n) bestehenden Wärmepo- 
tentials a^ b (b=l,...,4) innerhalb des in Abb . 1 dargestellten 
Rahmens substituieren. Nimmt man für die bei den einzelnen Ab- 
nehmergruppen b eingesetzten Energieumwandlungssysteme einen 
Wirkungsgrad an, so läßt sich die folgende Nebenbedingung 

formulieren: 




1 = 1 (b = l * , . . ,4 ) 
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2.3.2 Berücksichtigung der Betriebs zustande der Netze 

Die zum Transport der leitungsgebundenen Energieträger be- 
nutzten Netze sind im allgemeinen vermascht und kapazitäts- 
beschränkt. Da die Abnehmer gewisse Mindestpotentiale (Drücke 
bzw. Spannungen) bezüglich der gelieferten Energieträger for- 
dern und die Potentiale vom jeweiligen Betriebs zustand des 
Netzes abhängig sind, müssen diese Betriebs zustände im Mo- 
dell berücksichtigt werden. Ein Betriebs zustand heiße im fol- 
genden zulässig, wenn dieser Betriebs zustand technisch mög- 
lich ist und außerdem die Mindestpotentialbedingungen er- 
füllt. 



2. 3. 2.1 Berücksichtigung des Betriebs zustandes beim Energie- 
träger 1=3 (Schwachlaststrom) 

Eine nähere Erläuterung des Begriffs Schwachlas tstrom ist mit 
Hilfe von Abb. 2 möglich, die eine qualitativ typische Bela- 
stungskurve darstellt. Diese Belastungskurve ist neben ihrer 
Spitzenlast durch ihre Belastungstäler während sogenannter 



Schwachlastzeiten charakterisiert . 




Belastungskurve bei elektrischem Strom 
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Mit Schwachlaststrom (oder Nacht- und Tagesnachlades trom 
bei Raumheizung) wird nun die elektrische Leistung bezeich- 
net, die man zur Auffüllung der Belastungs täler ( zu in der 
Regel verhältnismäßig niedrigem Preis ), zusätzlich verkauft. 
Die in Abb . 2 schraffierte Fläche beschreibt dann die Schwach- 
lastarbeit . 



Setzt man voraus, daß der Betriebs zustand der elektrischen 
Netze auch während der Verteilung der jährlichen Spitzenlast 
zuläassig ist, dann kann man zur Sicherstellung eines zuläs- 
sigen Betriebs zustandes bei Schwachlaststrom folgendermaßen 
vorgehen : 



Man faßt die Zellen i(i=l,...,n) eines Versorgungsgebietes 

3 

so zu Zellengruppen ;j(j = l,...,J ) zusammen, daß man jeder die- 
ser Zellengruppen eine sekundäre Quelle mit der Quelleistung 

3 

Qj zuordnen kann. Zur Sicherung eines zulässigen Betriebs- 
zustandes genügt dann die folgende Nebenbedingung: 



CO 

1 

■J3 


L — 

b = 1 , 4 


6 

^ib 


z. 




mit 


3 

z j 


- 






i . 3 

3 


= 






j 3 


= 


wobei 


bei der 


Festlegung 


Erwärmung der 


Kabel 


für 


muß . 









( 6 ) 



Grundlast innerhalb der Zellengruppe j , 

Indexmenge für alle Zellen der Zellen- 
gruppe j, 

Anzahl der Zellengruppen, 



2. 3. 2. 2 Berücksichtigung des Betriebs zustandes bei den Energie- 
trägern 1=1,2 (Gas, Heizwasser) 

Bei Gas und Heizwasser sind die für Schwachlaststrom getroffe- 
nen Annahmen nicht gültig. Es müssen deshalb, um die Mindest- 
potentialbedingungen zu erfüllen, die bei den verschiedenen 
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Belastungsfällen an den Knotenpunkten der Netze auftretenden 
Drücke explizit im Modell berechnet werden. Zu diesem Zweck 
müssen Netzberechnungen in das Modell integriert werden. 

Zur Berechnung eines Netzes mit N Kanten und K Knotenpunkten, 
wovon E Einspeisestellen seien, kann man zur Bestimmung der 

N + E + K- E= N + K (7) 

unbekannten Flüsse und Drücke: 

a) K Gleichungen aus den Kirchhof f 1 sehen Knotenpunktsbedingungen 

mit N + E unbekannten Flüssen, 

b) N Gleichungen aus den Druckverlusten der N Einzelrohre des 

Netzes mit K - E unbekannten Drücken, 

c) die Belastung aller Knotenpunkte, 

d) die Rohrnetzdaten, 

e) die Einspeisedrücke 

heranziehen . 

Zur Aufstellung der einer integrierten Netzberechnung entspre- 
chenden Nebenbedingungen wird nun zunächst jedem Knotenpunkt 
eines Netzes eine Zelle i(i=l,...,n) zugeordnet. 

Netze zum Transport von Gas und Heizwasser lassen sich in die 
Klasse der ungerichteten, endlichen und einfachen Graphen ein- 
ordnen. Zur Berechnung der Betriebs zustände müssen auf den 
Kanten der Netze jedoch Bezugsrichtungen zur Orientierung der 
auf tretenden Flüsse gewählt werden. Diesem System von Bezugs- 
richtungen werde nun ein gerichteter Graph mit der Inzidenz- 
matrix 



(ts > Vy 



( 1 = 1 , 2 ) 



( 8 ) 
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zugeordnet, wobei j = Index für alle Kanten, 

s = Index für alle Knotenpunkte, 
K^ = Anzahl der Knotenpunkte im 
Netz des Energieträgers 1, 

N^ = Anzahl der Kanten im Netz des 
Energieträgers 1 . 



Damit führt die Formulierung der Kirchhof f * sehen Knotenpunkts - 
bedingungen auf die Nebenbedingungen: 



4 K (1=1,2), 

11 P 1 T + 1 n 1 _ 1 

r - - E Z_ t = -i <l-i = ~ z i 



*t 

r 



ib *b 



b = l 



j=l 



s: 3 



(9) 



( s = 1 , • . . ,K X ) , 



mit 



HJ- — <D . 



wobei F^ = konstanter Faktor zur Umrechnung der Flüsse 
in Wärmeeinheiten 



= Variable zur Darstellung des in der Kante j 
auftretenden Flusses beim Energieträger 1 



ist, und die Zuordnung 0 — <D definiert sei durch: 

a— ® : : = Zelle i ist Knotenpunkt s (lo) 

zugeordnet . 

Bei E 1 (1=1,2) Einspeisestellen lassen sich die folgenden Ne- 
benbedingungen formulieren: 

P s = P Es ( s= K 1 -E 1 -1, K^'-E 1 -2 , . . .‘jK 1 ) (11) 

( 1 = 1 , 2 ) , 

mit Pg ^ 0, wobei 

P Es = Konstante zur Darstellung des Einspeise- 
druckes am Knotenpunkt s für Energie- 
träger 1, 
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pg = Variable zur Darstellung des Betriebs- 
druckes am Knotenpunkt s für Energie- 
träger 1 . 

Nimmt man nun noch die Druckverlustgleichungen für alle in den 
Netzen auftretenden Einzelrohre in das Modell auf, dann ist 
der Betriebs zustand der Netze vollständig beschrieben. Ohne 
jetzt eine Aussage über die Form dieser Druckverlustgleichungen 
machen zu wollen, können dieselben formal folgendermaßen be- 
schrieben werden durch: 

Ps - Psii - f j <q j> » (12) 

wenn die Kante j eines Netzes mit den Knotenpunkten s und s + 1 
inzident ist. Mit (12) und (8) kann jetzt die folgende Neben- 
bedingung zur Bestimmung der Druckverluste im Einzelrohr for- 
muliert werden: 

Ps - Psii - u] f j ^ ■ 0 (1=1 ’ 2) » , (13) 

(j=l,...,N A ) . 

Um jetzt schließlich einen zulässigen Betriebs zustand sicher- 
zustellen, genügt es, für entsprechende Knotenpunkte die Ein- 
haltung gegebener Mindestpotentiale P m ^ n zu fordern, und die 
gesamte dem System zur Bedienung der Wärmepotentiale verfüg- 
bare Leistung zu beschränken. 



Also müssen noch die Ungleichungen 

Ps> P miis’ K 1 - E 1 

( 1 = 1 , 2 ) 

und 



(14) 



£ £ ^ - £ ** < 1 * 1 - 

i = l b = l u u 1 = 1 

mit u 1 e U 3- und Ju 3 "! = E 3- , 



2) (15) 
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wobei = Indexmenge aller Einspeisestellen 

des Energieträgers 1 

Q , = Konstante zur Darstellung der maxi- 
u 1 

malen Kapazität der Emspeises teile 

1 

u 

in das Modell ist, aufgenommen werden. Damit ist das gesamte 
Modell formal beschrieben. 



2 . 4 Struktur und Klassifikation des Modells 

Die Struktur des Modells ist blockartig, da die einzelnen 
Energieträger in voneinander unabhängigen Netzen transportiert 
werden . 

Weiterhin ist das Modell linear in allen Variablen, mit Aus- 
nahme derer, die zur Berechnung der Druckverluste in Einzel- 
rohren eingeführt wurden. Da diese in Gleichheitsnebenbedin- 
gungen auftreten, kann das Modell als i.a. nichtkonvexe Maxi- 
mierungsaufgabe mit linearer Zielfunktion gekennzeichnet wer- 
den. 

Die Dimension des Modells ergibt sich bei Berücksichtigung 

eines Gasnetzes mit 35o Knotenpunkten u. 5oo Kanten, 

eines Heizwassernetzes mit 6o Knotenpunkten u. 7o Kanten, 

eines Stromnetzes mit 6o Zellengruppen 

zu 3.092 x 4.180, wenn 400 Zellen gebildet werden. 
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3 . Praktische Berechnung des Modells 

Für die im Modell betrachteten Netze sind bezüglich der Ener- 
gieträger Gas und Heizwasser folgende Annahmen gültig 
[3, S. 96-loo ; 4] : 

a) die Strömung ist turbulent 

b) die Rohrwandungen verhalten sich wie im Übergangsbereich 
zwischen hydraulisch glattem und hydraulisch rauhem Ver- 
halten . 

Zur Berechnung des Druckabfalls Ap längs eines Einzelrohres, 
welches von einem Fluß q durchflossen wird, kann dann Gleichung 
(12) in folgender Form angegeben werden [3, S. A2; 4] : 

A p (q) = aq E ^ q) , (16 ) 

worin der Faktor q neben Rohrkonstanten das spezifische Gewicht 
des durchfließenden Mediums beschreibt, und der Exponent g (q) 
zunächst nicht explizit darstellbar ist [4l . Somit ist die 
Form der Druckabfallgleichung gemäß Gleichung (16) für eine 
praktische Berechnung des Modells nicht geeignet. 



3 . 1 Approximation der Druckabfallgleichung für Einzelrohre 

Die Gleichung (16) wird zur Darstellung des Druckverlustes für 
Einzelrohre im Modell linear approximiert. In Anbetracht der 
Dimension des Modells wurde eine Approximation durch lediglich 
zwei Geraden angestrebt. Die Zulässigkeit einer solchen Vor- 
gehensweise konnte durch umfangreiche Testnetzberechnungen 
[3 , S. A9 ffTJ nachgewiesen werden. 

Die Berechnung der approximierenden Geraden kann nach der Me- 
thode der kleinsten Quadrate erfolgen. Mit der Bestimmung der 
beiden Geraden hat man dann einen die Druckabfallgleichung 
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Ap (q) approximierenden Plygonzug f* (q) zur Verfügung, so daß 
das Modell mit den Methoden der Separablen Programmierung be- 
rechnet werden kann. 

Zur Darstellung des approximierenden Polygonzuges f* (q) werden 
die in Abb. 3 dargestellten Bezeichnungen gewählt. Führt man 
nun noch spezielle Variable 

A k , 0 £ £r 1 (k = l 4) (17) 

ein, so lassen sich alle Punkte des Polygonzuges f* (q) im 

it ^ 

Intervall 4? f (q)i= entsprechend der bei der A -Methode 

der Separablen Programmierung üblichen Art und Weise folgender- 







Abb. 3 

Approximierender Polygonzug der Druckabfall- 
gleichung für Einzelrohre und dessen konvexe 
Hülle 



maßen beschreiben: 

■ t . »; 

k = l 



( 18 ) 
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4 

mit \ = 1, A k 0 

k = l 

und a) nicht mehr als zwei positiv, 

b) nur benachbarte positiv . 

Für die Variable q im Intervall 4 q 4- q 4 gilt 
analog : 

4 

q = IZ q k *k * 

k = l 
4 

mit ) J\ k = 1,- A k ^ 0 

k = l 

und den Bedingungen (2o). 



Damit können jetzt die Nebenbedingungen (9) und O 
werden zu: 



t 



b = l 



y. 1 r 1 - 1 
y ib r b 



N 



- f1 L *.ä L ^ ■ -i 

j=i 



k = l 



( 1 = 1 , 2 ) , 

(s=l, . . . .K 1 ) , 



mit 



0 — 0 



und 



Ps Ps 



1 _ + 1 **1 1 ] 
+ 1 t sj Z f jk ^j> 



33 L jk ' v jk 

k = l 



= 0 ( 1 = 1 , 2 ) , 






(19) 

( 20) 

dann 

( 21 ) 
( 22 ) 

3) umgeformt 
> 

(9a) 

(13a) 
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wobei 



II *jj ■ > 0 



k = l 



und die Bedingungen (2o) erfüllt sein müssen. 



3 . 2 Bestimmung des Optimums 

Da die oben formulierte Maximierungsaufgabe im allgemeinen 
nicht konvex ist, wird man im allgemeinen ein lokales Optimum 
als Lösung dieser Maximierungsaufgabe erhalten. 

Es soll jetzt ein Verfahren skizziert werden, das die Bestim- 
mung des globalen Optimums bei den zum Test des Modells ver- 
fügbaren Daten in zwei Schritten ermöglichte. 



Der erste Schritt dieses Verfahrens besteht in der Berechnung 
einer oberen Grenze (Bound) für den Wert des globalen Optimums 
der zu lösenden' Maximierungsaufgabe. Zu diesem Zweck wird für 
jede Teilstrecke j ( j =1 , . . . ,N 1 ; 1 = 1,2) die Menge aller Punkte 
des die Druckabfallgleichung approximierenden Polygonzuges 
f* 1 (qj) durch deren konvexe Hülle (P^) ersetzt. Dazu ge- 
nügt es hier, einfach die Bedingungen (2o), die einen be- 
schränkten Basiseintritt bezüglich der speziellen Variablen 
fordern, fallen zu lassen. Damit erhält man dann eine 
"modifizierte" konvexe Maximierungsaufgabe, deren globales Op- 
timum mit den bei der linearen Programmierung üblichen Methoden 
berechnet werden kann. Der Wert des so bestimmten Optimums ist 
außerdem Bound der ursprünglichen Maximierungsaufgabe . 



Im zweiten Schritt des Verfahrens werden nun zunächst alle 
Flüsse (j=l,...,N^; 1=1,2) mittels der beim ersten Schritt 
erhaltenen Lösungswerte der speziellen Variablen 
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gemäß Gleichung (21) berechnet. Die so berechneten Flüsse 
erfüllen alle Knotenpunktsbedingungen sowohl der "modifizier- 
ten" als auch der ursprünglichen Maximierungsaufgabe, da diese 
Bedingungen linear sind. 

Anschließend werden zur Darstellung der eben berechneten Flüsse 
Werte für die speziellen Variablen so gesucht, daß 

diese die Bedingungen (2o) erfüllen. 

In der beim ersten Schritt gefundenen optimalen Lösung werden 
nun alle Werte der speziellen Variablen dieser Lösung durch 
jk ersetzt. Damit erhält man einen Vektor, der eine (i.a. nicht 
zulässige) Startbasis für die ursprüngliche Maximierungsaufgabe 
bildet . 

Ergibt die erneute Berechnung der Maximierungsaufgabe - unter 
Einhaltung der Bedingungen (2o) - einen Zielfunktionswert, der 
dem im ersten Schritt ermittelten Bound entspricht, so hat man 
das globale Optimum der ursprünglichen Maximierungsaufgabe ge- 
funden . 



4 . Approximat ions fehler 

Zur Sicherstellung eines zulässigen Betriebszustandes sind die 
Druckbestimmungen an den Knotenpunkten mittels der in das Mo- 
dell integrierten Netzberechnungen von ausschlaggebender Be- 
deutung . 

Um die Fehler, die bezüglich der Druckbestimmung an den Knoten- 
punkten der Netze bei der Berechnung des Modells auftreten^ab- 
zuschätzen, sind Netzberechnungen mittels konventioneller Ver- 
fahren Jjl , 2, 4, s] durchgeführt und deren Ergebnisse mit denen 
der in das Modell integrierten Netzberechnungen verglichen 
werden . 
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Dabei ergab sich» daß die Abweichung zwischen den im Modell 
bestimmten Drücken von den konventionell berechneten beim 
Gasnetz (1=1) an 94% aller Knotenpunkte nicht über 3 mm WS 
lagen. Diese Abweichungen sind für praktische Belange jedoch 
unerheblich, da sie innerhalb der Fehlergrenze praktisch 
realisierbarer Messungen liegen. Auch beim Heizwassernetz 
bewegen sich die Approximationsfehler innerhalb des für die 
Praxis maßgeblichen Toleranzbereiches ||, S. 162 ff]] . 



5 . Rechenzeiten 

Das oben erwähnte Zwei-Schritt-Verfahren wurde mit dem übli- 
chen £in-Schritt- Verfahren , d.h. die Bedingungen (2o) sind 
stets bindend, zur Lösung eines Testproblems verglichen. Die 
Dimension des hierbei verwendeten Testmodells, dem eine etwas 
andere Formulierung als die hier skizzierte zugrunde lag, er- 
gab sich zu 7oo Zeilen und 1267 strukturellen Variablen, wo- 
bei die Dichte der Matrix vom benutzten LP-Code zu o,56 be- 
rechnet wurde. (Zur Lösung wurde der LP-Code OPTIMA, Vers. 3.0 
benutzt, der -auf einer Rechenanlage CDC 6600 zum Ablauf ge- 
bracht wurde . ) 

Die Ergebnisse dieses Vergleichs sind in Tabelle 1 dargestellt 
Im Kopf dieser Tabelle ist die Central Processor- (CP-) Zeit, 
die von der Eingabe des zum Ablauf des LP-Codes notwendigen 
Steuerkartensatzes bis zur Ausgabe der Rechenergebnisse benö- 
tigt wurde, mit "gesaunte CP-Zeit" bezeichnet. 

Der Tabelle 1 kann entnommen werden, daß die Startbasis für 
Schritt zwei das Zwei-Schritt- Verfahrens in der Nähe einer 
optimalen zulässigen Lösung der ursprünglichen Maximierungs- 
aufgabe liegt, da in Schritt zwei lediglich 8 Iterationen bis 
zum Erhalten des Optimums benötigt wurden. 
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Verfahren 


Anzahl der 
Iterationen 
bis zur Lö- 
sung 


Gesamte 

CP-Zeit 

C se< 3 


CP-Zeit bis 
zum Aufbau 
der 1. Basis 

[secl 


CP-Zeit zur 
Durchführung 
der Itera- 
tionen 
see 


Ein-Schritt 


1034 


1334j49^ 


149^64" 


1159^104^ 


Zwei-Schritt 










Schritt eins 


1282 


696,33” 


149,154" 


521,648" 


Schritt zwei 


8 


216,319” 


157,292" 


5,244" 


Summe 


1290 


912*649” 


306^446" 


526*892^ 



Tabelle 1 

Rechenzeitvergleich zwischen Ein- und Zwei- 
Schritte Verfahren 



Vergleicht man die Rechenzelten beider Verfahren - es sind die 
gestrichelt gekennzeichneten Werte zu vergleichen -, so ergibt 
sich für das Zwei-Schritt-Verfahren eine wesentlich geringere 
Rechenzeit. 
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Entscheidungen einer Person unter Unsicherheit 
von R. Rabusseau und W. Reich, Bochum und Heidelberg 



Vorbemerkung 

Dieser Vortrag beruht auf einer Diplomarbeit von Herrn Reich und mir in Heidel- 
berg; an der mathematischen Fakultät. Die Diplomarbeit selbst ist dis Ausarbei- 
tung und Ergänzung der von Herrn Dr. W. Böge im Sorrmer-SemBster 1969 gehaltenen 
Vorlesung "Spieltheorie" . 

Ich möchte Sie daher bitten, Herrn Reich und mich als Stellvertreter von Herrn 
Dr. Böge anzusehen, mit dessen Einwilligung wir hier vortragen. 

Es handelt sich um die axiomatische Kennzeichnung konsistenter Strategie-Aus- 
wahlprinzipien und die Herleitung der Existenz einer persönlichen reellen Nutzen - 
funktion und deren Eindeutigkeit bis auf isotone affine Transformationen, sowie 
der Existenz und Eindeutigkeit einer persönlcihen Wahrscheinlichkeitsverteilung. 

Die Arbeit gehört in die BAYESsche Entschei dungs theorie, d.h. sie benutzt sub- 

% 

jektive Wahrscheinlichkeiten und wurde angeregt durch das Buch von L.J. SAVAGE: 
The Foundation? of Statistics. Allerdings werden die Axiome von L.J. SAVAGE teil- 
weise von W. Böge durch andere ersetzt. Anders als bei SAVAGE wurde nicht zuerst 
die Existenz einer persönlichen Wahrscheinlichkeitsverteilung, sondern zunächst 
aus den Axiomen (S K), (R K) und (A T) die Existenz und Eindeutigkeit bis auf 
isotone affine Transformationen eines bezüglich eines BOOLEschen Verbandes der 
Ereignisse additiven reellen Funktionais auf der Menge der Strategien, welches 
mit den Präferenzen verträglich ist, gefolgert. Aus weiterr \ Axiomen, nämlich 
(E) und (Z) wurde dann gefolgert, daB dieses Funktional ei i Integral bezüglich 
elne^ eindeutig bestimmten Wahrscheinlichkeitsmaßes ist. 
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1 . STRATEGIE* AUSWAHLPRINZIPIEN 

Eine Person P befinde sich in einem Entscheidungslonf likt . Wir 
gehen aus von folgenden Grundbegriffen: 



a) der Menge S der externen Schicksale (Natur), 

b) der Menge H der möglichen Handlungen von P, 

c) der Menge C der möglichen Konsequenzen für P, 
die sich ergeben aufgrund der von P ausgewählten 
Handlungen , 

d) die Konsequenzfunktion c : S * H *C 

e) eine Abbildung x : S ►H nennen wir eine Strategie , 

genauer Handlungsstrategie. Sie gibt an, wie P auf 

das externe Schicksal reagiert. 

. S 

Sei X c H die Menge der P zur Verfügung stehenden 
Strategien . 

f) Da jede Handlung h zu einer Konsequenz c(s,h) für die 
Person P führt, liegt es nahe, jetzt auch die sich für 
jede Handlung und für jedes externe Schicksal ergeben- 
de Konsequenz in den Strategien von P zu berücksichtigen. 

Die einer Strategie x« X zugeordnete Abbildung 

x* : S — ► C mit x (s) : - c(s, x(s)) 



heißt Konsequenz.strategie ■ 
x ist so definiert, daß das Diagramm 




kommutativ ist. 
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B 8 i 8 p i b 1 : 



P sei gezwungen, das Haus zu verlassen und habe die möglichen Handlungen 
mit N oder ohne 0 Regenschirm zu gehen. 

h - tu, o } 

Die möglichen externen Schicksale seien Regenwetter R oder trockenes Wetter T 

s ■ [R. T } 



x (1) R— -M x (2) R>— -*M 



5 

X : * H besteht aus folgenden Strategien (Abbildungen): 

R>-— *M 
T/* 

r v m 

• T‘*0 



T— *0 



x (3) R 



x C4) R\ 

T— M3 



Es ergeben sich die Konsequenzen: 



M 



t n 
P t 



x‘ 1J R. — »t 



Tr^ 



(3)^ 

R>-^n 
T •— *t 



t trocken 
n naß 
p peinlich 



x l2) * R 



t 

T •— *p 



x (4) ‘ R 



n 

T i — ►p 
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Zur Beschreibung von P würde genügen, wie P in einer gegebenen 
Situation aus der Menge X ihrer möglichen Strategien diejenigen 
aussucht, die sie tatsächlich anwendet. 

Wir nehmen aus Gründen der mathematischen Einfachheit an, daB P 
in jeder konkreten Situation nur aus einer endlichen Teilmenge Y 
von X auswählen kann. Sei Pot ß X die Menge der endlichen Teil- 
mengen von X. Diese Auswahl denken wir uns gegeben durch eine 
Abbildung (Strategieauswahlprinzip ) 

f : Pot X *Pot X 

e e 

Wir fordern, daß f dem sogenannten Schnitt konsistenzaxiom (SK) 
genügt : 

(i) /\ f M c M 

M « Pot g X 

(ii) /\ f N I 

0/M« Pot X 
e 

(SK) (iii ) A [Nnf M / 0=^»f(N/\M) - N*f M] 

M,N « Pot X 
e 

Sei f die Menge der (SK) -Auswahlprinzipien , ^ die Menge der Quasi- 
anordnungen (auch Rangordnungen) auf (in) X (d.h. vergleichbar ♦ 
transitiv). Dann gilt der 



SATZ 



Die Quasianordnungen einer Menge X und die 
(SK) -Strategieauswahlprinzipien auf X ent- 
sprechen sich umkehrbar eindeutig vermöge 



? ** 



A f ii 

M «Pot _X 



{x € M| 



A y s x} 
y . M 



f »f* *•< 

A y < x^ l.. - fr x « f (x,y } 
x,y « X Def . 



Wir gehen daher im folgenden immer von einer (festen) Quasianord 
nung auf X aus# 
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Bern, : Sei X : ■ {x^ t . x n } sins endliche Menge und £ eine Quasianordnung 

in X. Qann gilt 

2 

ist durch 3/2 n Angaben bestimmt (im Gegensatz zu n ) 

Für die Anzahl der Quasianordnungen S(n) in der Menge X ergibt sich 

n 

S(n) » £ S(n, k) , wobei 

k-1 



S(n, k) 



£ 






ni 



T 



1 £ n^ £ n - k + 1 

Aufgrund der Rekurs ions formal: 

S(n, k) - k [ S(n-1, k-1) + SCn-1, k)] 

wurden die SCn, k) und S(n) mit dem Computer errechnet. Es ergaben sich die 
folgenden Tabellen. 



S(2.k) 


S(3,k) 


S(4,k) 


S(5,k) 


1 


1 


1 


1 


2 


6 


14 


30 




6 


36 


150 






24 


240 
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1 



3 



13 



75 



541 





128 





V“ 


CM 


O 


o 


o 


o 


o 


o 


CD 


CD 


m 






CNI 


03 


CM 


CD 




CD 


a 


CD 


CD 


to 






O 


cn 


cn 


O 


<3- 


CM 


o 


to 


CO 


in 


% 
























O 




T~ 


in 


CD 


m 


cn 


in 


o 


cn 


CO 




x— 






LD 


T— 


CD 


CO 


m 




CM 


CM 


«3- 


K-M 








CO 


x — 


«3- 


CO 


CM 


cn 


to 


CM 


in 










LO 


co 


cn 


O 


to 


cn 


CM 














V— 


CM 


cn 


x— 




O 
























r— 




T— 


o 


o 


O 


O 


o 


CD 


o 


O 




r" 


2 




T— 


in 


CO 


CM 


CM 


co 


(M 


CO 




to 






m 


r* 




v- 


X~ 


*3- 


in 


CO 




CM 


























cn 






03 


to 




cn 


CO 


*— 


CM 












T 


CO 


m 


o 


CM 


cn 


to 




CO 










■c— 


CO 


cn 


cn 


*3- 


cn 




CD 


CO 




































V— 


CM 


T— 












*3- 


CO 


^r 


o 


CD 


CD 


o 






in 






in 


cn 


CM 


o 


CM 


CM 


CM 






tn 






CM 


r\ 


CO 


o 


in 




cn 






CO 








cn 


C3 


to 


*- 


*— 


o 






in 


CO 








*3- 


CM 


cn 










*3- 












T 


T_ 


' r ~ 








in 


CO 
















































2 


V 


to 


to 


o 


CD 


o 


o 








tn 






CM 


cn 


a 


o 


CM 










cn 


* 




x— 


CO 


■3- 


CO 


* — 


CD 








CM 


r\ 
































CO 


to 


in 


in 








r>* 












T— 












^3- 


co 
















































2 


























T — 


CM 


o 


o 


o 


o 










cn 


to 




CO 


«3- 


to 


o 


CM 










co 








cn 


in 


CO 


f\ 










to 


CO 










X— 












«3- 




S (14, k) 



129 



(M 


o 


o 


O 


o 


o 


O 


O 


o 


O 


o 


O 


o 


cn 


CO 


CM 


CO 


O 


CO 


o 


O 


o 


o 


o 


o 


o 


o 


«fr 


co 


03 


CO 


CM 


in 


CM 


«fr 




ID 


CO 


CM 


CO 


CM 


fr 


CO 


cn 






CO 


r- 


CM 


«fr 


rx 


CO 


rx 


CM 


*— 


o 




cn 


o 


o 


CO 


ix 


«fr 


IX 


in 


o 


(O 


cn 


cn 


rx 




rx 


«fr 


o 


r " 


cn 


rx 


cn 


«fr 


cn 


cn 


co 


CM 


O) 




«fr 


cn 


cn 


«fr 


Cn 


«fr 


in 


ID 


cn 


co 


co 


CO 


CM 






«fr 


o 


rx 




co 


cn 


CO 


«fr 


CT) 


in 


rx 


«fr 






CM 


03 


CD 


ID 


* — 


«fr 


in 




rx 


CD 


«“ 


cn 








«f 


in 


CO 


m 


cn 




rx 


CM 


CD 


rx 


T— 










«fr 


«fr 


«fr 


CD 


cn 


cn 


T“ 


CD 


co 


«fr 












(M 


CO 


CO 


rx 


CD 


CD 


in 




CD 
















<- 


CM 


CM 


▼“ 






CD 





«- o 


CD 


o 


o 


o 


O 


CD 


o 


o 


o 


o 


o 


v- 




cn 


in 


CM 


CM 


«fr 


CD 


«fr 


CD 


o 


CD 


a 


a 


CO 




T- 


rx 


rx 


«" 


CD 


cn 


CO 


rx 


o 


CO 


CO 


CO 


cn 


* 


03 


cn 


CD 


o 


«fr 


CD 


in 


in 


o 


o 


«fr 


o 


co 


cn 




CD 


CO 


CM 


«fr 


cn 


in 


CO 


CD 


«fr 


CM 


CM 


«fr 


v- 




in 


IX 


O 


cn 


rx 


cn 


o 


in 


IX 


V“ 


CD 


cn 


v ~ r 






O 


* — 




in 


CM 


CD 


CM 


IX 


CM 


rx 


GO 


CD 






CD 


o 


v” 


o 


cn 


in 


o 


cn 


CD 


CM 


in 










cn 


rX 


CO 


in 


«fr 




CD 


cn 


CM 


CO 












CD 


CO 


CD 


o 


CM 


IX 


rx 


CD 


CD 














CM 


ix 


cn 


«fr 


cn 


cn 




CM 





«- «fr 


CD 


«fr 


o 


CD 


CD 


CD 


CD 


o 


a 


CD 


1 




cn 


in 


CM 


CD 


«fr 


«fr 


«fr 


CD 


CD 


o 


CD 


in 




o 




CD 


CD 


«fr 


CO 


CD 


O 


o 


CO 


CD 


cn 


* 
























in 


CM 


«fr 


cn 


CD 


CO 


cn 


in 


CO 


CM 


«fr 


CO 


^ — 




«r- 






rx 


CM 


CM 


rx 


CD 


in 


CD 


CD 


o 


IX 


w 




in 


CD 


in 


o 


CD 


cn 


T— 


V" 


m 


CD 


CD 


























in 


CD 






«fr 


in 


cn 


CM 


x— 


cn 


IX 


«fr 


cn 












CD 


in 


CD 


x — 


GO 


03 


cn 


IX 


*-* 










«— 


cn 




«fr 


CD 




CD 


«fr 


cn 



m CD CO ID (M 





«- CD 


CD 


o 


a 


o 


CD 


O 


CD 


CD 


o 


cn 




«fr 


CD 


o 


CD 


«fr 


CO 


CD 


CD 


CD 


a 


IX 




O 


CD 


o 


CD 


CD 


«fr 


CD 


«fr 


a 


co 


in 


«r- 


CM 


r - 


CO 


rx 


CD 


«fr 


T— 


CD 


«fr 


CD 


CM 


T— 




IX 


cn 


co 


m 


cn 


o 


CM 


CO 




cn 






«~ 


«fr 


CD 


CM 


«fr 


o 


r— 


in 


cn 


CD 


CD 






cn 


R 


cn 


CM 


cn 


cn 


cn 


cn 


CM 








CM 


CM 


rx 




cn 


cn 


CM 












«P- 


cn 


«fr 


5 






CD 




32 766 65 534 131 070 



130 



o 


o 


CD 


o 


o 


O 


o 


o 


o 


O 


O 


o 


O 


CD 


CD 


in 


CO 


CM 


CO 


o 


CO 


o- 


o 


o 


O 


o 


o 


CD 


CD 


CD 


CD 


CD 


T— 


cn 


O' 


CD 


CO 


o 


o- 


O - 


O' 


o 


O 


CD 


CD 


CD 


O' 


CO 


CO 


x~ 


CD 


CM 


CO 


o 


CO 


CM 


o 


O 


CD 


CO 


O' 


cr> 


CD 


(N 


CM 


CM 


x — 


Csl 


CO 


x— 


CD 


CO 


CO 


CO 


cn 


ix 


o 


x — 


CD 


x— 


in 


m 


in 


CO 


rx 


IX 


CO 


IX 


IX 


o 


CO 


o- 


o 


CD 


CO 


CD 


ix 


ro 


CO 


x— 


IX 


o 


rx 


o- 


CO 


CM 


CO 


ro 


ro 


O- 


x — 


ro 


CO 


'It 


CO 


CO 


IX 


in 


CM 


CM 




CO 


ro 


xp 


CO 


cn 


o 


CD 


CD 


O' 


in 


CD 


in 


't 


O' 




CO 


CO 


o- 


CO 


CM 


CD 


CM 


in 


CM 


o 


o- 


ro 


CD 


rx 




x— 


rx 


CD 


CO 


ix 


X“ 


o- 


CO 


CD 


co 


in 


X — 


CD 


CD 






CO 


CO 


in 


x- 


CO 


o 


CO 


»X 


IX 


CD 


CD 


O' 


CO 








CM 


CD 


m 


rx 


CD 


CD 


CO 


in 


in 


CO 


in 


in 








T” 


CM 


CM 


in 


o 


CO 


O 


IX 


00 


CM 


o- 


in 










x— 


00 


o 


o 


O' 


IX 


in 


rx 


CM 


CO 


CO 














CO 


o 


O 


CD 


o 




CD 


CM 


















X— 


CM 


(M 


CO 


CM 


x— 







CO 


O' 


o 


a 


CD 


a 


o 


a 


o 


O 


O 


O 


a 


a 


x— 


CM 


o 


CO 


CO 


CO 


o 


CD 


o 


CD 


o 


CD 


a 


CD 


x~ 


CO 


CD 


rx 


T ~ 


CD 


CD 


CD 


o 


CO 


o 


CD 


a 


CD 


CD 


CO 


CD 


x~ 


O 


CO 


CD 


CM 


o 


CO 


CO 


*t 


o 


CO 


in 


IX 


CO 


CM 


CM 


x— 


o 


x— 


in 


CO 


CD 


CM 


to 


CO 


co 




CO 




CO 


CO 


(O 


rx 


CM 


CO 


CO 


O 


«“ 


co 


CM 


CO 


(M 


rx 


CO 


CO 


CO 


rx 


CD 


CO 


o- 


CO 


o- 


CD 


*3- 


(M 


O' 


o 


CO 


CM 


CO 


a 


CO 


CO 


CO 


x — 


CM 


co 




▼“ 


in 


x— 


CO 


CD 


x — 


X — 


rx 


CD 


*” 


CO 


CD 


IX 




O 


x— 


CD 


o 


in 


to 


CO 


CO 


CO 


CD 


co 


CD 


CM 






CO 


CO 


O' 


in 


o 


CD 


CO 


CO 


x— 


X— 


CM 


CM 








OD 


in 


ro 


CO 


CD 


a 


CO 


CO 


CO 


CD 


CD 








x— 


co 


co 


IX 


cn 


in 


CD 


x— 


O' 


CO 


CD 










x — 


to 


CD 


a 


in 


CM 


CO 


CM 


in 


CM 














CM 


rx 


X— 


CO 


a 


in 


x— 




















X- 


x— 


X- 









CO 


O 


CD 


o 


CD 


o 


CD 


O 


o 


o 


o 


O 


CD 


o 


CO 


CD 


to 


CM 


CD 


CD 


CD 


o 


CD 


o 


O 


CD 


CO 


CO 


O 


in 


CO 


CO 


CM 


CD 


o 


CD 


a 


o 


CD 


CD 


(M 


CD 


CO 


CO 


CO 


r - 


CD 


CD 


CM 


CD 


CO 


CO 


in 


O' 


CO 


in 


x— 


CD 


in 


CM 


co 


in 


O' 


rx 


ID 


CM 


in 


CD 


x~ 


CO 


in 


O' 


CO 


o 


to 


CO 


in 


CO 


«5t 


co 


CM 


CD 


CO 


o 


to 


O 


rx 


CD 


O' 


a 


O' 


X— 


x — 


CD 


CD 


in 


co 


CD 


CM 


CM 


CD 


o 


CM 


IX 




O 


CM 


CO 


a 


O' 


in 


CM 


o 


CO 


CD 


rx 


CO 






in 


CM 


o 


O' 


cn 


O 


to 


CD 


CM 


co 


IX 






CM 


O 


to 


CO 


in 


in 


in 


CD 


CO 


in 


8 








CO 


o 


rx 


co 


CD 


a 


CD 


CO 


X“ 










CM 


co 


O' 


in 


CD 


CD 


CO 


CD 


x— 














CM 


O' 


in 


in 


CM 







230 283 190 977 053 5 315 654 681 981 355 130 370 767 029 135 901 



262 142 524 286 1 048 574 

386 634 060 1 160 688 606 3 483 638 676 



131 



xp 


o 


o 


o 


o 


o 


o 


o 


C3 


(D 


CO 


o 


o 


a 


o 


o 


CD 


in 


(M 


o 


co 


CO 


CO 


o 


o 


o 


O 


O 


a 


o 


CD 


CD 


o 


o 


CD 


X— 


CO 


o 


xp 


o 


(VI 


o 


o 


CO 


CO 


O 


o 


CD 


CD 


o 


o 


o 


CD 




X** 


o 


CD 


CO 


Ps 


CO 


xp 


xt 


in 


xp 


CO 


o 


CM 


o 


a 


o 


CD 


CO 


CO 


CO 


IS 


CO 


T— 


o 


CO 


CO 


(VJ 


o 


CD 


o 


in 


CO 


o 


co 


xp 


CD 


is 


co 


CO 


o» 


(VI 


CO 


o 


Ps 


xp 


CD 


ID 


CO 


X— 


CO 


Xt 


o 


CO 


(M 


o 


o 


xp 


o 


CD 


Ps 


o 


CO 


(M 


in 


IS 


in 


CO 


CO 


CO 


co 


CO 


xp 


ps. 


CO 


o 


CD 


in 


x— 


XP 


o 


o 


CD 


CD 


xp 


in 


Xp 


co 


IS 


Ps 


CO 


CD 


ps 


CO 


CO 


CD 


X— 


Ps 


CO 


CD 


CO 


CO 


CD 


CO 




(VI 


ID 


x— 


cn 


in 


xp 


CO 


(VI 


xp 


CO 


CD 


CD 


O 


xp 


o 


rs 


CO 


CO 


rs 


Ps 


CO 


xtf- 


CO 


o 


Ps 


X— 


o 


xp 


x— 


CO 


CM 


X“ 


CO 


CO 


CO 


CO 


cn 


Ps 


C3 


Xtf- 


o 


CD 


CO 


in 


o 


(VI 


CO 


o 


O 


xt 


X” 


rs 


CO 


x— 


co 


CD 


CD 


CM 


v- 


O) 


o 


CO 


(VI 


xp 


CM 


Ps 


in 


CO 


CO 


CD 


Xp 


in 


IS 


CD 


(VI 


CD 




CO 


o 


CO 


CD 


in 


CO 


xt 


(M 


X“ 


CM 


o 




in 


CD 


CO 


'CD 


(D 








X“ 


CO 


CO 


Ps 


in 


CO 


CO 


CO 


CO 


CO 


IS 


CM 


in 


O) 


IS 






CO 


CO 


r- 


03 


CO 


in 


rs 


in 


CO 


CD 


xp 


in 


rs 


CM 


CM 


IS 








in 


X~ 


in 


Ps 


Ps 


rs 


rs 


xp 


CD 


xp 


CO 


in 


x— 


CO 


CO 










CO 


CO 


xp 


CO 


03 


(VI 


xp 


xp 


CD 


CO 


co 


\— 


Xp 


co 












xp 


x— 


in 


CO 


o 


CD 


X- 


Ps 


CO 


o 


CO 


CM 


Ps 














(VI 


Ps 


CD 


CO 


xp 


cn 


CO 


CO 


o 


CM 




IS 


















X— 


co 


in 


in 


XT 


(M 


X“ 






CO 




































CM 


O 


o 


o 


o 


o 


o 


cT 


o 


o 


o 


o 


o 


o 


o 


CD 


CD 




CO 


o 


o 


co 


CO 


o 


o 


o 


o 


o 


CD 


o 


o 


o 


CD 


o 


CD 




CO 


CO 


CM 


co 


in 


CO 


xf 


o 


co 


o 


CD 


o 


o 


o 


CD 


O 


CD 




x— 


CO 


in 


o 


CO 


CO 


03 


CO 


co 


o 


co 


CO 


xp 


CO 


CM 


co 


CM 






o 


CO 


CO 


03 


T~ 


CD 


o 


CO 


CD 


CM 


CO 


o 


CO 




co 


ro 




x— 


o 


*— 


CM 


in 


CO 


CD 


o 


rs 


o 


in 


rs 


CO 


rs 


CO 


xp 


co 




xp 


CM 


1 


CO 


CD 


o 


xt 


CD 


CD 


(VJ 


XT 


CO 


(O 


CD 


in 


CD 


co 




CO 


co 




in 


in 


to 


fs 


CD 




(M 


in 




xp 


CD 


(VI 


rs 


CD 




(M 


(VI 


is 


P'S 


T ~ 


xp 


rs 


X“ 


CD 


CO 


CD 


X- 


rs 


xp 


Ps 


co 


Xt 




CO 


o 


o 


CO 


o 


o 


in* 


(O 


X— 


CD 


(M 


o 


in 


o 


X“ 


CO 


o 




CD 


rs 


x— 


X“ 


xt 


CO 


in 


Ps 


CO 


x— 


x— 


xp 


Ps 


rs 


ro 


C3 


o 




x— 


(VI 


is 


o 


CO 


X— 


Ps 


xt 


in 


(M 


CO 


CO 


co 


Ps 


in 


CD 






CO 




Ps 


CD 


«t 


Ps 


Ps 


in 


X“ 


rs 


CO 


CM 


CO 


in 


xp 


in 


in 




rs 




X“ 


CD 


(VI 


CO 


in 


X“ 


CO 


o 


Ps 


rs 


X“ 


in 


01 


o 






co 






xt 


in 


CD 


co 


O 


Ps 


Ps 


(M 


(O 


ID 


x— 


in 


co 


(O 




in 








Ps 


CO 


in 


(VI 


in 


o 


X- 


xt 


ao 


CD 


o 


xp 


T _ 














CO 


v 


CO 


co 


«t 


x— 


cn 


CD 


Psl 


co 


CD 


CM 




o 












xt 


rs. 


T_ 


(N 


o 


CM 


X“ 


o 


Xp 


o 


X“ 




CO 














X“ 


in 


X“ 


03 


X“ 


CO 


X— 


Xp 


x- 






CM 


















X— 




CM 


x— 


X“ 










CD 


o 


o 


o 


o 


o 


o 


O 


C3 


O 


CD 


C3 


o 


CD 


o 


o 






CO 


xp 


o 


CO 


o 


p 


o 


o 


o 


O 


O 


o 


a 


o 


o 


a 






(M 


CO 


xp 


o 


o 


(VI 


xp 


xp 


xp 


CO 


«t 


<o 


o 


o 


o 


o 






CO 


is 


in 


CO 


co 


CO 


xp 


CO 


xp 


in 


o 


CO 


o 


co 


co 


CO 






in 


co 


CO 




co 


o 


Ps 


(M 


in 


03 


IS 


in 


o 


Xt 


CO 


(VI 






xt 


to 


<p 


rs 


co 


(VI 


in 


CO 


co 




«t 


X- 


ID 


«t 


CD 


rs 






5 


X“ 




o 


co 


x— 


X“ 


CO 


T” 


o 


CM 


o 


CD 


o 


co 


in 






co 


rs 


Ps. 


(VI 


X~ 




X“ 


O 


X“ 


xp 


O 


(VI 


(VI 


CM 


CD 


o 






in 


x— 


O) 


CO 


CD 


X” 


CD 


O 


rs 


CD 


CO 


X“ 


CD 


IS 




rs 






(M 


rs 


«t 


xt 


O 


xt 


Ps 


O 


xp 


CO 


CO 


IS 


ps 


rs 


co 


co 






CO 


CO 


Ps 


cn 


x— 


cn 


CM 


(VI 


co 


CO 


xt 


CD 


CM 


o 


rs 


Ps 






co 




*t 


CO 


(VI 


CD 


CD 


CO 


o 


in 


xp 


CO 


o 


CM 


X“ 


co 






CO 




CO 


CD 


in 


X” 


CO 


(VI 


CM 


co 


CD 


(VI 


co 


xp 


o 


(VI 






xt 






Ps 


CD 


(VI 


(VI 


Ps 


xp 


CO 


CD 


CO 


rs 


xp 


CM 


CD 






cn 








CD 


CO 


in 


CO 


CD 


rs 


CO 


C3 


X— 




Xp 


xp 






in 










Ps 


£ 


xp 


xt 




CO 


CO 


CD 


CD 




CO 






in 












cn 


cn 


CO 


o 


CO 


CO 


CO 


a 


in 








co 
















CO 


CO 


rs 


IS 




CM 










co 



co 




132 



Def . : A £ S KeiBt Zerlegungsmenge für X, falls eins der äqui 

valenten Bedingungen gilt: 



(i) 


/\ X, « 

x.y «X 


y S\A € X 




(ii) 


X A* X S\A - 


X 




( iii ) 


X A* X S\A 


X 




Erklärung: x : S 


i 

< 

< 

X 

X 

1 


— »H 



x A (a) - x(a) 

x A ist die Restriktion von x auf A. 

A 

X A ! " { X A* x€X } ' (x A* y S\A J A 8 ‘ X A 

(x A* y S\A } S\A 8 ‘ y S\A 

SATZ Das System Jr (X) : * {AfSi X ft * x g\/\ “ x } der Zerlegungs- 
mangen von X ist eine BOOLEsche Algebra, d.h. 

A, B €&X) B€#(X) , S«,£(X) 

A«£(X)^^S\A«.£(X) 

Sei£- eine BOOLEsche UnteralgsJjra vonJJ(X). 
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2. DAS RESTRIKTIONSKONSISTENZ-AXIOM (RK) 

Das Axiom (RK) lautet: 

Für alle A €&, alle u, u’€ X^, alle v, v 1 €. Xg^ gilt: 

(RK) 

u * V u ' X v 4 "'"fr U X V ' U ' X V * 



Def . : Für A, B eJS 



A ^ B 4 T7 =t> *A* X S\A * yß" *S\B 
x— *y Def. 



Wir sagen: A ist nicht schlechter als B für den Über- 
gang von x nach y. 



Bern. : ist eine Quasianordnung auf , die die folgende 

x—*y 

Bedingung der Addit ivität erfüllt: 



A Jb B , C 5b D 
x — *y x — *y 

A o C = B ^ 0 * 0 



A o C 5^ B kj D 

x — *y 



Def. : Ein System {A^, ..., A n } von Teilen von S heißt Zerlegung 

von S , falls gilt: 

(i) \!j A. - S 

1 1 

(ii) A A, n a. « 0 
ii'J J 

Wir zeigen anhand eines Beispieles, daß nicht jede additive Quasi- 
anordnung ^ auf & durch ein Maß u erzeugt wird, d.h. 

x— *y 

A A £ B4==b.*'(A) > ^ CB) 

A, Bt£ x— y 

Man muB allerdings für S eine mindestens 5-elementige Hengs wählen 1 ^ 



1) vgl. C. KRAFT, J. PRATT, A. SEIDENBERG: Intuitive Probability on 
Finite Sets. Annals of Math. Statistics 1959, 30, S. 406-419, 

A. KIRSCH: Gerechte lineare Ordnungen und Punktbewertungen, 
in: Der Mathematikunterricht 1969/1, S. 64-64. 
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Aus (RK) folgt das 

Vereinigungslemma (VL) : Seien x,ycX, , ..., A n , B^ , ..., 

B n £ 'ir mit A^ ^ B^ , i * 1 , ..., n. A.^ Aj * B^aBj * 0 

für i J* j. Dann gilt: 

1.) UA. £ UB. 

1 x — ►y 1 



2.) Es gebe ein i€ {1, ..., n} mit A. > B.. Dann gilt 

x — *y 1 

n & 

UA, > ÜB, 



1 r T 

1 1 x— *y i 1 



Bern.: (RK) 4~ -t> (VL) 



Divisionslemma (DL) : Seien m, n c N , x, y e X 



0 = A q , A,j , A m , A^n Aj - 0 für i / j 



0 * B q , B 1 , . .., B n , B^n Bj * 0 für i i* j 

seien Mengen aus & mit 

0 B A £ A, < ... £ A 

o ~ 1 ^ ^ m 

x — *y x — ►y x — *y 



0 - B <£, B, < ... < B 

o ^ 1 ~ ^ n 



x — *y x — ►y 



x — ey 



und 



U A. 4 U B. 

1 x — *y 1 



Dann gilt: 



m 

A A 

r» 0 s 



A . r r ü A, 6 u „Bi 1 

-0 1 1 x— ►y i-n-s*1 J 
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3. DAS ARCHIMEDISCHE UND TEILUNGSAXIOM (AT) 



9 ••• 9 



Axiom (AT) Zu x.y e X, A*„£ , A„>.,0 gibt ss eine Zerlegung 
——————— x-"*y 

[A v . .., A jQ von S, n € N f , mit A^ A für i - 1 
' n x-*y 

Hit Hilfe der bis jetzt formulierten Axiome kommt man fast bis 

zum ersten Hauptsatz. Nur an einer Stelle benötigt man das folgende 

stärkere 

Axiom CAT) Zu A 0, Bf/ gibt es eine Zerlegung 
x— *y 

[ B 1 , . . . , B ] £ £ von B mit B. ^ A für i * 1 , . . . , n. 
n x-*y 

rv 

Für B : - S erhält man CAT). Wir vermuten, daß aus CAT) auch CAT) 
folgt, sind aber bis jetzt nicht in der Lage, dies zu beweisen. 

Def . : Das Tripel CX, ,£ ) heißt (SB)-Raum ( SAVAGE -Bö GE- Raum ) , falls > 

Quasianordnung auf X ist, die den Axiomen CRK) und CAT) genügt. 



In (AT) bzw. CAT) werden für den Fall, daß es x^ y gibt, zum ersten Mal Exi- 
stenzen gefordert. 

Das Axiom CAT) mag zunächst recht künstlich aussehen. Denn während (SK.) und CRK) 
notwendig für die Existenz einpr Nutzenfunktion sind, ist das bei CAT) nicht der 
Fall. Wir zeigen in der ARbeit, daß das Axiom CAT) im Fall des erweiterten Re- 
genschirmbeispieles (Erweiterung durch Würfelwürfe) realisiert werden kann, und 
zwar auf realistische Weise. Dr. Böge stellte die Thesaäuf, daß diese Erweite- 
rung auch. tin jedem anderen Beispiel möglich ist. 
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4. DER 1. HAUPTSATZ 

OBf . : Eine Abbildung IM : X — iR heißt ^r- -additiv , wenn ea eine Familie 

(N a ) V x rR gibt , so daß für jede Zerlegung . 

. .., A n J C £_ von S gilt: 
n 

N(x) - T N Ä (x. ) V x e X 

iTi M i A i 

Dsf . : Eine Abbildung N : X — *IR heißt (person 1. ) Nutzenfunktion von P, falls 

1. N ist £ -additiv 

2. y x > y N(x) ^ N(y ) 
x#y c a 

Wir geben nun 3 Formulierungen des 1. Hauptsatzes. 

1 . Hauptsatz : Sei (X, > y£~) ein (SB)-Raum. Dann gilt: 



(i) Es gibt eine Nutzenfunktion N auf X 

C ii ) Falls S ? 0 ist, so gilt: 

N ist eine weiters Nutzenfunktion auf X 0-£> 

a V 0 N - a-N+ b 
be R 

Zusatz: In (ii) gilt N *=>£.” auch ohne die Vor. S / 0. 



1 . Hauptsatz 



Sei (X, Jb & ) ein (SB)-Raum. Dann gilt: 



(i) 

(ii) 



Es gibt eine Nutzenfunktion N auf X 

N ist eine weitere Nutzenfunktion auf X mit -additiver 

iss 

Zerlegung 

V V N. - a-N. + ß C A) 

a > 0 endl . additives M 
Maß 
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Def » : Eine Abbildung N : X — NR heißt eine Pseudonutzenfunktion auf X, falls gilt 

(i) N ist X -additiv 

(li) £ x t* H M - N ( y J 3 

1. Hauptsatz : Sei (X, > ,<£ ) ein (SB)-Raum. 

(i) Dann gibt es eine Nutzenfunktion N auf X 

Cii) Falls S ^ 0 ist, so gilt: 



N ist eine Pseudonutzenfunktion auf X 6^ i> 

V 

a*Q 
b «R 



N * a • N -i- b 
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5. DER 2. HAUPTSATZ 



Def. : Dis Konsequenzfunktion c : S* H — *C heißt zulässig (für t ), falls 

eine der beiden äqu. Bed. gilt: 

(i) x * y* — O x ~ y 

* 

(ii) Es gibt eine Quasianordnung >,* auf X mit 
x > y crt> x* £ y* 

B e heißt präferenzirrelevant für Y £ X^, wenn für alle 
x,y « Y gilt: 

x y —O x ß y 0 

B € heiße vollständig präferenzirrelevant für ¥ £ X^, wenn es 
präferenzirrelevant für eine geeignete Erweiterung Y' let (Y € Y* £ X^) 
die bezüglich eines geeigneten BQOLEschen Verbandes £ * c 
das Axiom (AT) erfüllt. 

Die Konsequenzfunktion c : S y H C heißt erschöpfend » wenn sie zu- 
lässig ist und wenn jedes B € £ ^ vollständig präferenzi irrelevant ist 
für jede endliche Menge konstanter Konsequenzstrategien in X^» für 
alle A * £ 

Axiom (E) Die Konsequenzfunktion c sei erschöpfend. 

Wir formulieren nun das letzte Axiom (Z). In ihm werden Existenzen gefordert. 

Def. : Ein System £ von Teilen t einer gegebenen Menge heißt abstrakter simpli- 

zialer Komplex , falls gilt: 

(i) Cr [ t>s *-+*'*£] 

Jedes t *£ ist ein endl. Teil der gegebenen Menge. 






Def.: 



Bef.: 



Def.: 



(ii) 
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Beispiele von simplizialen Komplexen: 



r, . 


- {0, 11} } 






C 2 • 


* U. {1}, {2}} 




© (D 


^3 ! 


■ Pot {l, 2} 




© © 


•Q ! 


- CPot {1, 2, 3} ) 


\ { 1 . 2 . 3 } 




*^5 ! 


- Pot {l, 2, 3} 






Zusammenhang 


saxiom (£) : 




Es gibt eine 


Zerlegung{D 1 , . . 


. , D }c / von S und 
mJ 


1 eine Überdeckui 


C * C 4| V e e e 


\j C m mit C i / 0 , 


i * 1 , • • • , m, so 


daß 


’• \ 


m Menge allerJJL 
U i 


-Treppenfunktionen 


mit Werten in C 


2 ‘ C ijK 


: » C ^ a Cj a 


t 0 =-t> 





3 T. t’« C ijk mit Y Dv < # d v Y^ . V «{i.J.k} 
3 e Der abstrakte simpliziale Komplex ^"mit 



g 



- {tg{i, .... m), r\ C. 4 0} 
ist 1 

ist zusammenhängend und einfach zusammenhängend. 



2 e Hauptsatz : Sei}£: * (X,ifc ,£) ein (SB)-Raum, der auch die 

Axiome (E) und (Z) erfüllt. Dann gibt es eine (eben 
falls mit N bezeichnete) Nutzenfunktion 

* 

N ? X ►F , die zu der Nutzenfunktion 

N s X HR in der Beziehung steht: 

N (x) » N (x*) V x« X 



Weiter gibt es ein endlich additives subjektives 
Wahrschein lie hkeitsmaB 

P : — -» [0 , l] und eine subjektive Nutzenfunktion 



F , so daB 
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N ■ j N*dP , d.h. 
a n 

/\ N (x*- ) ■ £ N (c. ) P (A. ) , wobai 
X" 1 1 1 

(A 1# . .., A n )£jC irgendeine Zerlegung von S ist mit 

x* (A^) £ { Cj j (solche Zerlegungen gibt es auch). 

Das Auswahlprinzip f der Person ist bestimmt durch 

f (Y)-fxHi A N (x* ) > N (y* ) für Y€ Pot X. 1 
c yeY 8 

*■ 

Durch f ist P eindeutig und N eindeutig bis auf einen posi 
tiven konstanten Faktor und eine additive Konstante bestinrmt 



d. h. sei N ■ J* N* dP. Dann gilt 

N* • j* N* ' dP .mit N' : » Menge der Nutzenfunktion auf 
- (X,* ,X) genau dann, falls 



V N* • • a • N*+- b 

a > 0 
b c IR 



N (x^ ) ist also gleich dem Erwartungswert der subjektiven Nutzenfunktion bezoger 
auf die Konsequenzstrategie ** bezüglich des subjektiven Wahrschein lichkeitsmas- 
ses P. Wir nennen diesen Erwartungswert auch den mittleren persönl iehen Nutzen 

4h 

von x . 

Der 2. Hauptsatz besagt also: Der Entscheidende wählt aus' jeder endlichen Menge 
von Handlungsstrategien diejenige aus, deren mittlerer persönlicher Nutzen maxime 
ist CBernoulliprinzip) . 

Er handelt so, als würden die externen Schicksale einer Menge A mit der obje 
tiveren Wahrscheinlichkeit f?(A) ein treten, obgleich völlige Unsicherheit hinsicht 
lieh des Eintretens der externen Schicksale besteht. 

Der Entscheidende ist also in der Lage, aufgrund seiner Reflexionen über sein ei- 
genes Auswahlprinzip f eine subjektive Wahrscheinlichkeit P(A), die mit f in Ein- 
klang steht, für alle Mengen A festzulegen. 
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6. PRAKTISCHE BEDEUTUNG OER ARBEIT FOR OIE ENTSCHEIDUNGSTHEGRIE 



Das Ziel der Arbeit let die Herleitung des BERNOULLIschen Prin- 
zips für subjektive Wahrscheinlichkeiten. 

Darüber hinaus möchte die Arbeit nicht nur dem Leser ein Krite- 
rium in die Hand geben« um Entscheidungsprobleme zu lösen« son- 
dern auch ein Hilfsmittel« um einer zweiten Person bei der Fest- 
legung der Präferenzordnung ihrer Strategien zu helfen. Damit 
werden wir in der Lage sein, für eine zweite Person, die die 
Strategien anwenden soll, die notwendigen Entscheidungen zu tref- 
fen. (Das benötigt man in der typischen beratenden Situation des 
Mathematikers im Beruf.) 

Um die Präferenzordnung der Strategien angeben zu können, brau- 
chen wir von der Person Angaben. Im Abschnitt 7 werden wir an 
einem Beispiel sehen, wie stark die Anzahl der Angaben zurück- 
geht. wenn wir die Ergebnisse der Arbeit benutzen. Mit dem 
2. Hauptsatz, dem BtRNQULLI-Prinzip für subjektive Wahrschein- 
lichkeiten. haben wir dann ein Hilfsmittel, das uns die Beschrei- 
bung der Präferenzordnung mit relativ wenig Angaben ermöglicht. 

Zu dem in der Literatur weit verbreiteten noch weniger Angaben 
benötigenden Minimax-Prinzip werden wir ein Beispiel geben, das 
zu einem nicht plausiblen Resultat führen wird. Mit dem 2. Haupt- 
satz jedoch werden wir im selben Beispiel ein Resultat erhalten, 
das der Leser auch ohne mathematische Überlegungen erwartet hätte. 

Somit ist die Arbeit mit dem Minimax-Prinzip nicht allgemein ver 
träglich und möchte auch keine Modifikation desselben sein. 
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Bayesian Reliability Demonstration 

by G. J. Schick and T. M. Dmas, Los Angeles and Culver City 



ABSTRACT 



A Bayesian approach to reliability demonstration testing is 
described and differences between the Bayesian viewpoint and 
the commonly employed classical approach are highlighted. 
A procedure for selecting a specific inverted gamma probability 
density to characterize the prior distribution of the MTBF of 
electronic hardware is developed and a table of Bayesian 
demonstration plans for a practical range of input parameters is 
provided. In addition, procedures for implementation of the 
plans and two illustrative examples are given. Finally, two 
commonly employed classical plans are compared to a Bayesian 
plan illustrating the efficiency of the latter in terms of 
demonstration test time requirements. 
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INTRODUCTION 

Recent publications [1, 2, 3, 4, 5, 6, 10] 1 have shown that the aerospace industry is becoming 
increasingly interested in the advantages generated when reliability demonstration tests of electronic 
systems are based on Bayesian concepts as opposed to classical concepts. In particular, the major 
advantages include: 

A. Less test time, generally, is required to demonstrate a reliability requirement, thus, at 
times, reducing costs dramatically. 

B. The Bayesian approach takes advantage of prior knowledge. 

C. The Bayesian interpretation of test results appears to be more appealing to engineers. 

These advantages have been discussed in considerable detail in the reference literature. 

The objective of this paper is to describe a Bayesian approach to reliability demonstration, provide a 
table of Bayesian based demonstration plans, and specify a procedure for implementing the plans. 

From the statistical viewpoint, reliability demonstration tests are sampling schemes wherein the 
hardware or system is tested by simulating mission conditions in order to determine in a 
probabilistic sense if the reliability requirement has been achieved. Thereafter, the resulting test 
data (sample) and preselected decision rule permit an accept or reject decision. Prior to the 
consideration of the Bayesian approach, reliability demonstration test decision rules have been 
largely based on the Neyman-Fearson theory of hypothesis testing. In many cases a demonstration 
test is specified in terms of total operating time, T, that the system must experience, and a 
maximum allowable number of failures, a. The elements T and a are determined by the following 
preselected factors: 

A. Underlying statistical model in the form of a probability density function of 0. 

B. Minimum acceptable MTBF, 0 , . 

C. Specified MTBF, 0 O , where 0 x < 0 O • 

D. Consumer’s risk, y 2 • The probability of accepting systems with a true MTBF equal to the 
minimum acceptable MTBF, 0 j . The probability of accepting systems with a true MTBF 
less than 0 x will be less than or equal to y 2 . 

E. Producer’s risk, y x . The probability of rejecting systems with a true MTBF equal to the 
specified MTBF, 0 O . The probability of rejecting systems with a true MTBF greater than 
0 O will be less than or equal to y x 2 



1 Numerals in brackets refer to References at end of paper. 

2 AH definitions given here are in agreement with MIL-STD-781B[ 1 1]. 
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Subsequently, minimum T and the associated acceptable number of failures, a, are determined such 
that they satisfy the inequalities 

Pr (Acceptance 1 0 = 0! ) = Pr (f < a 1 0 = 0j ) < y 2 (1) 

and 

Pr (Rejection 1 0 = 0 O ) = Pr (f > a I0 = 0 o )<7i (2) 

where the probabilities are determined by operating on the probability density function and f 
represents the number of failures observed. For example, if the applicable statistical model is 
assumed to be the Poisson distribution, which is quite common in practice, the above inequalities 
would be written 

a f 

Pr(f<a 10 = 0,) = £ e -T / 0 i (T/0 , ) <y 2 (3) 

f=0 

and 

°° f 

Pr (f > a 10 = 0 O ) = £ e- T / 0 o (T/0 O ) /f! < 7l (4) 

f=a+l 

Minimum T and associated a are obtained by simultaneous solution of (3) and (4). Upon concluding 
the test, the customer considers the MTBF requirement to be satisfied and accepts the system if the 
number of failures in T hours is less than or equal to a. Otherwise the system is rejected. 

Wald’s [12] sequential decision approach wherein a third decision, i.e., to continue testing if no 
accept or reject decision can be made, has reduced the expected test time to reach a decision when 
compared to the method just described. Wald’s theory eventually resulted in the well known 
MIL-STD-781 (AGREE) [11] demonstration test plans now imposed on the aerospace industry in 
all serious reliability programs. However, due to increasing reliability requirements and the 
increasing complexity of aerospace hardware, the test time needed to demonstrate reliability has 
reached a point where the cost of such tests has become prohibitive in spite of the efficiencies of 
AGREE test plans. Thus, the interest in Bayesian approach and its associated advantages listed in 
the introductory paragraph. 

BAYES THEOREM 

Let H l5 H 2 . . ., H n represent a mutually exclusive and exhaustive collection of hypotheses. 
Suppose that an event S exists and the conditional probabilities P(S/H i ) are know. Also suppose 
that the probabilities P(H i ) are known. P(H|) is termed the prior probability that Hj is true. Then 
for the discrete case the following conditional relationship is known as Bayes’ Theorem: 



P(H:) (P(S I H:) 

P(H i IS) = ,P(S)*o 
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where 



n 



P(S)= £ P(Hj) P(S I Hj) 
i=l 



P(H^I S) is the posterior probability that is true upon observing a sample statistic, S. In the 
continuous case Bayes* Theorem may be expressed äs 



„»IS, fffl 1 * 5 " 1 » 

f oo 

g(0) h(S I0)d0 



where g(0) is the prior probability density function with the MTBF 0, as the random variable, 
P(0 IS) is the posterior probability density function with 0 as the random variable and h(S|0) 
represents the conditional sampling distribution of S. The theorem, as expressed above, permits 
probabilistic statements regarding 0 upon observing the statistic S. Thus, observed information (test 
data) is combined with the pripr density via the sampling distribution to form the posterior density. 
The conditional sampling distribution in Bayes’ Theorem having the same interpretation as the 
sampling distribution in the classical approach presents no new problem. On the other hand, the 
prior density, represented by the density, g(0), in Bayes* Theorem does require special attention. 

RELIABILITY DEMONSTRATION-BAY ESI AN APPROACH 

Insofar as reliability demonstration is concerned, within the Bayesian context, the consumer is 
interested in being assured that the product received, i.e., the accepted product, has an MTBF no 
less than a minimum acceptable value, 0j. On the other hand, it is in the best interest of the 
producer that the accept-reject decision criteria be structured to assure that the true MTBF of the 
product is indeed less than the specified 0 O if the product is rejected. Here, 0! < 0 O - Stated slightly 
differently, the consumer requires the probability to be 1 — ß* that the true MTBF of the product 
be greater than 0 , given that the product has been accepted by the reliability demonstration f est. 
Concurrently, the producer requires that the probability be 1 — a* that the true MTBF of the 
product be less than the specified value, 0 O , given the product has been rejected by the reliability 
demonstration test. Symbolically this is equivalent to 

Pr (0 > 0j I Acceptance) > 1 - ß* 



and 



Pr (0 < 0 O I Rejection) < 1 - a* 

Note that the preceding inequalities refer to probabilistic statements regarding 0 after an accept or 
reject decision has been made. 
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The foregoing may be reduced to a set of mathematical inequalities in view of Bayes’ Theorem. Let 
the prior distribution of the MTBF be known, agreed upon by the consumer and producer and 
represented by the density g(0). The interpretation of the conditional sampling distribution h(S|0) 
is identical to that of the sampling distribution employed in the classical approach. I.e., S is an 
observed statistic being a function of f failures, T test time, and 0 is the MTBF. Then the 
requirements of the consumer and the producer may be written as follows: 

a 

h [ f , T 1 0 ] d0 



/•tf i 

J 0 



Pr (0 <0j I f < a, T) = 



E 

f=0 



/; 



g(0) 



a 

E 

f=0 



■<ß* 



(5) 



h [f, Tl0] d0 



/ 6 °° 

° g(0) E h[f,T|0] d0 

, 0 f=a+l 

Pr (0 < 0 O I f > a, T) = > 1 — a* (6) 

OO 

g(0) E h[f,Tl0|]d0 

f=a+l 

We let the observed number of failures and time; f and T, respectively, be the observed statistics. 
Thus, we replace h(S|0) with h [f, TI0] where the latter may be represented by a discrete 
distribution such as the Poisson for a practical application. 

The summations in the above expressions provide the conditional probability that 0 < f < a and f > 
a + 1, respectively, as required by the left hand sides. It is clear from examining Eqs 5 and 6 that the 
MTBF, 0 is a random variable where 0 < 0 < «> and a* and ß* are the Bayesian producer’s and 
consumer’s risks, respectively. This is a point of departure from the classical approach wherein the 
risks concern the probability of passing or failing the test itself as implied by Eqs 1 and 2. Note that 
h[f, T I 0] may be a continuous density if necessary, requiring only that the summation signs in Eqs 
5 and 6 may be replaced by integral signs. 

Upon stating the inputs 0 0 , 6 x a* and ß*, simultaneous solution of Eqs 5 and 6 will generate many 
feasible solutions (T, a). Moreover, the minimum T, say T m , and associated a may be determined for 
the given inputs. Thus, the decision rule is to test the product for a total of T m hours; if a or less 
failures are observed the product is accepted, otherwise, it is rejected. Note that this decision rule is 
the same as in the classical case and upon varying the inputs it is possible to develop a table of 
reliability demonstration plans. 

The major differences between the classical and Bayesian approaches are now evident. First, the 
prior distribution concept does not exist in the classical interpretation. In the Bayesian framework 
it is of paramount importance. This implies that if efficiencies are to be realized through Bayesian 
methodology careful attention must be given to the selection of the prior distribution. Second, the 
interpretation of the consumer’s and producer’s risks is fundamentally different. In the classical case 
the consumer wishes to minimize the probability that the reliability demonstration test accepts a 
bad product. The producer wishes to minimize the probability that the reliability demonstration 
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rejects a good product. In the classical approach, given the inputs 7i > 7i > » and the underlying 

model, all of which generate the decision rule, the relative frequency with which the accept or reject 
decision is made reduces to a dependence upon the true, but unknown, value of 0. 

In the Bayesian framework the consumer wishes to have maximum assurance that the product that 
has been accepted and he has received will meet his MTBF requirement for his application; while 
the producer wishes maximum assurance upon product rejection that the MTBF is indeed less than 
the value 0 O . The decision rule is generated by the inputs a*, ß*, 0 lf 0 O , and the conditional 
sampling distribution as in the classical case, however, thereafter, it depends entirely on the prior 
distribution employed. It is evident, therefore, that the classical decision rule is oriented toward the 
value of the MTBF going into the demonstration test. The Bayesian viewpoint is oriented toward 
having both the producer and the consumer be primarily concerned with the probabilistic assertions 
they may make after the reliability test has been conducted and a decision regarding the product 
has been reached. The classical view might be meaningful for the producer, for he is, of course, 
greatly interested in not having a product rejected and, therefore, attempts to enter a reliability test 
with as high an MTBF as possible within scheduling and budget constraints. However, the classical 
approach appears to be less meaningful for the consumer since his primary concern is to have high 
assurance, within scheduling and budget constraints, that the product he receives will meet his 
requirement. 

The third major difference is of an obvious mathematical nature. In the classical case, the true 
MTBF, 0, is an unknown constant whereas in the Bayesian case, 0 is a random variable having 
customary properties including a probability density function describing it. This permits the 
different interpretations discussed earlier and again points out the importance of the prior 
distribution. 

THE PRIOR DISTRIBUTION AND PARAMETER SELECTION 

Application of Bayes’ Theorem to reliability demonstration requires that all objective and subjective 
MTBF information, be integrated to produce a prior density of the MTBF before the demonstration 
test is conducted. Subsequently, to bridge the gap between this information and a convenient 
mathematical form, one of the idealized statistical probability density functions is “fitted” to the 
information. A number of density functions are candidates for such a fit. Selection of a particular 
one can be based on statistical tests for goodness of fit or other more subjective considerations. The 
density function selected to represent the prior density herein is the inverted gamma. 

The inverted gamma prior distribution may be written as 

g(6) = jafy [r(ß)(0 - A)P +1 J | e~ a ^ e - A >> 6 > 0 < 7 > 

where 0, the MTBF, is a randonr variable 

0>O 

a, 0>O 



A >0 
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Note that a and ß are the scale and shape parameters, respectively. Figure 1 illustrates some of the 
family of prior densities that can be represented by the inverted gamma. All densities shown have 
origins at A = 0. If A > 0, the shapes remain the same, then each density would be translated to the 
right by a distance, say A, and it is implied that the MTBF cannot be less than A. 

The inverted gamma was selected for the following reasons: 

A. There is good reason to believe that the distribution of the operational^ true) MTBF of 
electronic systems tends toward positive skewness, [5] and [9]. With this in mind, the 
inverted gamma is particularly suitable since upon properly selecting its parameters it 
conveniently fits positively skewed densities. Also, under certain parameter conditions, 
the inverted gamma approaches symmetry. 

B. The inverted gamma is mathematically tractable within Bayes’ Theorem. 

C. Parameter selection for the inverted gamma is relatively simple. 

PARAMETER SELECTION 

Specifying values for oc and ß , the scale and shape parameters, respectively, of Eq 8 is tantamount to 
selecting a particular probability density function from the inverted gamma family. A procedure for 
specifying a and ß is developed by recognizing that if the value of any two fractile points of an 
inverted gamma density are known, the distribution can be completely specified. If we let the 
predicted MTBF, 0*, be the median, i.e., the 50th fractile, Fq $q, that cuts off the lower 50 percent 
of the prior distribution it is implied that the predictor of Fq believes that the probability is 
one-half that the true MTBF is less than Fq^q. A second fractile point may be determined by 
subjectively selecting an MTBF value such that the odds are 3: 1 that the true MTBF is less than the 
value selected. The latter value, Fq 75 , cuts off the lower 75 percent of the distribution. Since the 
cumulative distribution of the inverse gamma with the origin at zero is 



G{x; a, (3) = f* [ l/r(/3)] oß e-*/ 6 (ließ* 1 dö = £ U/Ü) e -0 '/* ( a/x )> (8) 

J 0 j =0 

it follows that substituting x = Fq ^q and x = Fq 75 in Eq 8 yields 
^^0.50» a ’ 01 = ® -50 



and 



G[F 0 . 75 ;a,/3] =0.75 




FREQUENCY 
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Figur» 1. Variow Inverted Gamm» Demitiw with Zero Origin and Median at 0 « 100 
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For a selected integer value of ß , expression ( 8 ) provides unique solutions Zj and Z 2 where 



z. _ol / F 0.50 



z 2 - «/Fq.75 



( 9 ) 



Furthermore, the following ratio will be unique: 

Q = z,/z 2 

= (a/F 0 5q)/(«/Fo.75) 

= F 0.75^ F 0.50 

Thus, for a selected value of ß, a unique value is obtained for Q. If various ß values are considered, 
then corresponding values for Q can be generated. Moreover, if a table of (0j, Qj) is constructed, it 
would be necessary to merely calculate the ratio Qj to select the corresponding 0 j from the table. 
Once ß is determined, oc is easily calculated from either expression in Eq 9. Here, the first expression 
will be used. Since the subscript is no longer necessary, 

a=ZF 0.50 (10) 

Therefore, the table can be augmented with corresponding values of c*j for given 0j to provide a 
table of (Qj, j3j, «j). The table that has been constructed (Table 1) is (Qj, 0j, Zj) since the 
proportionality factor Z simplifies the implementation. Table 1, developed with the aid of tables by 
Pearson [ 8 ], presents a range of Q considered useful for practical applications. If the origin is 
shifted to A > 0, all preceding arguments are the same except that (F m — A) is substituted for all 
fractile points, F m , and A is selected subjectively. A detailed parameter selection and evaluation 
procedure is provided in a subsequent section. 

PRIOR DISTRIBUTION EVALUATION 

Once a specific density has been selected by virtue of the calculated Q, a useful graph of the 
cumulative distribution may be generated by substituting ß and a where a = ZFq $q into Eq 8 , and 
setting X = Fq 5Q, Fq 75 and other desired values of F m , a graph of the probability density 
function may be constructed. Either or both of these graphs may be examined to judge if the prior 
density selected does, in fact, represent the prior beliefs of the producer and consumer. Thus, a 
method of interpreting a and ß is available. If the graph is acceptable the inequalities represented by 
Eqs 5 and 6 may be implemented directly. If the selected prior distribution is not acceptable, new 
values for Fq $q and/or Fq 75 must be determined which, in turn, will modify Q and, of course, a 
and ß. The new cumulative distribution and the probability density function may then be 
constructed and examined. Eventually the prior distribution will be acceptable and the process will 
terminate. It is of interest to note that as Q becomes large, ß decreases, causing a more positive 
skewed distribution. When Q decreases, ß increases, forcing the prior toward symmetry. More than 
two points can, of course be use'd for the construction of the distribution function. This is discussed 
in detail in [7, 13, and 14], 
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Table 1 

SELECTED RATIOS OF Q FACILITATING THE ESTIMATION OF ß AND 
Z FOR THE INVERSE GAMMA PRIOR DISTRIBUTION 



Q 


ß 


Z 


1.746 


2 


1.679 


1.548 


3 


2.674 


1.449 


4 


3.672 


1.386 


5 


4.671 


1.344 


6 


5.671 


1.312 


7 


6.671 


1.288 


8 


7.670 


1.268 


9 


8.670 


1.252 


10 


9.670 


1.238 


11 


10.669 


1.226 


12 


11.669 


1.215 


13 


12.669 


1.207 


14 


13.669 


1.199 


15 


14.667 
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FUNDAMENTAL FORMS 



Upon setting the prior distribution equal to the inverted gamma with specific values for a and ß and 
defining the conditional sampling distribution as Poisson, inequalities Eqs 5 and 6 become 



a 



f e V[l/(0-A)]0 +1 ^ (l/f!)e-<^ T >/< ö - A )[T/(0-A)] f dö 
J A 



[l/(0-A]0 +1 ^ (l/f!)e-(« +T )/( e -A)[ T /0-A)] f d0 
f=0 



( 11 )- 




oß [ 1/0— A)]^ + ' (l/fOe^^T^^IT/CO-A)]*^© 

f=a+l 



( 12 ) 



If A - 0, i.e., the origin is at zero, performing the integration in Eqs 11 and 12 under the 
appropriate transformations yields 



a V 

£ (R/S) f (V!/f!) £ e- S /° (S/D) J /j! 

f=0 j=0 

<ß* (13) 

£ (R/S) f (V!/f!) 

f=0 



0-1 a V 

(0-0! ^ (l/j!)e“ W W> -(W/S)^ ^ (R/S) f (V!/f!) ^ (1/j!) e — S S* 

j=0 f=0 j=0 

>1-«* (14) 

(0—1)! — (W/S)^ (R/S) f (V!/f!) 
where f=0 

S = ZP + R 
W= ZP 



V= ß+f- 1 
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and the following proportionality factors have been defined and inserted for convenience 

D= 0 i/0 o (discrimination ratio) 

p = e*ie 0 

R = T/0 O 

Z = Ot/d*. 

Thus, for the A = 0 case, given a*, ß*, 0 O and 6 x where 0j <0 O getting Eq 8 represent the prior density 
of 0 with specified a and ß and letting the Poisson distribution represent the conditional sampling 
distribution, Eqs 13 and 14 may be solved for minimum R and the associated a. Upon determining 
values for R and a, the following application and decision rule may be invoked: Test the subject 
system for total of T = R0 O hours. If a or less failures are observed, the system is accepted as having 
fulfilled the pertinent consumer and producer reliability requirements. If the number of failures 
observed is greater than a, the system is rejected. 

By varying the inputs a*, ß*, 0 O , 6 x and selecting specific priors, reliability demonstration test plans 
may be generated by Eqs 13 and 14. Such plans were constructed for practical ranges of the input 
variables and are presented in Table 2. If A > 0, Eqs 13 and 14 can be modified easily and Table 2 
may be employed upon applying a simple adjustment. 

IMPLEMENTING THE PLANS 

The Bayesian risks ol* and ß* are specified by producer and consumer negotiation as are 0 O and 0! . 
Routine reliability predictions provide 0 *, thus, P = 0*/0 O and the discrimination ratio is 
D = 0 1 /0 O - In order to select a specific demonstration plan it remains only to select the specific 
inverted gamma prior density, i.e., a and ß. An iterative routine for selecting a and ß when the origin 
is assumed to be zero (A = 0) follows: 

A. Let F m denote the m ** 1 fractile point on the prior distribution such that Pr(0 < F m ) = m. 
Let Fq 5 q and Fq 75 be two such points. Set the predicted MTBF, 0*, equal to Fq ^q, 
i.e., the median. 

B. Subjectively select an MTBF value such that the odds are 3: 1 that the true MTBF is less 
than the value. Denote the selected value Fq 75 . 

C. Calculate the ratio Q = Fq 75 /FQ $q. 

D. Enter Table 1 at the nearest value of the calculated ratio, Q, and read the values for ß and 
Z directly. Note that a = Z0*. 

E. Substitute ß and a (where a = Z0*) into G[x; a, 0]. Set x = Fq 5 Q, FQ 75 and other 
desired values of F m and graph the resulting cumulative distribution functions (or 
probability density function if desired). If the resulting graph satisfies the subjective 
beliefs of the producer and consumer, go to Item F. If the graphed prior density is 
unsatisfactory select new values for Fq and/or Fq^ and return to Item C. 
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F. Enter Table 2 with a*, 0*, ß , Z, P and D, and read R and a directly. 

The minimum required test time is T m = R0 O and the acceptance number is a failures. 

If it is determined beforehand that it is more meaningful that the origin of the inverted gamma prior 
is at some value A > 0, the preceding iterative routine applies with two modifications. First, all 
references to the F m must be replaced by (F m — A) and 0 l5 0 O , and 0* must be replaced by 
(0! — A), (0 O — A), and (0* — A), respectively. Second, the origin of the graph of the cumulative 
distribution function (Item E) must be labeled A. 

EXAMPLES 

Example 1 

Initial conditions: The prior density is represented by the inverted gamma with A = 0, 0j = 700 
hours, 0 O = 1,000 hours; and- (a*, 0*) = (0.1, 0.1). 

A. Let 0* = Fq 5 q = 1,000 hours. 

B. Let Fq 75 = 1,450 hours. 

C. Q = 1,450/1,000= 1.45. 

D. From Table 1 , (ß, Z) = (4,3.672). 

E. Let the resulting cumulative prior distribution and prior density function, shown in 
Figures 2 and 3, respectively, be acceptable to the consumer and producer. 

F. Enter Table 2 with (a*, 0*, 0, Z, P, D) = (0. 1, 0.1, 4, 3.672, 1.0, 0.7) which yields 
(R, a) = 5.05,7). 

Interpretation: Test the system for R0 O = 5,050 hours. If 7 or less failures are observed accept the 
system, otherwise reject it. 

Example 2 

Initial conditions: Same as in Example 1 except that A = 100 hours, 0j = 770 hours, 0 O = 1,060 
hours and (a*, 0*) = (0.1, 0.1). 

A. Let 0* = Fq 5Q - A = 960 hours. 

B. Let Fq 75 - A = 1 ,300 hours. 

C. Q = 1 ,300/960 = 1 .354. (Note D = 670/960 = 0.7) 

D. From Table 1 , (0, Z) = (6, 5.67 1). 
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E. Let the resulting cumulative prior distribution and prior density function be acceptable to 
the consumer and producer. 

F. Enter Table 2 with (a*, 0*, 0, Z, P, D) = (0. 1, 0. 1, 6, 5.671, 1.0, 0.7) which yields 
(R,a) = (4.68,7). 

Interpretation: Test the system for R0 O = 4,961 hours. If 7 or less failures are observed accept the 
system, otherwise reject it. 

COMPARISONS AND SUMMARY 

A typical example of the relatively smaller test time requirement due to the Bayesian approach as 
opposed to two common classical plans over a range of P = 6*/d 0 for risks of ( 7 j , y 2 ) ~ (0.1, 0.1), 
(a*, 0) = (0.1, 0.1) and a discrimination ratio D = 6 l /6 0 = 0.8 is exhibited in Figure 4. Although the 
plans are not strictly comparable due to the different interpretation of risk, Figure 4 shows the 
respective test time requirements if the inputs were interpreted in the sense necessary to apply each 
method. The single sample plan was generated by Eqs 3 and 4 and provides a single test time 
requirement for all P. The sequential test’s expected test time is from an AGREE type of plan and 
here it is necessary to interpret 0* of P = 0*/0 o as the classical “true” MTBF. The discrimination 
ratio for the later case would be quoted as 1.25:1 in AGREE parlance. The savings due to the 
Bayesian approach are evident from Figure 4. Additional investigations have shown that a decrease 
(increase) in the Bayesian risk or an increase (decrease) in D will effect an increase (decrease) in the 
test time requirement. This relationship parallels classical sensitivity and was expected. 

It is significant to note that in many cases no testing is required by the Bayesian plans, the product 
is judged acceptable, and both producer and consumer requirements are satisfied. This result is due 
to the influence of a “good” prior distribution. That is, the shape of the prior density is such that 
the requirements are fulfilled without additional information from a demonstration test. This 
further points out the importance of prior distribution selection. Some significant work in the 
development of useful prior distributions for the MTBF of electronic hardware has been conducted 
by Feduccia and Klion [5] and Schafer, et al [9]. Also of interest is the work of Lin and Schick 
[7] , and Winkler [ 13] , [ 14] . These efforts tend to indicate that the all important prior distribution 
concept is indeed capable of meaningful, practical interpretation. Thus, the savings in reliability 
demonstration test time requirements typified in this paper are a very real possibility. 
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Table 2 

RELIABILITY DEMONSTRATION TEST PLANS BASED ON THE 
INVERTED GAMMA PRIOR DISTRIBUTION 
(R = UPPER TABLE VALUE, a = LOWER TABLE VALUE) (Continued) 



0 = 6, Z = 5.7 a* = 0.05, ß* = 0.05 



D 


0.8 


0.7 


0.5 


0.8 


P 










1.7 


t 


t 


t 


t 


1.5 


t 


t 


t 


t 


1.2 


22.76 


8.21 


t 


t 




29 


14 






1.1 


25.64 


10.97 


t 


t 




31 


16 






1.0 


26.94 


12.82 


t ' 


t 




31 


17 






0.9 


28.05 


13.86 


1.31 


t 




31 


17 


4 




0.8 


26.67 


14.10 


2.83 


t 




28 


16 


5 




0.7 


23.55 


13.52 


3.66 


t 




23 


14 


5 




0.5 


11.73 


7.61 


2.96 


t 




7 


4 


1 





tNo Testing Required. 



0 = 13, Z= 12.7 


a* = 


0.05,0* = 


0.05 


D 


0.8 


0.7 


0.5 


0.3 


P 










1.7 


t 


t 


t 


t 


1.5 


t 


t 


t 


t 


1.2 


15.25 


t 


t 


+ 




25 








1.1 


27.68 


t 


t 


t 




36 








1.0 


33.70 


11.00 


t 


t 




40 


17 








36.08 


14.83 


t 


t 




40 


19 






0.8 


34.89 


16.09 


t 


t 




36 


18 








27.66 


14.19 


1.68 


t 




25 


13 


2 




0.5 


9.22 


7.28 


3.39 


t 




0 


0 


0 
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Table 2 

RELIABILITY DEMONSTRATION TEST PLANS BASED ON THE 
INVERTED GAMMA PRIOR DISTRIBUTION 
(R = UPPER TABLE VALUE, a = LOWER TABLE VALUE) (Continued) 



0=3, 


Z= 2.7 


a* = 0.1,0 








0=4 


\,Z = 3.7 


a* = 0.1,0 


* = 0.1 




D 


0.8 


0.7 


0.5 






D 


0.8 


0.7 


0.5 


0.3 


P 












P 










1.7 


t 


t 


t 


t 




1.7 


t 


+ 


t 


t 


1.5 


3.54 


t 


t 


t 




1.5 


t 


t 


t 


t 




7 




















1.2 


6.92 


3.51 


t 


t 




1.2 


7.98 


2.85 


t 


t 




9 


6 










11 


6 






1.1 


8.14 


4.07 


t 


t 




1.1 


8.78 


3.77 


t 


t 




10 


6 










11 


6 






1.0 


7.95 


4.56 


t 


t 




1.0 


9.49 


5.05 


t 


t 




9 


6 










11 


7 






0.9 


8.40 


5.00 


0.86 


t 




0.9 


9.39 


5.05 


0.35 


t 




9 


6 


2 








10 


6 


2 




0.8 


8.10 


4.77 


1.27 


t 




0.8 


9.22 


5.59 


1.10 


t 




8 


5 


2 








9 


6 


2 




0.7 


6.99 


4.49 


1.63 


t 




0.7 


8.22 


5.44 


1.60 


t 




6 


4 


2 








7 


5 


2 




0.5 


4.54 


3.79 


1.80 


0.26 




0.5 


4.34 


3.56 


1.51 


t 




2 


2 


1 


0 






1 


1 


0 




tNo Tes 


ting Required 
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Table 2 

RELIABILITY DEMONSTRATION TEST PLANS BASED ON THE 
INVERTED GAMMA PRIOR DISTRIBUTION 
(R = UPPER TABLE VALUE, a = LOWER TABLE VALUE) (Continued) 



0=6, Z = 5.7 


a* = 0.1,0 


* = 0.1 






ß ~ 


13, Z= 12.7 


a* = 


0.1,0* = 0.1 




D 


0.8 


0.7 


0.5 


0.3 




D 


0.8 




0.5 


0.3 


P 












P 










1.7 


t 


t 


t 


t 




1.7 


t 


t 


t 


t 


1.5 


t 


t 


t 


t 




1.5 


t 


t 


t 


t 


1.2 


7.13 


t 


t 


t 




1.2 


t 


+ 


t 


t 




11 




















1.1 


9.85 


2.81 


t 


t 




1.1 


7.19 


t 


t 


t 




13 


6 










13 








1.0 


10.99 


4.68 


t 


t 




1.0 


13.31 


t 


t 


t 




13 


7 










17 








0.9 


11.96 


5.65 


t 


t 




0.9 


15.52 


5.60 


t 


t 




13 


7 










17 


8 






0.8 


11.32 


5.84 


0.55 


t 




0.8 


14.24 


6.83 


t 


t 




11 


6 


2 








13 


7 






0.7 


9.77 


5.90 


1.06 


t 




0.7 


8.60 


4.95 


0.03 


t 




8 


5 


1 








4 


2 


0 




0.5 


4.59 


3.66 


1.8 


t 




0.5 


7.90 


6.12 


2.56 


t 




0 


0 


0 








0 


0 


0 




tNo Testing Required. 
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Table 2 



RELIABILITY DEMONSTRATION TEST PLANS BASED ON THE 
INVERTED GAMMA PRIOR DISTRIBUTION 
(R = UPPER TABLE VALUE, a = LOWER TABLE VALUE) (Continued) 



ß = 


3, Z= 2.7 


a* = 0.2, j3 


* = 0.2 






ß= 


4, Z = 3.7 


a* = 0.2, ß* = 0.2 




D 


0.8 


0.7 


0.5 


0.3 




D 


0.8 


0.7 


0.5 


0.3 


P 












P 










1.7 


t 


t 


t 


t 




1.7 


• 


t 


t 


t 


1.5 


t 


t 


t 


t 




1.5 


t 


t 


f 


t 


1.2 


0.95 


t 


t 


t 




1.2 


0.22 


t 


f 


t 




2 












2 








1.1 


1.44 


0.48 


t 


t 




1.1 


1.74 


f 


t 


t 




2 


2 










3 








1.0 


1.84 


1.09 


t 


t 




1.0 


2.39 


0.84 


t 


t 




2 


2 










3 


2 






0.9 


2.20 


1.04 


t 


t 




0.9 


2.32 


1.47 


t 


t 




2 


1 










2 


2 






0.8 


2.54 


1.37 


t 


t 




0.8 


2.77 


1.44 


t 


t 




2 


1 










2 


1 






0.7 


2.21 


1.13 


0.27 


t 




0.7 


2.52 


1.30 


0.19 


t 




1 


0 


0 








1 


0 


0 




0.5 


2.09 


1.66 


0.81 


t 




0.5 


2.58 


2.03 


0.93 


t 




0 


0 


0 








0 


0 


0 




tNo Testing Required. 


















164 



Table 2 

RELIABILITY DEMONSTRATION TEST PLANS BASED ON THE 
INVERTED GAMMA PRIOR DISTRIBUTION 
(R = UPPER TABLE VALUE, a * LOWER TABLE VALUE) (Continued) 



0=6, Z=5.7 


a* = 0.2,0 


* = 0.2 




D 


0.8 


0.7 


0.5 


0.3 


P 










1.7 


t 


t 


t 


t 


1.5 


t 


t 


t 


t 


1.2 


t 


t 


t 


t 


1.1 


0.99 


t 


4 . 

t 


t 




3 








1.0 


2.25 


t 


f 


t 




3 








0.9 


3.07 


1.27 


t 


t 




3 


2 






0.8 


3.13 


1.52 


t 


t 




2 


1 






0.7 


2.36 


1.57 


t 


t 




0 


0 






0.5 


3.49 


2.70 


1.12 


t 




0 


0 


0 





0= 13, Z= 12.7 


a* = 


0.2, 0* = 0.2 




D 


0.8 


0.7 


0.5 


0.3 


P 










1.7 


t 


t 


t 


t 


1.5 


t 


t 


t 


t 


1.2 


t 


t 


t 


t 


1.1 


t 


t 


t 


t 


1.0 


0.82 


t 


t 


t 




3 








0.9 


3.80 


t 


t 


t 




4 








0.8 


3.28 


1.00 


t 


t 




1 


0 






0.7 


3.85 


2.26 


t 


t 




0 


0 






0.5 


6.39 


4.80 


1.62 


t 




0 


0 


0 





tNo Testing Required. 



165 



REFERENCES 

1. H. Balaban. A Bayesian Approach to Reliability Demonstration. Annals of Reliability and 
Maintainability Conference, 8th Conference, 1969, P 497-504. 

2. A. J. Bonis. Bayesian Reliability Demonstration Plans. Annals of Reliability and 
Maintainability Conference, 5th Conference, 1966, P 861-873. 

3. A. M. Breipohl, R. R. Prairie, and W. J. Zimmer. A Consideration of The Bayesian Approach in 
Reliability Evaluation. I.E.E.E. Transactions in Reliability, Vol. 14, October 1965, P 107-1 13. 

4. C. H. DeWitt, R. H. Myers, K. L. Wong, and E. G. N. Yem. Reliability Prediction and 
Demonstration for Airborne Electronics. RADC Contract Number F30602-67-C-0221, April, 
1968. 

5. A. J. Feduccia and J. Klion. How Accurate are Reliability Predictions? Proceedings, Annual 
Symposium on Reliability, 1968, P 280-297. 

6. C. W. Hamilton. Bayesian Procedures and Reliability Information. Second Annual Aerospace 
Reliability and Maintainability Conference, 1963, P 278-283. 

7. C. Y. Lin and G. J. Schick. On-Line Console-Aided Assessment of Prior Distributions. Ninth 
Reliability and Maintainability Conference, Detroit, Michigan, 22-20 July 1970. 

8. Karl Pearson. Tables of the Incomplete Gamma Function. University Press, Cambridge, 
England, 1946. 

9. R. E. Schaefer, J. Collins, M. L. Luden, et al. Bayesian Reliability Demonstrations: Phase 
1 — Data for the A Priori Distribution. RADC-TR-69-389, Final Technical Report, February, 
1970. 

10. R. E. Schaefer and N. D. Singpurwalla. A Sequential Bayes Procedure for Reliability 
Demonstration. Annals of Reliability and Maintainability Conference, 8th Conference, 1969, 
P 507-514. 

11. Reliability Tests: Exponential Distribution. U. S. Department of Defense, MIL-STD-781B, 
Washington, D. C., Government Printing Office, 1967. 

12. A. Wald. Sequential Analysis. John Wiley and Sons, Inc., New York, 1947. 

13. R. L. Winkler. The Assessment of Prior Distributions in Bayesian Analysis. Journal of the 
American Statistical Association, Vol. 62, No. 319, September 1967, P 776-800. 

14. R. L. Winkler. The Qualification of Judgement: Some Methodological Studies. Journal of the 
American Statistical Association, Vol. 62, No. 320, December, 1967, P 1 105-1 120. 




166 



Lösungsvorschläge für stochastische Zielprogramme 
von W. Rödder, Aachen 



Obersicht : Für ein stochastisches Zielprogramm der Form 

Gx + V - y* 

s . d . Ax ‘- } b 
x * 0 

werden deterministische Ersatzprobleme vorgeschlagen, in- 
dem 

1. ) die Wahrscheinlichkeit der Zielerreichung maximiert bzw. 

2. ) die Wahrscheinlichkeit der Nicht-Zielerreichung minimiert 

wird . 

Es werden Familien von Verteilungen für V angegeben, bei dene: 
diese Aufgaben zu konvexen bzw. konkaven Optimierungsproblemen 
führen. Ferner lassen sich einfache hinreichende Bedingungen 
für die Existenz von optimalen Lösungen aufstellen. 

Einführung in die deterministische Zielprogrammierung 

Untersuchungen über die Kriterien, an denen Erttscheidungsträge] 
in Unternehmungen ihre Entscheidungen messen, führten in den 
letzten Jahren zu der Hypothese, daß ein Modell, in dem die 
Entscheidungen mittels einer eindimensionalen "Zielgröße" ge- 
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ordnet werden (Gewinn, Erlös u.a.), nicht in jedem Fall 
dem realen Sachverhalt gerecht wird. 

Vielmehr zeigt z.B. Heinen ( HL ) , da ß der Unternehmer 
sich simultan von mehreren Zielen leiten läßt. In einem 
mathematischen Modell spiegelt sich diese Tatsache in einer 
vektorwertigen Zielfunktion (F) wider. 

Sie gibt den funktionellen Zusammenhang zwischen Ent- 
scheidungen x und Zielgrößen F(x) an. 

Wir wollen i.f. unterstellen, daß sich der Entscheidungs- 
träger hinsichtlich der einzelnen Zielkomponenten 
F^(x) : i = 1, ..., r Anspruchsniveaus gesetzt hat; 

diese Anspruchsniveaus seien bezeichnet mit 
f % 3t. x 

» • • • » y r > = : y • 

Sind die möglichen Entscheidungen gewissen Restriktionen 
g^(x) ( * J 0: k = 1, ..., m unterworfen, läßt sich der 

Imperativ 

’’Triff Deine zulässige Entscheidung so, daß das Anspruchs- 
niveau y möglichst gut erreicht wird” formal schreiben als: 

F (x) * y* 

s.d. g^(x) 0 k = 1, ..., m (1) 

Der Ausdruck F (x) ^ y kann konkretisiert werden in eine 
Abweichungsfunktion oder allgemeiner eine Präf erenzordnung , 
mittels der der Entscheidende die Zielerreichungsgrade F(x) 

3t 

bzgl. der Anspruchsniveaus y ordnete 

Hierzu finden sich einige Vorschläge bei Charnes und 

Cooper TOC! S. 215 ff, Ijiri [ Ul S. 55 und Rödder L^O] 

S. 50 ff. 

Unterstellt man in (1), daß F(x) und g^(x) linear sind, er- 
hält man ein lineares Zielprogramm der Form 
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Gx i y* 

s . d . Ax f -? b 
x i 0 



( 2 ) 



(Hierbei sei: x € JR n , y t |R r , bf p m , G eine 
r x n - Matrix und A eine m x n - Matrix). 

Als einfaches Beispiel für (2) mag folgender Sachverhalt 
dienen (vgl. CIJJ s. 45 ) : 

In einer Firma arbeiten 2 Gruppen mit verschiedenen 
Produktionsraten ; 

Gruppe 1 stelle pro Stunde 1 Einheit 
Gruppe 2 " " "0,6 Einheiten her. 

Pro Tag sollen insgesamt - falls möglich - mindestens 
12 Einheiten hergestellt werden, und jede Gruppe sollte 
möglichst höchstens 8 Stunden arbeiten. 

In das Modell (2) gefaßt ergibt das: 




s.d. x ± * 0 i = 1,2 
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Stochastische Lineare Zielprogrammierung 

In der Literatur zur stochastischen linearen Programmierung 
wurde hinreichend diskutiert, daß in einem Modell der Form 

Max cx 

s.d. Ax b (3) 

x - 0 



die Komponenten des Vektors (c, A, b) gewöhnlich nur 
Schätzungen der sich in der Zukunft realisierenden Koeffi- 
zienten sind. Liegen die Informationen über diese Koeffi- 
zienten in Form von Wahrscheinlichkeitsverteilungen vor, 
beschreibt ein stoehastisches L.P. -Modell 

Max cx 

s.d. Ax b 

x - 0 

(wobei c, A, b) ein Zufallsvektor ist) 
besser den realen Sachverhalt als (3). 

Diese Überlegungen verlaufen völlig analog für lineare 
Zielprogramme, die bei stochastischem (G, A, b) die 
folgende Form haben: 

Gx * y* 

s.d. Äx -<*> E (1 

X äs 0 

Einschränkend soll i.f. angenommen werden, daß (4) für den 
Fall zu lösen ist, daß die Restriktionen deterministiBfch 
sind und die Zielformulierung 
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^ _ 

Gx + V - y ist, wobei G fest und V ein r-dimensionaler 

Zufallsvektor ist. 

Insgesamt ergibt sich also das Problem 
Gx / V - y* 

s . d . Ax = b ( 4a ) 

x 0 



1 • Maximierung der Wahrscheinlichkeit der Zielerreichung 

Die Aufgabe (4a) bedarf einer Interpretation, da der 

Imperativ "Wähle di-ine zulässige Entscheidung so, daß 

_ £ 

der Zufallsvektor Gx + V möglichst > y ist” 
zunächst einmal sinnlos ist. 

Wir nehmen i.f. an, der Entscheidungsträger habe eine 
konkrete Nutzenkonzeption, und zwar ordne er solchen 
Ergebnissen von 

Entscheidung x und 

Realisation der "Umwelt" V° 

den Nutzen 1 zu, falls Gx + V° ^ y* 

und 0 sonst , 

Die Maximierung des erwarteten Nutzens ergibt 
dann das deterministische Entscheidungsproblem 

Max P { Gx t V * y* ] 
s.d. Ax b 



x 



0 
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Eine solche Nutzenkonzeption ist z.B. sinnvoll, wenn eine 
Nicht-Zielerreichung in auch nur einer Komponenten 
existenzgefährdend für den Entscheidungsträger ist. 

Bei unternehmerischen Entscheidungen könnten folgende 

3t . ... 

Beispiele genannt werden: Die sind Liquiditäts- 

forderungen oder sie sind in einem Produktionsmodell die 
Auftragsmengen eines stark dominierenden Kunden, bei 
deren Nichterreichen ein Abbruch der wirtschaftlichen 
Beziehungen droht! 

Sind F^(w) die Verteilungsfunktionen der stetig ver- 
teilten, unabhängigen Zufallsvariablen V., gilt: 

r 1 

P [Gx + V i y* j = P {(Gx + V). * y ± * j = 

i=l 

r 

IT P [ V. ^ y.* - (Gx).. J 
i = l 

Unter derv genannten Voraussetzungen ergibt sich für (5) 

• 3t 

also mit - (Gx)^ und nach Logarithmierung der 

Zielfunktion : 
r 

Max 1 ln (1 - F. (z. )) 

* — - l l 

i = l 

s.d. Xi = y* - (Gx) i i e [l, r 1 

Ax ) b 



r 

II (1 - F i (y i *-(Gx) i » 
i=l 



X 



0 



( 6 ) 
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Wir geben i.f. hinreichende Bedingungen für die Konkavität 
der Zielfunktion und die Existenz einer optimalen Lösung 
von ( 6 ) an . 



Satz 1 : Sind die (z^ stetig differenzierbar und 

sind F^ 1 ( z^ ) = f^z^ Polya - Typ 2 Dichten , 

ist (6) eine konkave Maximierungsaufgabe. 



Bew. 



dz. 



ln (1 - F, 



(z^) = 



f . ( Zi ) 
nrnTT 



Nun gilt nach CbMp] , S. 378: 



sind f i (z i ) Poly.a-Typ 2 Dichten , so ist 
f. (z.) 

„ 1 _ r monoton nichtfallend; daraus folgt 
1_F i (2 i ) 

die Behauptung 



Für eine weite Klasse von Verteilungen sind die Be- 
dingungen dieses Satzes erfüllt. So haben z.B. folgende 
Verteilungen Polya-Typ 2 Dichten. 

Die Gamma verte il ung , die Weibullverteilung (vgl. 

LBMPJ) , S. 378)mit den Dichten. 

x ( x t f _1 <*? <- xt) / r (co a>o, t * o) 
Xoct * -1 (-Xt*> (X>0, t * o) 

für«*- - 1 . 

die Normalverteilung, die F-Verteilung und die ^-Verteilung 
(vgl. LKARl , S. 124 f). 

Wir wenden uns jetzt der Existenz einer optimalen Lösnng 
von (6) zu: Da die Zielfunktion nach oben beschränkt ist, 




muß nur garantiert werden, daß es überhaupt ein 
zulässiges x gibt, für das alle ln (1 - f^(z^))> - 00 «ind 
nauer gilt folgender 

Satz 2 : Ist K : = ( x : Ax £ ^ ) b, x * 0 ] kompakt und 

F\ ( z^ ) stetig für it fl, ..., rj , so ist 

hinreichend für die Existenz einer optimalen 
Lösung von (6), daß es ein x £ K gibt, so daß 
für jedes i € {, 1 , * r] die Wahrscheinlich- 
keit, das Ziel zu erreichen, positiv ist« 

Bew. : Da ln (1 - ( z ^)) * 0 (also nach oben 

beschränkt), ist hinreichend für die Existenz 
einer optimalen Lösung von (6), daß 

1 - Fi (z.) >0 i € ll, ...» r] bzw. 

r 

TT (1 - W > 0 

i=l 

Hat z.B. V auf den ganzen R T eine positive Dichte, ist 
die Bedingung des Satzes erfüllt. 

Beispiel einer speziellen Verteilung der : 

Ist V. auf t - -X ,*•) ( X. > 0) it { 1, . . . , r ] 

X ^ V ^ 

exponentiell verteilt, ergibt sich für (6) folgendes 
Optimierungsproblem : 
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-Min H y± 
i=l 

s.d. y i = (l + X ^ (y^* - (Gx ) ^ ^ + 

Ax * b 
x * 0 

Dies ist ein lineares deterministisches Optimierungs- 
problem, das z.B. mit dem Simplexalgorithmus gelöst 
werden kann. 



2 . Minimierung der Ziel-Nichterreichung 

Die Aufgabe (4a) kann eine andere Interpretation als 
unter 1. erfahren. 

So unterstellen wir dem Entscheidungsträger jetzt 
fofcg4nde Nutzenkonzeption: 

Einer Entscheidung x und der Realisation der "Umwelt” G° 
ordne der Entscheidungsträger den Nutzen 0 zu, falls 
G°x<y*, und sonst den Nutzen 1. 

Das besagt, der Entscheidende ist zufrieden, sobald auch 
nur in einer Komponenten das Ziel erreicht wird. Diese 
Bewertung der Zielvariablen ist dann sinnvoll, wenn die 
Ziele miteinander völlig substituierbar sind. 

Bei den gleichen Voraussetzungen und Bezeichnungen wie 
unter 1 ergibt sich folgendes Minimierungsproblem: 



1 Wobei für 1 Q |R 



1 + = max (0, 1) 
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r 

Hin £1 in (z^) 

i=l 

s.d. z. = y .* - (Gx) . 

i i l 

(7) 

Ax ^ b 
x = 0 



Analog zu 1. läßt sich wieder eine hinreichende Bedingung 
für die Existenz einer optimalen Lösung zu (7) angeben: 

1 

Satz 3 : Sind die F i ( z stetig, K kompakt, so ist 

hinreichend für die Existenz einer optimalen 
Lösung von (7), daß es kein x € X gibt, so 
daß das Ziel y ^ für irgendein i für jede Reali- 
sation von V. erreicht wird. 

(Der Beweis erübrigt sich). 

Das Problem (7) ist ein konkaves Minimierungsproblem, falls 
die V. i€ [l , ..., r] gamma- ,weibull- oder normalverteilt 
sind. Bei den ersten beiden Verteilungen folgt diese Tagsache 

aus f 1 (z.) * 0 ( ln F. (z. ) * 0), bei der 

1 dz. 1 1 

1 

Normalverteilung zeigt man die Behauptung mittels einer 
Diskussion der 2. Ableitung 



1 



siehe Satz 2 
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Zusammenfassung 

Die Arbeit verfolgt zwei Ziele. Zunächst wird unter Verwen- 
dung des Dynamischen Programmierens eine neue Herleitung des 
Theil* sehen Satzes über dynamische Sicherheitsäquivalente ge- 
geben. Sodann wird ein Satz hergeleitet, der die Existenz dy- 
namischer Sicherheitsäquivalente auch für nichtquadratische 
Kriterien nachweist. 

1 . Einleitung 

Der Begriff dynamisches Sicherheitsäquivalent wurde bereits 
1957 von H. Theil T7] geprägt. Grobgesprochen besagt er, daß 
in gewissen stochastisch dynamischen Optimierungsproblemen 
die auftretenden stochastischen Prozesse durch bestimmte de- 
terministische Größen, die sogenannten dynamischen Sicher- 
heitsäquivalente, ersetzt werden können. Oder präziser, Theil 
konnte den folgenden Satz beweisen [8j : "Die optimale Politik 
eines durch eine lineare Zustandstransformationsgleichung und 
ein quadratisches Kriterium beschriebenen (dynamisch) stocha- 
stischen Optimierungsproblems ändert sich nicht, wenn man die 
auftretenden Zufalls variablen durch gewisse bedingte Erwar- 
tungswerte ersetzt." Diese bedingten Erwartungswerte nennt 
man dynamische Sicherheitsäquivalente . 

Vir wollen uns diesen Satz etwas veranschaulichen. Gegeben 
sei das Rückkopplungs-Regelungsproblem der Fig. 1.1 




Fig. 1.1 
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Dann besagt der Theil'sche Satz folgendes: Für den Fall, daß 
die Regelstrecke eine lineare Transformation beschreibt und 
das Kostenkriterium quadratisch in der Zustands- und Regel- 
größe ist, kann die stochastische Störung des Systems durch 
dynamische Sicherheitsäquivalente ersetzt werden, d.h. durch 
bedingte Erwartungswerte der Störvariablen unter der Bedin- 
gung der bis zum gegenwärtigen Zeitpunkt erhältlichen Infor- 
mation. 

Noch deutlicher wird die Bedeutung des Theil' sehen Theorems 
für den Unternehmensforscher, wenn man den obigen Regelungs- 
vorgang als einen Produktionsglättungsvorgang interpretiert. 
In diesem Falle spezialisiert sich Fig. 1.1 zu Fig. 1.2 



stoch. 

Nachfrage 




-* lagern 



Produktion 



produzieren ^ 



Fig. 1.2 



Lagerbestand v 

| 7 



Sollbestand 



Hier besagt der Satz von Theil, daß unter der per definitio- 
nem (backlog-Fall) gegebenen linearen Bilanzgleichung und 
bei Vorliegen eines quadratischen Kostenzusammenhangs der 
stochastische Nachfrageprozeß für die Optimierung durch eine Folge 

bedingter Erwartungen ersetzt werden kann. Diese bedingten Er- 
wartungswerte sind aber nichts anderes als Erwartungsprognosen 
der Nachfrage. 

Damit wird aber die Bedeutung des Nachweises der Existenz 
dynamischer Sicherheitsäquivalente (zumindest innerhalb der 
Lagerhaltungstheorie) klar ersichtlich. Zeigt doch der Satz 
von Theil, daß bisher nur sichergestellt ist, daß das in 
weiten Bereichen der Lagerhaltungstheorie meist unreflektiert 
verwandte Verfahren der Nachfrageprognostizierung lediglich 
im Falle quadratischer Kostenabhängigkeiten nicht suboptimal 
ist. Quadratische Kriterien treten aber in Problemen der 
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Unternehmensforschung nicht allzu häufig auf. Es wäre daher 
von größtem Interesse, den Theil' sehen Satz in der Weise zu 
erweitern, daß auch nichtquadratische Kostenabhängigkeiten 
auf dynamische Sicherheitsäquivalente führen. 

Wir wollen in dieser Arbeit zwei Probleme behandeln. In einem 
ersten Teil soll zunächst eine neue Ableitung des Satzes von 
Theil mit Hilfe des Dynamischen Programmierers gegeben werden. 
Theil selbst hat zur Ableitung seines Satzes ein quadratisches 
Variationsverfahren benutzt, das die Inversion einer i.a. 
hochdimensionalen Matrix erforderlich macht. Es ist daher 
nicht uninteressant, diesen wichtigen Satz nochmals unter 
Verwendung des Dynamischen Programmierers abzuleiten. Der 
zweite Teil der Arbeit ist dann einer Erweiterung des Satzes 
auf nichtquadratische Kriterien gewidmet. 

2. Ein neuer Beweis des Satzes von Theil über die Existenz 
dynamischer • Sicherheitsäquivalente . 

A) Das Modell 

Wir betrachten das durch die 3 Be stimmungs großen a) System, 
b) Störung und c) Kostenkriterium festgelegte stochastisch 
dynamische Optimierungsproblem. 

a) System 
Es seien 

x k : Zustand z.Z. k (k = 1 , ... ,N) 
x Q : reellwertig vorgegebener Anfangs zustand 
u^ : Steuerung z.Z. k (k = 0, ... ,N-1) 
r k : Störung z.Z. k (k = 1, ... ,N) 

Diese Variablen seien verknüpft durch die lineare Zustands- 
transformationsgleichung 

x k = x k-i + u k-i + r k c 2 - 1 ) 

Ferner hänge die Steuerung von sämtlichen realisierten (der 
Messung zugänglichen) Zuständen ab 

u k-1 = Vk-I^k-I» ’ x 1’ x o } (2.2) 

Die Gin. (-2.1) und (2.2) bilden das ’’System" . Es sei noch- 
mals zum besseren Verständnis des durch (2.2) beschriebenen 
Rückkopplungsvorganges und insbesondere im Hinblick auf die 




Fig. 2.1 



Verzögerung ■ 



Y k-1 ^ x k-1 ’ * * • * X 1 ’ x o } 



b) Störprozeß 

Der Störprozeß {r^.} (k = 1, ... ,N) sei eine extern vorgege- 
bene Folge von Zufallsvariablen, für die die bedingten Erwar- 
tungsire rte 



E (r k | r j-i’ x o3 

existieren mögen. 



(k > j > i > 0) 



Si) Kostenkriterium 

Das Kostenkriterium sei gegeben durch den quadratischen Aus- 



druck 



= + 1 



(2.3) 



Hierbei bedeutet f eine positive Konstante und o( einen Dis- 
kontfaktor. 

Das Optimierungsproblem kann nun folgendermaßen formuliert 
werden: Es ist eine Folge u°, ... ,u^_^ derart zu finden, 
daß der Erwartungswert der Kosten minimiert werde. Unsere Auf- 
gabe wird nun sein, zu zeigen, daß die Folge {u£} nicht von 
dem Störprozeß fr^.} , sondern lediglich von gewissen bedingten 
Erwartungswerten der r^ abhängt. Zunächst noch eine Bemer- 
kung. 

Bemerkung 

Das soeben definierte Modell ist auch unter den linear-qua- 
dratischen Modellen von sehr einfacher Struktur. Es ist aber 
einsichtig, daß das im folgenden benutzte Verfahren der 
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Dynamischen Programmierung ohne Schwierigkeiten auch auf 
multivariate Modelle mit etwas reichhaltigerer linearer 
(stabiler) Zustands transformation und komplizierterer qua- 
dratischer Kostenfunktion angewandt werden kann.*) 

Daß wir gerade dieses Modell gewählt haben, liegt an seiner ein- 
fachen Interpretierbarkeit als Produktionsglättungsmodell. 

Für ein solches Modell wäre x^. zu identifizieren mit dem 
(positiven bzw. negativen) Lagerbestand am Ende der Periode 
(k-1,k], r^ stellt (nach Multiplikation mit -1) die Nach- 
frage während dieser Periode dar (kumuliert in k) , und u^.^ 
beschreibt den Produktionsentscheid am Ende der Periode 
(k-2,k-1] . Die Zustancftransf ormationsgleichung (2.1) ist die 
bekannte Bilanzgleichung der Lagerhaltungstheorie, und die 
Kostenfunktion (2.3) beschreibt diskontierte quadratische 
Lager- und Produktions(abweichungs) kosten während der N Ent- 
scheidungsperioden. 

B) Lösung des Optimierungsproblems 

Zur Lösung des Optimierungsproblems mit Hilfe der Dynami- 
schen Programmierung geht man aus von dem folgenden System 
von Funktionalgleichungen 



f k+F x 2 = min E { x k+1 + ? 2u k + * f k+2 ( ' x k+l1 r k’ r k-1’ ••'’ r 1’ x ol 



(k = 0,1, ...,N-2) 



mit der Anfangsgleichung (k = N-1) 

^N^ X N-1^ = m ^* n ^{ X N + ? U N-1 1 r N-1 ’ * * * ,r 1 ,X ol 

U N-1 

(Bemerkung: Zur Vereinfachung der Notierung haben wir in der 
Wertfunktion f^ + ^ die Abhängigkeit von der ’’Vorgeschichte" 
l r k ,r k-l’ ...’rpxj nicht gekennzeichnet.) 

Damit reduziert sich das obige Optimierungsproblem auf die 
Lösung der Funkt ionalgle ichungen . Diese Lösung existiert und 
ist eindeutig. Man erhält nach einer Reihe von Rechenschrit- 
ten, die im Prinzip in [4] zu finden sind, der folgenden 
Lösungs algor ithmus : 

*) Eine solche Erweiterung wurde als Teil einer Dissertation am 
Lehrstuhl für Operations Research der FU Berlin bereits durch- 
geführt . 
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a) Optimale Politik 

u °_k = S ;T-k (x K-k + ; (k = 



? « s i;_ 



k+1 



J N-k 



? - 1 



; (k = 2,...,N) 



o k := r (N-k+1 1 N-k, . . . ) + — ; (k = 2,...,N) 

2W N-k+1 



‘TJ-k 



^ f T N-k+1 l r N-k’ * * * ,r 1 ,x d 



T, t 



N-k+1 



2 <? clS I . i ^ k+1 Q N _ 



k+1 



W- 



N-k+1 # 



: = i - ^as. 



N-k+1 



r( N-k+1 |N k,...) E { r ]\T_}£+^ I * * * ,r 1 ,x o } 

und speziell 

r(k|0, ...) := E { r k |x Q } 

Die Anfangswerte des Algorithmus sind gegeben durch: 
1 



D N-1 ~ „ «2 

1 + Q 



Q n .1 s= r(N| N-1 , . . . ) 



b) Wertfunktionen %«k+1 ^ x N-k^ 

^N-k+I^N-k^ = " ?^* S N-k x N-k + T N-k x N-k + %-! 



(2.4) 

( 2 . 5 ) 

( 2 . 6 ) 

( 2 . 7 ) 

( 2 . 8 ) 

( 2 . 9 ) 

( 2 . 10 ) 
( 2 . 11 ) 

( 2 . 12 ) 

( 2 . 13 ) 



(2.14) 
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N-k 



: = ot pN-k+1 S N-k^N-k + W N-k+1 CT ’rr. ( ' N-k+1 * N_k ’ 



• •) + 



+ 


T r (N-k + 1) + %_ k ] 


(2.15) 


II 

1 


E f R N-k+1 1 r R-k’ * * * ,r 1 ,x 0 } 


(2.16) 


T r (N-k+1) 


: = E { T N-k+1 r N-k+1 1 r N-k » ' * * ’ r 1 ’ x o J 


(2.17) 


CT r:r (N-k+1 N-k,...) := E ( r N-k+1 1 r N-k’ * * * ,r 1 ,X of 


(2.18) 


mit dem Anfangswert 

%-1 : = Sj^CNlN-l,...) + a rr (N|N-1,...) 


(2.19) 



Betrachten wir -zunächst den Algorithmus für die optimale Po- 
litik, Setzt man sukzessive die Gin. (2.5) ••• (2.15) in Gl. 
(2.4) ein, so zeigt sich, daß die optimalen Entscheidungen 
u° lediglich von bedingten Erwartungswerten der Form (2.10) 
abiSagen. Das heißt aber, daß es für die Bestimmung der opti- 
malen Politik iu°} lediglich auf die durch (2.10) definierten 
dynamischen Sicherheitsäquivalente ankommt und nicht auf die 
volle stochastische Struktur des Störprozesses. M.a.W. man 
bekäme dieselbe optimale Politik, wenn man von vornherein 
den stochastischen Störprozeß durch seine dynamischen Sicher- 
heitsäquivalente ersetzt hätte. Damit ist der Theil’sche Satz 
bewiesen. 

3 . Ein Satz über die Existenz dyn amischer Sicherheitsäqui- 
valente bei nichtquadratischen Kr iterien 

Wir wenden uns nun der Erweiterung des Theil’ sehen Satzes 
auf nichtquadratische Kriterien zu. In jüngster Zeit wurden 
in verschiedenen Richtungen Anstrengungen unternommen, zu 
weiteren Aussagen über dynamische Sicherheitsaquivalente zu 
gelangen. Zu erwähnen sind etwa die Arbeiten von Ziemba [9], 
Midier [3]" und Malinvaud [2]. Ziemba und Midier weisen die 
Existenz dynamischer Sicherheitsäquivalente nach, falls die 
Kostenfunktion linear vom Zustand abhängt und weitere weniger 
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einschneidende Voraussetzungen erfüllt sind. Die lineare Ab- 
hängigkeit von den Zustandsvariablen (bei nichtbeschränktem 
Zustandsraum) stellt aber für die meisten Probleme der Unter- 
nehmensforschung eine nicht vertretbare Einschränkung dar. 
Malinvaud dagegen führt für ein relativ allgemeines nicht- 
lineares und nichtquadratisches Modell eine Sensitivitäts- 
analyse durch. Er untersucht die Frage, wie die optimale Poli- 
tik eines deterministischen Modells sich ändert, wenn gewisse 
Parameter des Modells in definierter Weise stochastisch wer- 
den. Liese Analyse setzt allerdings u.a. eine zweifache Dif- 
ferenzierbarkeit der Kostenfunktion nach den Zustands varia- 
blen und das Fehlen statischer Nebenbedingungen voraus, bei- 
des Voraussetzungen, die in den Anwendungen der Unternehmens- 
forschung nur selten gegeben sein dürften. 

Wir wollen nun einen weiteren Satz über dynamische Sicher- 
heitsäquivalente ableiten. Auch dieser Satz enthält noch eine 
stark einschränkende Voraussetzung, ist aber zumindest den 
Problemen der Lagerhaltungstheorie, auf die -wir unser beson- 
deres Augenmerk richten, besser angepaßt als die zuvor er- 
wähnten Untersuchungen. Drei wesentliche Voraussetzungen bzgl. 
der Modellkonstituenten a) , b) und c) werden wir fordern müs- 
sen 

VI: Das System sei linear. 

V2: Der Störprozeß sei ein Gaußprozeß. 

V3: Das Modell befinde sich in einem "eingescffiingenen Zustand” 
(steady state) und lasse ein Kriterium zu, das mittlere 
Kosten pro Periode beschreibe. 

Voraussetzung VI ist zweifellos am einschneidensten. Sie be- 
sagt, daß nicht nur (wie im Satz von Theil) die Zustandstrans- 
formationsgleichung (2.1) als linear vorausgesetzt wird, son- 
dern auch die Funktion in Gl. (2.2). Wie wir wissen, ist 
z.B. für Lagerhaltungsprobleme vom backlog Typ die Zustands- 
transformationsgleichung ex definitione linear; kritisch ist 
daher lediglich die geforderte lineare Abhängigkeit der Poli- 
tik von den Zustandsvariablen. Wir werden darauf zum Schluß 
noch einmal zu sprechen kommen. 

Die Voraussetzung V2 ist in vielen Fällen nicht kritisch. 

Wie Box und Jenkins [1] nachgewiesen haben, lassen sich viele 
ökonomische Zeitreihen als Realisierungen (i.a. höher kor- 
relierter und häufig stationärer) Gaußprozesse interpretieren. 




las 

V3 schließlich besagt, daß sich das Kostenkriterium als Er- 
wartungswert einer vom Zustand und der Steuerung abhängigen 
Funktion f(x,u) darstellen läßt 

P := E{f(x,u)j ( 3 . 1 ) 

wobei F also nicht mehr von der Periodennummer abhängt. Ein 
solcher "eingeschwungener Zustand" wird in vielen Fällen schon 
nach wenigen Entscheidungsperioden erreicht. (Im Falle einer 

quadratischen Durchschnittswertkostenfunktion K = y 2 u^ 

kann man zeigen, daß schon nach* wenigen Perioden der asymptoti- 

Wert lim E{K} = Efx^+f^u 2 } = O ' 2 + o 2 * o ' 2 erreicht wird, 

P xx j uu 

wobei und <T UU respektive die Varianzen von Zustand und 

Steuerung sind C 5 J ) • An f sind, wie wir sogleich sehen wer- 
den, noch weitere allerdings nur schwache Forderungen zu 
stellen. 

Der Beweis für die Existenz dynamischer Sicherheitsäquiva- 
lente kann nun sehr rasch geführt werden. 

Aufgrund von V2 und VI folgt, daß Zustand und Steuerung eben- 
falls normalverteilt sind, so daß 

P := Ejf(x,u)} = F( 0 ^, CT^, M*, ( 3 . 2 ) 

lediglich von den ersten und zweiten (zentrierten) Momenten 
der Zustands- und Steuervariablen abhängt. Ohne wesentliche 
Beschränkung der Allgemeinheit wollen wir der Einfachkeit 
halber a n nehmen, daß F nicht von den Mittelwerten und 
und der Kovarianz abhängt. Häufig wird F von der 

Kovarianz sowieso nicht abhängen, da Zustand und Steuerung 
in der Kostenfunktion separiert auf treten, d.h. 

f (x,u) = L(x) + P(u) . ( 3 . 3 ) 

Damit wird aus (3.2) 

P = p (° r L’°'uE • (3.4) 

Eine notwendige Bedingung für die Existenz eines Minimums 
von F ist nun das Verschwinden der 1. Variation von F (nach 

der die lineare Politik bestimmenden Gewichts folge) 
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Hierbei wird vorausgesetzt, daß die Ableitungen 
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existieren. Ferner wollen wir fordern, daß 

8F 3F 

— >0 und — >0 -und beschränkt ist, so daß der Quotient 



3 G 



xx 



3 <5 



uu 



3F 



0" 



da 



uu 






(3.6) 
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positiv und beschränkt ist. 

Man beachte, daß die Forderungen nach Existenz der Ableitungen 
von F nach den Varianzen weit weniger einschränkend ist als 
Malinvauds Forderung nach zweimaliger Differenzierbarkeit 
von f nach dem Zustand. Unstetigkeiten in f sind dabei durch- 
aus zugelassen. Dadurch wird es möglich, für f auch solche 
Funktionen zu verwenden, die geeignet sind, mögliche Beschrän- 
kungen im Zustands- und/oder Ent s che idungs raum zu beschreiben. 
Die Positivität der Ableitungen dürfte in den allermeisten 
Fällen aus ökonomischen Gründen sichergestellt sein, 
o 

Mit 0 wird nun aus (3*5) 



& G 2 



e 2 6 ot 



= 0 



(3.7) 



Dieser Ausdruck stellt ein bemerkenswertes Ergebnis dar. 

Er ist nämlich nichts anderes als die erste Variation eines 
quadratischen Kriteriums mit dem "reduzierten Funktional“ ([6], Kap. 7) 



+ 



0 2 CT 2 

uu 



(3.8) 



Damit ist die Existenz dynamischer Sicherheitsäquivalente 
aber fast schon bewiesen. Denn für ein durch (2.1), V2 und 
(3.8) definiertes Optimierungsproblem weist man leicht nach 
[6] , daß die optimale Politik gegeben ist durch 
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u k = v ?rX x k + £ z^(e)r(k+1+klk. . .)] , (3.9) 

k* = 0 

wo y Konstanten sind, 

2(6 ) f - W — < 1 

1 2B d 

ist und r(k+1+k| k. . . ) die durch (2.10) definierten dynami- 
schen Sicherheitsäquivalente darstel^n. D.h. auch hei nicht- 
quadratischem Kriterium existieren unter den gegebenen Vor- 
aussetzungen dynamische Sicherheitsäquivalente, wobei man 
o 

sich 0^ lediglich über (3.6) zu verschaffen hat. Treten da- 

2 2 

bei mehrere 0 -Werte auf, so hat man dasjenige 0 >0 zu ver- 
wenden, das auf die niedrigsten Kosten führt. 

Wir fassen das Ergebnis noch einmal zusammen in dem folgenden 
Satz : 

Unter den Voraussetzungen VI bis V3 und der zusätzlichen 

2 2 

Bedingung, daß die Ableitungen von F nach und flf 

p XX LlUI 

positiv und (bzgl. 0*^) beschränkt sind, existieren dyna- 
mische Sicherheitsäquivalente. 

Hierbei ist die einzige Voraussetzung, die wirklich kritisch 
ist, die Forderung VI nach Linearität der Politik. 

Für quadratische Kriterienist allerdings, wie man (2.4) ent- 
nimmt, diese Linearitätsvoraussetzung erfüllt. Es ist daher 
anzunehmen, daß für eine große Klasse nichtquadratischer Kri- 
terien eine lineare Politik zumindest eine gute Approximation 
darstellt. Somit erwiese sich z.B. das häufig unreflektiert 
angewandte Verfahren der Nachfrageprognostizierung als in vie- 
len Fällen durchaus gangbar. 
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Stochastische Programmierungsmodelle 
als Vektormaximumprobleme 
von W. Dürr, Regensburg 
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In der Mathematischen Programmierung sind die Möglichkeiten der Optimie- 
rung genau einer Zielfunktion die hauptsächlichen Forschungsziele . In den 
Anwendungsgebieten der Mathematischen Programmierung, z.B. in der Be- 
triebswirtschaftslehre, ist fast ausschließlich das Gewinnmaximierungs- 
prinzip die Grundlage optimaler Entscheidungen. In neuerer Zeit hat das 
Problem der mehrfachen Zielsetzung (Zielkonflikt) oder Vektormaximumpro- 
blem zunehmend Beachtung gefunden. 

Definition 1: 

Es sei X C $ n ; z k (x) seien K Zielfunktionen, also z k : $ n ? 

(k * 1, . . . ,K) . Die Aufgabe 

z l(x) 
t Z K <*>/ 

2 

heißt Vektormaximumproblem. 

Die Anführungszeichen bringen zum Ausdruck, daß nicht nur eine Ziel- 
funktion zu maximieren ist. Es ergibt sich die Frage, welche Lösungen 
x c X als " optimal” bezeichnet werden können. Dies sind sicherlich nur 
solche x, die gemäß der folgenden Definition (funktional) effizient 
sind. 

Definition 2: 

% e X heißt genau dann (funktional) effizient bezüglich X und 
z 1 (x),...,z K (x), wenn kein x* e X mit der Eigenschaft 

* Vgl. z.B. Kuhn-Tucker [1951], Karlin [1959], Bod [196.3], Gutenberg 
[1966], Geoffrion [1967], Johnson [1968]. 

2 

Entsprechend wird ein Vektorminimumproblem definiert. 



c" £ 



(VMP) "max” }z(x) | x e X 1 



mit z (x) * 
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und 



z k (x') > z k (?) (k = 1.....K) 

z k , ( x* ) > z k ,(|) für mindestens ein k' e {1, 



,K} 



existiert. 

Bei einem Entscheidungsproblem können nicht effiziente Lösungen vor vorn- 
herein außer Betracht bleiben. Der Entscheidungsträger wird sicherlich nur 
ein effizientes x e X wählen, allerdings bleibt die Präge, welches effi- 
ziente x er vorziehen wird. Er kann ein effizientes x wählen, das Lösung 
eines " Kompromißprogramms " ist. Es bieten sich verschiedene Möglichkei- 
ten für derartige Kompromißprogramme an 1 , z.B. 

K 

(KP1) max { l t, z,(x) | x e X} 

k=l K K ’ 

mit t^ > 0 (k * 1, . . . ,K) . 



(KP2) max (z^x) | x e X} 

unter den Nebenbedingungen 

z 2 (x) > d 2 
z K (x) > d R . 

Die Frage, inwieweit Lösungen von (KP1) und (KP2) effizient sind, soll 
hier nicht diskutiert werden , desgleichen nicht, inwieweit die Menge 
aller effizienten Lösungen, die als "vollständige Lösung" von (VMP) be- 
zeichnet sei, durch (KP1) und (KP2) bestimmt werden kann. Es sollen im 
folgenden nur gewisse Modelle der stochastischen Programmierung als Pro- 
gramme der Art (KP1) und (KP2) erkannt werden, was dann die Interpretation 
als Komp romißpro gramme gewisser Vektormaximumprobleme erlaubt. 

2. Stochastische Programmierung 
Gegeben sei das lineare Programm 



1 



Vgl. Geoffrion [1965], Dinkelbach [1971], Dürr [1971]. 




